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 是重复的，stage 是我正在研究的果蝇的生命阶段（L = 幼虫，A = 成虫），line 是遗传系的编号，temp 是饲养温度，surv 是存活的比例。

我想要做的是，使用 R 中的 lme4 包，拟合 3 路交互模型（线性混合模型）来运行 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 之间建立一个三向交互模型，其中 line 是随机的，rep 是分组变量

4

1 回答 1

2

许多类似问题的重复，但例如这个。

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

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

这是否真的是你需要的，我不知道。

于 2017-04-20T20:29:33.203 回答