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)

非常感谢！

Answer 1

希望这个答案不会太晚！

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

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

您的第二个和第三个模型是等效的。 lme4软件包中的Model 3也可以写成：

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

我发现此链接可能对您比较nlme和lme4软件包有很大帮助