r - 查找导致低 p 值的条目

Question

在 R 中有 2 个向量

u <- c(109, 77, 57, 158, 60, 63, 42, 20, 139, 15, 64, 18)
v <- c(734, 645, 1001, 1117, 1071, 687, 162, 84, 626, 64, 218, 79)

我想测试 H:u并且v是独立的，所以我运行卡方检验：

chisq.test( as.data.frame( rbind(u,v) )  )

并得到一个非常低的 p 值，这意味着我可以拒绝 H，这意味着它u不是v独立的。

但是当我输入

chisq.test(u,v)

我得到一个 0.23 的 p 值，这意味着我可以接受 H。我应该选择这两个测试中的哪一个？

此外，我想在这些向量中找到导致这个低 p 值的条目。任何想法如何做到这一点？

score 0 · Accepted Answer

检验统计量使用标准化残差平方和。您可以查看这些值以了解特定值的重要性

m = chisq.test(u, v)

residuals(m)
m$stdres

1 回答 1