我在为我的结果生成紧凑的字母显示时遇到问题。 我已经运行了一个方差分析,然后运行了 Tukey 的 HSD 来为每对生成 p 值,但我不知道如何(或者如果可能吗?)将字母分配给这些 p 值以显示哪些对彼此重要。
csa.anova<-aov(rate~temp*light,data=csa.per.chl)
summary(csa.anova)
TukeyHSD(csa.anova)
这运行了我需要的测试,但我不知道如何为每个 p 值分配字母以显示哪些对是重要的。
最佳答案
查找更多详情 here .
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 package 创建于 2021-06-03 (v2.0.0)
关于r - 如何为成对的 TukeyHSD 生成紧凑的字母显示,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/57737598/