用Tukey检验校正效应

Question

library(lme4)
library(emmeans)
m <- lmer(angle ~ recipe*temp + (1|replicate), data=cake)

我对每种食谱的简单效果（不是成对比较）感兴趣，因此基于此post

out.emtrends <- emtrends(m, pairwise~recipe, var="temp")
emtr <- as.data.frame(out.emtrends$emtrends)
tvalues <- emtr$temp.trend/emtr$SE
dfs <- emtr$df
pvalues <- 2*pt(-abs(tvalues), dfs)

现在，我想使用Tukey检验来校正我的p值；这是可能的还是Tukey的检查仅适合于治疗之间的差异？

Answer 1

使用pairwise时，调整实际上默认为为土耳其。如果只是通过

显示它，您将在对比结果下方的注释中看到它

out.emtrends

但是，有一个例外，那就是被比较的家庭中只有两个估计值。在这种情况下，只有一个比较。在那种情况下，没有多重性，因此不需要对多重性进行任何形式的调整。