全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
2631 2
2016-12-15
同一题,计算1:1配对资料的卡方值时候R语言程序及结果
> y<-c(1287,61,7,2)
> cc<-matrix(y,nrow = 2,ncol = 2)
> cc
     [,1] [,2]
[1,] 1287    7
[2,]   61    2
> mcnemar.test(cc)

        McNemar's Chi-squared test with continuity
        correction

data:  cc
McNemar's chi-squared = 41.309, df = 1, p-value =
1.3e-10

SAS程序及结果
data ex3;
do r=1 to 2;
do c=1 to 2;
input f@@;
output;end;end;
cards;
1287 7 61 2
;
proc freq;
tables r*c/chisq expected nopercent
nocol agree;
weight f;
run;



                                                     McNemar 检验
                                                 ---------------------
                                                 统计量 (S)    42.8824
                                                 自由度              1
                                                 Pr > S         <.0001

手工计算与SAS结果一致,求老司机解答,是否我R程序出错了?

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2016-12-15 23:05:59
可能问题太low了,木有人回答,在我翻找资料后也有所发现,这篇帖子给了我一个很好的启发:
https://bbs.pinggu.org/forum.php?mod=viewthread&tid=4195825&from^^uid=286968
于是我重新看了遍 mcnemar.test的参数信息,将程序改为
mcnemar.test(cc,correct = FALSE)
于是结果就和SAS一样了
McNemar's Chi-squared test

data:  cc
McNemar's chi-squared = 42.882, df = 1, p-value =5.813e-11

感叹一下,还是R牛啊!!!
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2020-6-20 22:40:59
你的SAS程序没有指定频数,少个weight语句
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群