次のような一連のデータがあります。
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)
しかし、その後、私の三者間相互作用はなくなりました。
基本的に、ライン、ステージ、および一時の間の 3 方向の相互作用モデルを作成するのに助けを求めています。ラインはランダムで、担当者はグループ化変数です。