使用Tukey's测试校正效果

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

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