Как получить точное значение p теста Краскела-Уоллиса (например, с 3 группами) в R?
Пример данных:
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))))
Я попробовал пакет coin
, используя функцию kruskal_test
kruskal_test(dv ~ group, data = df,distribution= "exact")
Хотя выдает ошибку:
Error in .local(object, ...) : ‘object’ is not a two-sample problem
Если я поменяю "exact"
на "approximate"
, он запустится, но это не точное распределение...
Есть предположения?
combn(levels(df$group), 2, FUN = function(x) kruskal_test(dv ~ group, data = subset(df, group %in% x), distribution = 'exact'), simplify = FALSE)
- person akrun   schedule 23.05.2020stats::kruskal.test(dv ~ group, data = df)
? - person duckmayr   schedule 23.05.2020kruskal.test
иks.test
. Я тоже добился успеха с вашимkruskal.test
подходом. - person Ian Campbell   schedule 23.05.2020exact
... - person Sinval   schedule 23.05.2020