Cochran q测试记忆包

我对“RVAideMemoire”软件包中的Cochran q测试有一些问题

我的数据如下：

library(RVAideMemoire)
a <- c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 1L, 0L, 0L, 
  0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
  1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L)
b <- c("0", "1", "0", "1", "1", "0", "0", "1", "1", "0", "1", "1", 
  "0", "1", "0", "0", "1", "1", "0", "0", "1", "1", "0", "0", "1", 
  "1", "0", "1", "0", "0", "1", "1", "0", "0", "1", "1", "0", "0", 
  "1", "1", "0", "0", "1", "0", "0", "1", "1", "0", "0", "1", "0", 
  "0", "1", "1", "0", "1", "1", "0", "0", "1", "1", "0", "0", "1", 
  "1", "0", "0", "1", "1", "0", "1", "0", "0", "1", "1", "0", "0", 
  "1", "1", "0", "0", "1", "1", "0", "1", "0", "0", "1", "1", "0", 
  "0", "1", "1", "0", "0", "1", "1", "0", "1", "0", "0", "1", "1", 
  "0", "1", "1", "0", "0", "1", "1", "0", "1", "0", "0", "1", "1", 
  "0", "0", "1", "1", "0", "0", "1", "1", "0", "1", "0", "1", "0", 
  "1", "0", "1", "0", "1", "0", "1", "0", "0", "1", "0", "1", "0", 
  "1", "1", "0", "1", "0", "1", "0", "0", "1", "0", "1", "0", "1", 
  "1", "0", "1", "0", "0", "1", "0", "0", "1", "0", "1", "0", "1", 
  "0", "1", "0", "1", "0", "1", "0", "1", "0", "1", "0", "1", "0", 
  "0", "1", "0", "1", "0", "1", "1", "0", "1", "0", "1", "0", "0", 
  "1", "0", "1", "0", "1", "1", "0", "1", "0", "1", "0", "0", "1", 
  "0", "1", "0", "1", "1", "0", "1", "0", "1", "0", "1", "0", "1", 
  "0", "1", "0", "1", "0", "1", "0", "1", "0", "0", "1", "0", "1", 
  "0", "1", "1", "0", "1", "0", "1", "0", "0", "1", "0", "1", "0", 
  "1", "1", "0", "1", "0", "1", "0", "0", "1", "0", "0", "1", "1", 
  "1", "0", "1", "0", "1", "0", "1", "1", "0", "1", "0", "1", "0", 
  "1", "0", "1", "0", "1", "0", "1", "0", "1", "0", "0", "1", "0", 
  "1", "0", "1", "1", "0", "1", "0", "1", "0", "0", "1", "0", "0", 
  "1", "1", "0", "1", "0", "1", "0", "0", "1", "0", "1", "0", "1", 
  "1", "0", "1", "0", "1", "0", "1", "0", "1", "0", "1", "0", "1", 
  "0", "1", "0", "1", "0", "0", "1", "0", "1", "0", "1", "1", "0", 
  "1", "1", "0", "0", "1", "0", "0", "1", "0", "1", "0", "1", "0", 
  "1", "1", "0", "1", "0", "1", "0", "1", "0", "1", "0", "1", "0", 
  "1", "0", "1", "0", "1", "0", "0", "1", "0", "1", "0", "1", "1", 
  "0", "1", "0", "1", "0", "0", "1", "0", "1", "0", "1", "1", "0", 
  "1", "0", "1", "0", "0", "1", "0", "1", "0", "1", "1", "0", "1", 
  "0", "1", "0", "1", "0", "1", "1", "0", "1", "0", "1", "0", "1", 
  "0", "1", "0", "1", "0", "0", "1", "0", "1", "0", "1", "0", "1", 
  "0", "1", "0", "1", "1", "0", "1", "0", "1", "0", "0", "1", "0", 
  "1", "0", "1", "1", "0", "1", "0", "1", "0", "0", "1", "0", "1", 
  "0", "1", "1", "0", "1", "0", "1", "0", "0", "1", "0", "1", "0", 
  "1", "0", "1", "0", "1", "0", "1", "1", "0", "1", "0", "1", "0", 
  "1", "0", "1", "0", "1", "1", "0", "1", "0", "1", "0", "0", "1", 
  "0", "1", "0", "1", "1", "0", "1", "0", "1", "0", "0", "1", "0", 
  "1", "0", "1", "0", "1", "0", "1", "0", "1", "1", "0", "1", "0", 
  "1", "0", "0", "0", "1", "0", "1", "0", "1", "1", "0", "1", "0", 
  "1", "0", "0", "1", "0", "1", "0", "1")
c <- c("1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
  "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
  "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
  "1", "1", "1", "1", "2", "2", "2", "2", "2", "2", "2", "2", "2", 
  "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", 
  "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", 
  "2", "2", "2", "2", "2", "2", "2", "3", "3", "3", "3", "3", "3", 
  "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", 
  "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", 
  "3", "3", "3", "3", "3", "3", "3", "3", "3", "1", "1", "1", "1", 
  "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
  "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
  "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
  "1", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", 
  "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", 
  "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", 
  "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "3", "3", "3", 
  "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", 
  "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", 
  "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", 
  "3", "3", "3", "3", "3", "3", "1", "1", "3", "1", "1", "1", "1", 
  "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
  "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
  "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
  "1", "1", "1", "1", "1", "1", "2", "2", "2", "2", "2", "2", "2", 
  "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", 
  "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", 
  "2", "2", "2", "2", "2", "2", "2", "3", "3", "3", "3", "3", "3", 
  "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", 
  "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", 
  "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", 
  "3", "3", "3", "1", "2", "2", "2", "2", "2", "2", "2", "2", "3", 
  "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", 
  "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", 
  "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", "3", 
  "3", "3", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", 
  "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", 
  "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", "2", 
  "2", "2", "2", "2", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
  "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
  "1", "1", "2", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
  "1", "1", "1", "1", "1", "1", "1", "1")

cochran.qtest(a~b|c)

---

这让我很困惑。据我所知，Cochran q是McNemar测试的扩展。McNemar将所有内容组合在一个2x2矩阵中，因此我认为我现在组织的方式是正确的？

也许您需要在块级别上汇总您的响应。比较

表（事实，块）

和

表（b，c）

的结果，然后思考错误消息的措辞。恐怕我现在更困惑了。我想要现在的样子：一张6x2的桌子。。。我不认为将所有值拆分为单独的列有什么意义。我认为这个测试是McNemar测试的扩展，它也将所有的值组合成一个2x2矩阵？

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)

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