我可以使用 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)