-1

我对 R 不是很擅长，我正在尝试将这段代码整合到一起，而这段代码并没有完全按照我的意愿进行。我真的很感激这方面的任何帮助！

我想在我的 ggplot 箱线图中的各个方面对治疗组进行 TukeyHSD 测试。不过，目前，我的图在图中的所有箱线图中应用了一个 TukeyHSD，这会产生大量的分组，如图所示：

我目前的情节

正如我所提到的，最好让 TukeyHSD 在单独的深度分隔“0”方面运行，然后是“5”方面，然后是“30”方面。这可以通过修改我一直使用的代码来实现吗？

data1 <- read.delim(file="clipboard")

data1$Treatment <- as.factor(data1$Treatment)
data1$Depth <- as.factor(data1$Depth)

model<- aov(MBC~Treatment*Depth, data=data1)
model
library(emmeans)
library('multcomp')
cld_dat = as.data.frame( cld(emmeans(model,~Depth*Treatment),
                               Letters = letters ) )

ggplot(data1, aes(x=Treatment, y=MBC, fill=Treatment)) +
  geom_boxplot() + 
  ylab("MBC") +
  ggtitle("Melinis") +
  facet_wrap(~Depth,ncol=3) +
  geom_text(data = cld_dat, aes(y = 140, label = .group))

还有一个问题，如果可能的话：我将如何添加另一个 y 变量“CB”作为第二行，与第一行变量“MBC”的方式相同？

感谢您的任何建议！

        Treatment Depth    MBC        CB
4

1 回答 1

1

如果我正确理解哪个因素具有这些方面级别，那么您需要的是

 cld(emmeans(model, ~ Treatment | Depth))
于 2020-09-07T15:21:43.833 回答