После значительного теста Крускала я хотел бы выполнить несколько тестов по окончании c с помощью тестов Вилкоксона в соответствии с запланированными сравнениями. Я не хочу проводить парные тесты Уилкоксона и контрольную группу. Раньше я делал это с ANOVA с функцией glht, включающей матрицу контраста «K»: glht(lm1, linfct=mcp(x=K)), но я не знаю, как это сделать с непараметрическими c тестами.
В kruskal.test(formula = y~x, data = data) x включает 4 фактора (уровни: A, B, C и D). Я хотел бы выполнить несколько пост-хо c тестов Уилкоксона в соответствии с запланированными сравнениями: A против B + C, B + C против D и A против D.
Я пытался использовать wilcox_test функция в пакете rstatix, но она не работает.
wilcox_test(data = data, formula = data$y~data$x, p.adjust.method ="bonferroni", comparisons=list(c(c("B", "C"), "D"), c(c("B", "C"), "A"), c("D","A"))))
R дает мне скорректированные p-значения для сравнений: B против C, A против B и A против D.
Не могли бы вы помочь мне с этой функцией? Или вы знаете другой способ выполнить то, что я хотел бы сделать?
Спасибо