r - 在 R 中，如何测试仅包含某些级别的因子的混合模型中的交互作用？

Question

我的 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交互的重要性？day1day2

标签： rlme4mixed-models