如何使用R语言进行卡方检验
Posted
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了如何使用R语言进行卡方检验相关的知识,希望对你有一定的参考价值。
参考技术A我们想知道喝牛奶对感冒发病率有没有影响,以下为数据统计的四格表:
从上文中我们通过卡方计算公式手动计算出P在0.1到0.5,说明牛奶对感冒发病率没有影响。
下面我们通过R计算一下具体p直:
x<-matrix(c(43,28,96,84),ncol=2,nrow=2)
chisq.test(x)$p.value
此外出来2X2联表,还可以计算一维表(自由度为个数),例如下面:
判断5种品牌啤酒的爱好者有无显著差异:
P值越大,支持原假设的证据就越强,给定显著性水平α(取0.05), 当P值小于α时,就拒绝原假设。
卡方检验结果怎么看呀,
卡方检验试用条件1.随机样本数据;
2.卡方检验的理论频数不能太小.
两个独立样本比较可以分以下3种情况:
1.所有的理论数T≥5并且总样本量n≥40,用Pearson卡方进行检验.
2.如果理论数T<5但T≥1,并且n≥40,用连续性校正的卡方进行检验.
3.如果有理论数T<1或n<40,则用Fisher’s检验.
上述是适用于四格表.
R×C表卡方检验应用条件:
1.R×C表中理论数小于5的格子不能超过1/5;
2.不能有小于1的理论数.
若不符合R×C表的卡方检验.可以通过增加样本数、列合并来实现.
统计专业研究生工作室原创,请勿复杂粘贴 参考技术A 一个是看卡方值,另一个是看卡方值对应的p值,也就是sig值,尤其是sig,如果sig0.05的话则无充分理由认为实际与期望有差异或者是两变量相关。
以上是关于如何使用R语言进行卡方检验的主要内容,如果未能解决你的问题,请参考以下文章
R语言构建仿真列联表并进行卡方检验(chisq.test):检验两个分类变量是否独立输出期望的列联表
R语言使用gmodels包的CrossTable函数计算获取二维列联表CrossTable函数可以设置输出比率(行比率列比率单元格比率)指定小数位的个数执行卡方检验Fisher精确检验等
R语言使用pwr包的pwr.chisq.test函数对卡方检验进行效用分析(power analysis)在已知效应量(effect size)显著性水平效用值的情况下计算需要的样本量
R语言使用coin包应用于分类变量独立性问题的置换检验(permutation tests)使用普通卡方检验chisq.test函数和置换近似卡方检验chisq.test函数检验分类变量的独立性