我有一组看起来像这样的数据:
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是饲养温度,而存活是存活的比例 .
我想要做的是,使用R中的lme4包,适合三向交互模型(线性混合模型)来运行ANOVA . 我原来的型号:
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之间 Build 一个3向交互模型,其中line是随机的,rep是一个分组变量
1 回答
重复了很多类似的问题,但是this .
您要问的3向交互模型的代码是:
这是否真的是你需要的,我不知道 .