4

我想知道 R 中是否存在执行以下操作的包/函数:“比例同质性卡方检验的事后成对比较”(或等效项),此处描述: http:// epm.sagepub.com/cgi/content/abstract/53/4/951

我的情况只是在 2 x X 矩阵上进行 chi 测试。我发现了一个差异,但我想知道哪些列对差异“负责”。

谢谢, 塔尔

4

3 回答 3

3

“卡方检验”通常生成为单个单元格偏差与“预期”的平方和 = 行和列总和的乘积除以总和。因此,可以将单个单元格对总和的贡献与 1 df 的卡方临界值进行比较。修改 chisq.test() 代码以返回单元格卡方是一项相当简单的任务。我刚刚补充说:

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

到最后的结构调用。他们不会得到 print()-ed,但您可以将结果分配给一个对象并使用:

 obj$cell.chisq
于 2010-03-20T15:15:12.863 回答
1

请参阅fifer包了解功能 chisq.post.hoc()

于 2015-01-16T21:12:01.947 回答
0

不确定它是否适用于您的问题,但我阅读了一篇文章,该文章使用以下内容执行成对 chi-sqr 测试

rcompanion::pairwiseNominalIndependence()
于 2017-10-08T03:17:26.340 回答