Cochran q测试RVAideMemoire包

时间:2017-02-08 17:41:52

标签: r

我在" 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矩阵中的所有内容组合在一起,所以我假设我现在组织它的方式是正确的吗?

0 个答案:

没有答案