1

結果のコンパクトな文字表示を生成するのに問題があります。ANOVA に続いて Tukey の HSD を実行して、各ペアの p 値を生成しましたが、これらの p 値に文字を割り当てて、どのペアが互いに有意であるかを示す方法 (または可能かどうか) がわかりません。

csa.anova<-aov(rate~temp*light,data=csa.per.chl) summary(csa.anova) TukeyHSD(csa.anova)

これで必要なテストが実行されますが、各 p 値に文字を割り当ててどのペアが有意であるかを示す方法がわかりません。

4

3 に答える 3

2

詳細はこちらをご覧ください。

mod <- lm(Sepal.Width ~ Species, data = iris)

mod_means_contr <- emmeans::emmeans(object = mod,
                                    pairwise ~ "Species",
                                    adjust = "tukey")

mod_means <- multcomp::cld(object = mod_means_contr$emmeans,
                           Letters = letters)

library(ggplot2)

ggplot(data = mod_means,
       aes(x = Species, y = emmean)) +
  geom_errorbar(aes(ymin = lower.CL, 
                    ymax = upper.CL), 
                width = 0.2) +
  geom_point() +
  geom_text(aes(label = gsub(" ", "", .group)),
            position = position_nudge(x = 0.2)) +
  labs(caption = "Means followed by a common letter are\nnot significantly different according to the Tukey-test")

reprex パッケージ(v2.0.0)により 2021-06-03 に作成

于 2021-06-03T07:33:08.597 に答える