私は R で lmer を使用して混合モデルに取り組んでおり、コーディングに少し固執しています。私は、ショウジョウバエのオスとメスの適応度を、3 つのブロックにわたって 35 の近交系 (遺伝子型) から測定しました。
私の応答変数は「フィットネス」で、n=10 人/性別/ライン/ブロックでテストされています。
性別は固定、ブロックはランダム、ブロック内にネストされたラインはランダムです。私は主にセックスとラインの相互作用に興味があります。したがって、私のモデルは次のようになります
m1<-lmer(FitnessCured~Sex+(1|Block/Line)+(1|Block)+(1|Sex:Line),noNAdata)
相互作用の重要性をテストしたい場合、Sex:Line私の計画は、上記のモデルを相互作用のないモデルanovaと比較し、2 つのモデルを比較するために使用することです。
e.g. m2<-lmer(FitnessCured~Sex+(1|Block/Line)+(1|Block),noNAdata)
anova(m1,m2)
しかし、私が疑問に思っているのは、 Sex:Line 相互作用 (ランダム効果として含まれる) の重要性をテストしている場合、R は Line が Block 内にネストされていることを知っているのでしょうか???
Block 内にネストされた Sex by Line 間の相互作用を指定するにはどうすればよいですか??
それは次のようなものでしょうか
m1T<-lmer(FitnessCured~Sex+(1|Block/Line)+(1|Block)+(1|Sex:Block:Line)
任意の考えをいただければ幸いです。以下に私のデータのサンプルを含めました
Block Line Sex FitnessInfected FitnessCured
2 1 2 M 1.4573 0.2215
3 1 2 M 1.1551 1.1379
4 1 2 M 1.4573 1.1379
7 1 2 M 1.4573 0.4108
9 1 2 M -1.5648 1.1379
11 1 2 F -0.2669 -1.2473
12 1 2 F 0.2785 -1.2473
13 1 2 F -0.5396 -1.2473
14 1 2 F -0.5396 0.4602
15 1 2 F 1.8237 -1.2473
16 1 2 F 0.7330 0.4965
17 1 2 F 1.5511 -1.2473
18 1 2 F -0.5396 1.4774
19 1 2 F 1.0966 1.1868
20 1 2 F -0.5396 -1.2473
21 1 3 M 1.2054 0.7162
22 1 3 M 1.2585 0.3146
24 1 3 M -1.5648 0.2672
26 1 3 M -0.8932 -0.8615
27 1 3 M 0.5047 1.1379
28 1 3 M 0.7704 1.1379
29 1 3 M -1.5648 -1.7689
31 1 3 F -0.5396 0.6782
32 1 3 F -0.5396 -1.2473
33 1 3 F -0.5396 1.0778
34 1 3 F -0.5396 -1.2473
35 1 3 F -0.5396 -1.2473
36 1 3 F -0.5396 0.7145
37 1 3 F -0.5396 0.7508