r - ggplot2 boxplot の個々のファセットに対して個別の TukeyHSD を実行して注釈を付ける方法は?

Question

私はRがあまり得意ではなく、このコードをまとめようとしていますが、私が望んでいる方法がうまくいきません。これについて何か助けていただければ幸いです。

私のggplotボックスプロットの個々のファセットで治療グループ間でTukeyHSDテストを実行したいと思います。ただし、現在、私の図では、図のすべての箱ひげ図に 1 つの TukeyHSD を適用しているため、図でわかるように膨大な数のグループ化が行われています。

私の現在のプロット

前述したように、TukeyHSD を個別の Depth で区切られた "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))

これが可能であれば、もう 1 つの質問: 最初の行変数 "MBC" と同じように、別の y 変数 "CB" を 2 行目に追加するにはどうすればよいですか?

ご提案ありがとうございます。

        Treatment Depth    MBC        CB