我有一组看起来像这样的数据：
rep stage line temp surv
1 L 149 18 0.6
2 L 149 18 0.7
3 L 149 18 0.25
1 A 149 18 1
2 A 149 18 1
3 A 149 18 1
1 L 149 25 0
2 L 149 25 0.2
3 L 149 25 0.3
1 A 149 25 1
2 A 149 25 1
3 A 149 25 1
1 L 208 18 0.6
2 L 208 18 0.4
3 L 208 18 0.55
1 A 208 18 1
2 A 208 18 1
3 A 208 18 1
1 L 208 25 0
2 L 208 25 0.05
3 L 208 25 0.05
1 A 208 25 1
2 A 208 25 0.857142857
3 A 208 25 0.7
当rep复制时，阶段是我正在使用的果蝇的生命阶段（L =幼虫，A =成虫），线是遗传系的数字分配，temp是饲养温度，而surv是比例幸存下来。
我想要做的是，使用R中的lme4包，适合运行ANOVA的3向交互模型（线性混合模型）。我原来的模特：
surv_3w.aov<-lmer(surv~stage*line*temp + (1|rep), data=dat_3w)
有效，但我想把线视为随机效果。我认为我正确地将rep视为分组变量(1|rep)，但我不确定。
我试过这个模型：
surv_3w.aov<-lmer(surv~stage*temp*(1|line) + (1|rep), data=dat_3w)
然后我的三方互动消失了。
基本上，我要求帮助在line，stage和temp之间建立一个3向交互模型，其中line是随机的，rep是一个分组变量
答案 0 :(得分：0)
重复了很多类似的问题，但例如this。
您要问的3向交互模型的代码是：
surv_3w.aov<-lmer(surv~stage*line*temp + (1 + line |rep), data=dat_3w)
这是否真的是你需要的，我不知道。