如何比较R中的箱形图？

Question

我有三组不同的数据集。 A，B和C. 每组有3列：

A： a，b和c

B：

d，e和f

C：

g，h和i

所以，

在这里，我尝试将我的数据作为示例进行演示：

A:

a <- seq(from=-1, to=4, by=0.5); b <- seq(from=-1, to=8, by=0.5); c <- seq(from=-2, to=3, by=0.5)

B:

d <- seq(from=-4, to=1, by=0.5); e <- seq(from=-3, to=2, by=0.5); f <- seq(from=-1, to=5, by=0.5)

C:

g <- seq(from=-1, to=1, by=0.5); h <- seq(from=-2, to=3.5, by=0.5); i <- seq(from=-4, to=0.5, by=0.5)

我为他们制作了箱图。

boxplot(a,b,c, col=c(4:6))
boxplot(d,e,f, col=c(4:6))
boxplot(g,h,i, col=c(4:6))

答

B：

C：

我的问题是，如何将所有三个箱图中的第1组（'深蓝'a，d和g）相互比较，看看它们之间是如何相互关联的？

对所有组2（'lightblue'b，e和h）和组3（'pink'c，f和i）执行相同的操作。

我尝试使用Tukey的测试（每组的前5个案例）：

TukeyHSD(aov(lm(a[1:5] ~ b[1:5]+c[1:5])),ordered=T,conf.level=0.95)

但是我收到了一个错误：

Error in TukeyHSD.aov(aov(lm(a[1:5] ~ b[1:5] + c[1:5])), ordered = T,  : 
  no factors in the fitted model
In addition: Warning message:
In replications(paste("~", xx), data = mf) : non-factors ignored: b[1:5]

我知道如何使用cut with：

制作因子

cut(a, seq(from=-5, to=4, by=0.2))

我的问题是我不想为我的数据提供范围;相反，我想将第一个箱图中第1组的整个数据与其他两个箱图进行比较。

提前致谢，