Доверительные интервалы двухстороннего и трехстороннего взаимодействия с использованием emmeans

Question

Это всего лишь общий вопрос о получении доверительных интервалов для взаимодействий в emmeans. Я прочитал все распространенные учебные пособия, но не могу понять, как это сделать для двусторонних и трехсторонних взаимодействий. Вот пример для трехстороннего взаимодействия.

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

Контрасты взаимодействия, создаваемые другим ответом, могут быть тем, что вы хотите, но они включают контрасты контрастов в трехсторонней модели. Когда присутствуют взаимодействия, часто то, что люди хотят, - это простые сравнения - сравнения одного фактора, в то время как другие остаются фиксированными. Они доступны через

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

Answer 2

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