在R中是否存在卡方检验的配对PostHoc比较？

Question

我想知道在R中是否有一个包/函数来执行下面描述的："Post Hoc Pair-Wise Comparisons for the X-Square Test of the X-Square Test of Proportions“(或它的等价物)：http://epm.sagepub.com/cgi/content/abstract/53/4/951

我的情况是在一个2乘以X的矩阵上做一个卡氏测试。我发现了一个差异，但我想知道哪个列是造成这种差异的“责任”。

谢谢，塔尔

Answer 1

“卡方检验”通常是单个单元格偏离“预期”的平方和=行和列和除以总和的乘积的总和。因此，可以将单个细胞对总和的贡献与具有1D.F.的卡方临界值进行比较。修改chisq.test()代码以返回单元格X平方是一项相当简单的任务。我刚加了一句：

cell.chisq = (x - E)^2/E,

到结尾处的结构调用。它们不会得到print()-ed，但您可以将结果赋给一个对象并使用：

 obj$cell.chisq

Answer 2

请参阅chisq.post.hoc()函数的fifer包

Answer 3

不确定它是否适用于您的问题，但我阅读了an article，它使用以下代码执行成对的chi-sqr测试

rcompanion::pairwiseNominalIndependence()