Мой фрейм данных R выглядит примерно так:
subject <- factor(rep(c(1,2,3,4,5,6),each=12))
dep <- runif(12*6)
f1 <- factor(rep(c("Female","Male","Female","Male","Male","Female"), each=12))
f2 <- runif(12*6)
f3 <- factor(rep(c("day1","day2","day3","day4"),each=3, times=6))
data <- data.frame(sub=subject, dep=dep, f1=f1, f2=f2, f3=f3)
m <- lmer(dep ~ f1*f2*f3 + (1|sub), data=data)
Это случайные данные, поэтому нет ничего существенного, но, скажем, существует значительное f1*f2*f3 взаимодействие.Теперь, чтобы понять, откуда происходит это взаимодействие, я хотел бы выполнить более простые взаимодействия, включающие только некоторые уровни f3.Например, как я могу проверить значимость взаимодействия f1*f2*f3, включающего только day1 и day2?