R语言笔记 独立性 检验

> mytable1 <- xtabs(~ Treatment+Improved,data=Arthritis)

> mytable2 <- xtabs(~ Improved+Sex,data=Arthritis)

> mytable1

Improved

Treatment None Some Marked

Placebo   29    7      7

Treated   13    7     21

> mytable2

Sexa

Improved Female Male

None       25   17

Some       12    2

Marked     22    6

> chisq.test(mytable1)

Pearson‘s Chi-squared test

data:  mytable1

X-squared = 13.055, df = 2, p-value = 0.001463

> chisq.test(mytable2)

Pearson‘s Chi-squared test

data:  mytable2

X-squared = 4.8407, df = 2, p-value = 0.08889

Warning message:

In chisq.test(mytable2) : Chi-squared近似算法有可能不准

在结果中,患者接受的治疗和改善的水平看上去存在着某种关系(p < 0.01)。而患者性别

和改善情况之间却不存在关系(p > 0.05) 。这里的p值表示从总体中抽取的样本行变量与列变

量是相互独立的概率。由于的概率值很小,所以你拒绝了治疗类型和治疗结果相互独立的原假

设。由于的概率不够小,故没有足够的理由说明治疗结果和性别之间是不独立的。

产生警告信息的原因是,表中的6个单元格之一(男性 - 一定程度上的改善)有一个小于5的值,

这可能会使卡方近似无效。

时间: 2024-08-08 17:29:50

R语言笔记 独立性 检验的相关文章

R语言笔记

循环: repeat { your code if( ) { break } } while (test_expression) { statement } for ( i in v) { print(i) } 函数: function_name <- function(arg_1, arg_2, ...) { Function body } 函数的返回值是要评估的函数体中的最后一个表达式 pair与map: list() #Create a hash object h = hash() #ad

R语言笔记5--读数据

1.读文本文件数据 (1)先设置工作目录,把文本文件放于该目录下 备注:在记事本里写完数据后,按一下回车,负责在R语言中出现错误 (2)读剪贴板 文本或EXCEL的数据均可通过剪贴板操作 (3)读excel文件数据 方法1:先把excel另存为空格分隔符的prn文件格式再读     方法2:安装RODBC包,再通过ODBC读

R语言之t检验

在R中,t检验是通过t.test()命令实现的,通过设置该命令中的选项可以完成不同的t检验.首先介绍该命令的常用选项: 一.单样本t检验 通过设定mu选项来进行,如> t.test(data2,mu=5) One Sample t-test data:  data2t = 0.25482, df = 15, p-value = 0.8023alternative hypothesis: true mean is not equal to 595 percent confidence interv

R语言笔记 提醒程序员注意的一些事项

经验丰富的程序员通常会发现R语言的某些方面不太寻常.以下是这门语言中你需要了解 的一些特性. 对象名称中的句点(.)没有特殊意义.但美元符号($)却有着和其他语言中的句点 类似的含义,即指定一个对象中的某些部分.例如,A$x是指数据框A中的变量x. R不提供多行注释或块注释功能.你必须以#作为多行注释每行的开始.出于调试目的, 你也可以把想让解释器忽略的代码放到语句if(FALSE){... }中.将FALSE改为TRUE 即允许这块代码执行. 将一个值赋给某个向量.矩阵.数组或列表中一个不存在

R语言笔记 R中设置图形参数--函数par()详解

R有着非常强大的绘图功能,我们可以利用简单的几行代码绘制出各种图形来,但是有时候默认的图形设置没法满足我们的需要,甚至会碰到各种各样的小问题:如坐标轴或者标题出界了,或者图例说明的大小或者位置遮挡住了图形,甚至有时候默认的颜色也不能满足我们的需求.如何进行调整呢?这就用到了"强大"的函数par().我们可以通过设定函数par()的各个参数来调整我们的图形,这篇博文就是对函数par()的各个参数进行一下总结. 函数par()的使用格式如下: par(..., no.readonly =

R语言笔记之语法篇

R语言基础之语法 分支结构Ifelse和 Ifelse函数 循环结构 for循环 while循环 repeat循环 breaknext 函数和包 1函数 1创建 2查看函数代码 3导入R文件 2包的安装和加载 1包的安装 2包的加载 向量化计算 1普通的加减乘除 2apply家族 1 applyxMARGIN2 2 lapply列表函数名函数其他参数 3 sapply列表函数名函数其他参数 4 mapply列表函数名函数其他参数 5 tapply向量因子因子列表函数函数名 3split 排序 数

R语言笔记之数据篇

R语言杂七杂八 与R语言有关的应用工具 探索性数据分析 统计推断 回归分析 机器学习-分类问题 R与Rstudio的获取与安装 包package一种扩展R基本功能的机制集成了众多函数 获取包 导入包libraryname 获取帮助 R语言特点 R语言语法基础之数据篇 R语言中的数据 R语言支持的数据类型 基本数据结构 一维数据类型 向量 vocter 1创建 2提取子集 因子 factor 1创建 2提取子集 二维数据类型 矩阵 matrix 1生成矩阵 1matrix方法 2 修改dim属性来

R语言笔记4--图形分析

  1.对x1进行直方图分析,绘制直方图hist()       2.探索各科成绩的关联关系,散点图绘制函数plot()            3.列联表分析,列联函数table(),柱状图绘制函数barplot() 4.饼图,饼图绘制函数pie()     5.箱尾图 箱子的上下横线为样本的25%和75%分位数 箱子中间的横线为样本的中位数 上下延伸的直线为尾线,尾线的尽头为最高值和最低值 异常值            6.箱线图      7.相图 每个观测单位的数值表示为一个图形 每个图的

R语言笔记 计算描述性统计量的函数

其中包括Hmisc.pastecs和psych.由于这些包并未包括在基础安装中. > library(Hmisc) > describe(mtcars[vars]) mtcars[vars] 3  Variables      32  Observations ----------------------------------------------------------------------------------------------------------------------