首页 文章

随机和固定效果之间的3种相互作用

提问于
浏览
0

我有一组看起来像这样的数据:

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 回答

  • 0

    重复了很多类似的问题,但是this .

    您要问的3向交互模型的代码是:

    surv_3w.aov<-lmer(surv~stage*line*temp + (1 + line |rep), data=dat_3w)
    

    这是否真的是你需要的,我不知道 .

相关问题