我在" RVAideMemoire"中遇到了Cochran q测试的一些问题。封装
我的数据如下:
response <- c(0,1,1,0,0,1,0,1,1,1,1,1,0,0,1,1,0,1,0,1,1,0,0,1,0,1,1,0,0,1)
fact <- gl(3,1,30,labels=LETTERS[1:3])
block <- gl(10,3,labels=letters[1:10])
cochran.qtest(response~fact|block)
不幸的是,我总是在cochran.qtest(a~b | c)中得到错误&#34;错误: 每个等级的'b&#39; b&#39;在每个街区&#34;
当我尝试在软件包解释中提供的示例时,运行正常:
table(fact, block)
block
fact a b c d e f g h i j
A 1 1 1 1 1 1 1 1 1 1
B 1 1 1 1 1 1 1 1 1 1
C 1 1 1 1 1 1 1 1 1 1
table(b,c)
c
b 1 2 3
0 90 90 90
1 90 90 90
我尝试使用as.numeric,as.factor更改我的参数,但是错误仍然存在,这让我很沮丧。
我真的很感谢你的帮助!
编辑: 根据评论我检查了以下内容:
{{1}}
这让我很困惑。据我所知,Cochran q是McNemar测试的延伸。 McNemar将2x2矩阵中的所有内容组合在一起,所以我假设我现在组织它的方式是正确的吗?