如何在 R 中获得 Kruskal-Wallis 检验的精确 p 值？

Question

如何在 R 中获得 Kruskal-Wallis（例如 3 组）检验的准确 p 值？

数据示例：

df <- data.frame(
    dv = c(0.80, 0.83, 1.89, 1.04, 1.45, 1.38, 1.91, 1.64, 0.73, 1.46,
           1.15, 0.88, 0.90, 0.74, 1.21),
    group = factor(rep(c("A", "B", "C"), c(5, 5, 5))))

我使用kruskal_test函数尝试了coin包

kruskal_test(dv ~ group, data = df,distribution= "exact")

虽然产生了错误：

Error in .local(object, ...) : ‘object’ is not a two-sample problem

如果我将"exact" 更改为"approximate"，它会运行，但它不是确切的分布...

有什么想法吗？

Answer 1

您收到错误的原因是您只能准确计算两个样本问题的分布。

来自help("kruskal_test")：

...分布可以通过蒙特卡洛重采样来近似，或者通过将分布分别设置为“近似”或“精确”来精确计算单变量双样本问题。