nlme パッケージの gls() を使用して、mod1 をランダム効果なしでビルドできます。次に、AIC を使用して mod1 を、ランダム効果を含む lme() を使用して構築された mod2 と比較できます。
mod1 = gls(response ~ fixed1 + fixed2, method="REML", data)
mod2 = lme(response ~ fixed1 + fixed2, random = ~1 | random1, method="REML",data)
AIC(mod1,mod2)
lme4 パッケージの gls() に似たものはありますか?これにより、ランダム効果なしで mod3 をビルドし、ランダム効果を含む lmer() を使用してビルドされた mod4 と比較できますか?
mod3 = ???(response ~ fixed1 + fixed2, REML=T, data)
mod4 = lmer(response ~ fixed1 + fixed2 + (1|random1), REML=T, data)
AIC(mod3,mod4)