我想知道 R 中是否存在执行以下操作的包/函数:“比例同质性卡方检验的事后成对比较”(或等效项),此处描述: http:// epm.sagepub.com/cgi/content/abstract/53/4/951
我的情况只是在 2 x X 矩阵上进行 chi 测试。我发现了一个差异,但我想知道哪些列对差异“负责”。
谢谢, 塔尔
我想知道 R 中是否存在执行以下操作的包/函数:“比例同质性卡方检验的事后成对比较”(或等效项),此处描述: http:// epm.sagepub.com/cgi/content/abstract/53/4/951
我的情况只是在 2 x X 矩阵上进行 chi 测试。我发现了一个差异,但我想知道哪些列对差异“负责”。
谢谢, 塔尔
“卡方检验”通常生成为单个单元格偏差与“预期”的平方和 = 行和列总和的乘积除以总和。因此,可以将单个单元格对总和的贡献与 1 df 的卡方临界值进行比较。修改 chisq.test() 代码以返回单元格卡方是一项相当简单的任务。我刚刚补充说:
cell.chisq = (x - E)^2/E,
到最后的结构调用。他们不会得到 print()-ed,但您可以将结果分配给一个对象并使用:
obj$cell.chisq
请参阅fifer包了解功能 chisq.post.hoc()
不确定它是否适用于您的问题,但我阅读了一篇文章,该文章使用以下内容执行成对 chi-sqr 测试
rcompanion::pairwiseNominalIndependence()