如何用Kruskal-Wallis检验比较R中的不同组

Question

我有如下所示的数据：

​

​

我想用Kruskal试验比较A，B，C，D组化合物1的浓度.例如，化合物1，A组和B组，比较它们的浓度。应对同一化合物内的不同群体进行比较。组比较应如下所示。

​

​

我尝试了下面的代码

# Performing Kruskal-Wallis test
kruskal.test(Group~Concentration, data = df) )

但是它使用所有的组A数据(包括compound1，compound2，.)。我想比较一下一个化合物，如何修改代码？

Answer 1

由于您没有使用dput或类似的方法提供示例数据，所以我用示例数据组成了一个数据框架。

您可以使用函数kruskal.test在R中运行Kruskal Test

set.seed(123)
df <- data.frame(group = sample(x = LETTERS[1:4], size = 20, replace = T),
                 compound = rep(1:4, 5),
                 concentration = rnorm(20, mean = 0.4, sd = 0.2))


kruskal.test(concentration ~ group, data = df)

结果是

> kruskal.test(concentration ~ group, data = df)

    Kruskal-Wallis rank sum test

data:  concentration by group
Kruskal-Wallis chi-squared = 0.5648, df = 3, p-value = 0.9044

此外，请查看?kruskal.test中的函数文档。