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软件包