r - 比较 nlme 和 lme4 之间随机截距和斜率无关的脚本

Question

因此，对于随机混合效果，我正在制作两个包之间的脚本比较列表。

对于独立的随机截距和斜率，如果我在 lme4 包中使用以下代码，那么 nlme 中对应的脚本是什么？

model1 <- lmer(y~A + (1+site) + (0+A|site), data, REML = FALSE)

另外，对于嵌套混合效果，它以与上述不同的方式计算随机效果，我的脚本是否正确？

model2 <- lme(y~A, random = ~1+site/A, data, method="REML")

和

model3 <- lmer(y~A + (1|site) + (1|site:A), data, method=FALSE)

太感谢了！

score 0 · Accepted Answer

我希望这个答案还为时不晚！

对于您的第一个模型，版本nlme将是：

model1 <-  lme(y ~ A , 
       random = list(A = pdDiag(~time)),
       data=data)

您的第二个和第三个模型是等效的。封装中的模型 3lme4也可以写成：

 model3 <- lmer(y~A + (1|site/A), data, method=FALSE)

我找到了这个链接，它可能会对您进行比较nlme和lme4打包有很大帮助

1 回答 1