如何得到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

如果我更改"approximate"的"exact"，它会运行，但它不是确切的发行版...

有什么想法吗？

Answer 1

你得到错误的原因是因为你只能精确地计算两样本问题的分布。

来自help("kruskal_test")

对于单变量两样本问题，

...the分布可以通过蒙特卡罗重新采样来近似，或者通过将分布分别设置为“近似”或“精确”来精确计算。