如何测试所有样品组合的显着差异？

Question

我正在分析R中9个图的数据集，并希望将它们相互比较。我的响应变量不是正常分布的。

现在我的问题：

要测试的是9 + 8 + 7 + ... + 2 + 1 = 45对组合。 R可以自动执行此操作吗？如果有，怎么样？ 我的愿望输出将是一个盒须图，其中我的9个图在x轴上，y轴上的响应变量和图上方的小写字母表示显着差异。

提前致谢！

Answer 1

这应该让你开始：

#some data
x <- rlnorm(100, mean=1:4)
DF <- data.frame(x=x, g=c("a", "b", "c", "d"), stringsAsFactors=FALSE)

#pairwise Mann-Whitney-U-test
pairwise.wilcox.test(DF$x, DF$g, p.adjust.method = "bonferroni")

#   Pairwise comparisons using Wilcoxon rank sum test 
#
#data:  DF$x and DF$g 
#
#  a       b       c     
#b 0.0016  -       -     
#c 6.3e-09 0.0020  -     
#d 1.9e-13 2.0e-08 0.1823
#
#P value adjustment method: bonferroni