Rでのlmer / lme混合モデルの仮定の確認
3つの異なるタスクで30人の男性と30人の女性をテストする繰り返しデザインを実行しました。男性と女性の行動がどのように異なり、それがタスクにどのように依存するかを理解したいと思います。これを調査するためにlmerとlme4の両方のパッケージを使用しましたが、いずれかの方法の仮定を確認しようとしています。私が実行するコードは lm.full <- lmer(behaviour ~ task*sex + (1|ID/task), REML=FALSE, data=dat) lm.full2 <-lme(behaviour ~ task*sex, random = ~ 1|ID/task, method="ML", data=dat) 相互作用のない単純なモデルと比較し、anovaを実行して、相互作用が最良のモデルであるかどうかを確認しました。 lm.base1 <- lmer(behaviour ~ task+sex+(1|ID/task), REML=FALSE, data=dat) lm.base2 <- lme(behaviour ~ task+sex, random= ~1|ID/task), method="ML", data=dat) anova(lm.base1, lm.full) anova(lm.base2, lm.full2) Q1:これらのカテゴリカル予測子を線形混合モデルで使用しても大丈夫ですか? Q2:結果変数(「振る舞い」)がそれ自体(性別/タスク間)自体に正規分布する必要はないことを正しく理解していますか? Q3:分散の均一性を確認するにはどうすればよいですか?単純な線形モデルでは、を使用しますplot(LM$fitted.values,rstandard(LM))。plot(reside(lm.base1))十分に使用していますか? Q4:正常性を確認するには、次のコードを使用しますか? hist((resid(lm.base1) - mean(resid(lm.base1))) / sd(resid(lm.base1)), …