У меня есть следующая модель
sub <- c(rep(c(1,2,3,4,5,6), each=2, times=4),rep(c(7,8,9,10,11,12), each=2, times=4))
y <- c(12,13,25,24,50,49,75,66,75,77,100,97,23,25,54,56,30,32,55,51,75,73,75,77,50,52,55,57,56,53,48,44,49,48,25,22,99,93,94,93,75,66,30,33,30,26,25,12,25,22,12,17,50,55,50,49,75,77,95,88,48,49,35,36,50,55,55,59,60,67,50,59,25,22,21,21,46,48,67,69,50,59,60,59,22,23,19,20,30,29,50,55,75,69,78,88)
sex <- rep(c("F","M"),each=48)
day <- rep(c("d1","d2","d3","d4"), each=12, times=2)
percent <- c(10,10,14,14,28,28,62,62,83,83,97,97,12,12,30,30,49,49,76,76,77,77,98,98,7,7,23,23,32,32,55,55,75,75,100,100,25,25,50,50,61,61,62,62,70,70,95,95,2,2,14,14,25,25,50,50,68,68,99,99,10,10,25,25,60,60,66,66,82,82,98,98,10,10,27,27,47,47,61,61,76,76,89,89,6,6,24,24,33,33,73,73,75,75,93,93)
dat <- data.frame(sub,y,sex,day,percent)
mod <- lmer(y ~ sex*day*percent + (1|sub), data = dat)
Теперь я хотел бы проверить некоторые контрасты, исправляя количество выполненных контрастов, поэтому я сделал
emmeans(mod, pairwise~sex|day, adjust="tukey")
$contrasts
day = d1:
contrast estimate SE df t.ratio p.value
F - M 0.14332647 5.430965 21.70 0.026 0.9792
day = d2:
contrast estimate SE df t.ratio p.value
F - M -1.81136569 5.386365 21.63 -0.336 0.7399
day = d3:
contrast estimate SE df t.ratio p.value
F - M -0.09560327 5.338141 21.52 -0.018 0.9859
day = d4:
contrast estimate SE df t.ratio p.value
F - M 19.38179381 5.521346 22.07 3.510 0.0020
Однако контрасты не корректируются для нескольких тестов, хотя я указал опцию adjust. Почему это так?