r - 使用 Tukey 检验校正效果

Question

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

我对每个食谱的简单效果（不是成对比较）感兴趣，所以基于这篇文章

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

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