nlme 中的多种随机效果

Question

我需要在我的

lme

模型中再包含一个随机效应。然而，这两个随机效应不相关，因此不能嵌套。

我不知道如何为此编写代码。我尝试过以下方法：

modelA <- lme(FT~ Init.Age + Sex + Status, random= (1|Patient) + (1|Site), data = data, na.action = na.omit)

modelB <- lme(FT~ Init.Age + Sex + Status, random= list (~1|Patient, ~1|Site), data = data, na.action = na.omit)

仅型号 B 有效。但它给了我与

random= ~1|Patient

完全相同的结果，我认为这不可能是正确的。结果如下：

感谢您的帮助！ 莉尔

Answer 1

这是

nlme

上一个流行的学习机会，它本质上假设嵌套随机效应。要通过两个因子指定交叉随机截距，我们需要

lme(FT ~ Init.Age + Sex + Status, random = list(one = pdBlocked(
  pdIdent(~ 0 + Patient), pdIdent(~ 0 + Site))), 
  data = data |> transform(one = factor(1)))

请注意，我们必须将

Patient

和

Site

编码为因子而不是整数。有关更多详细信息，请参阅如何在 nlme 与 lme4 中指定不同的随机效应？ 和 https://stats.stackexchange.com/questions/274777/modelling-random-struct-in-lmer-and-nlmelme/613085 .