使用 emmeans 的 2 路和 3 路交互的置信区间

Question

这只是一个关于在 emmeans 中获得交互置信区间的一般问题，我已经阅读了所有常见的教程，但我无法理解如何为 2 路和 3 路交互做到这一点。 这是 3 向交互的示例。

X= continous variable   
A= 3 levels, nested in L  
L= 2 levels  
G= 2 levels, crossed with L (and thus also A)  
ID= participants nested in G

然后我用随机截距拟合了一个线性混合模型：

fit1 <- lmer(X~G*A*L+(1|ID),data=df)
emms=emmeans(fit1,specs=pairwise~A*L*G)

我想要的比较使用以下方法完美显示：

contrast(emms,interaction="pairwise")$emmeans

但是，它们缺少置信区间，请注意运行：

summary(contrast(emms,interaction="pairwise")$emmeans,infere=T) or
confint(contrast(emms,interaction="pairwise")$emmeans) 

在emmeans数据上不起作用，它只是给出了不同级别的emmeans和置信区间，而不是对比。 这些函数对对比数据起作用，但不显示三向交互。

更新：答案我终于想通了：

confint(contrast(emmeans(fit1,~A*G*L),interaction=c("pairwise")))

Answer 1

confint(contrast(emmeans(fit1,~A*G*L),interaction=c("pairwise")))

Answer 2

另一个答案产生的交互对比可能是您想要的，但这些对比包括三向模型中的对比对比。 当存在交互时，人们通常想要的是简单的比较——比较一个因素而其他因素保持不变。 这些可以通过

pairs(emms[[1]], simple = “each”)