我想使用两个独立的嵌套随机效果来拟合随机效果模型。我可以使用R中的lmer包轻松完成此操作。以下是：
model<-lmer(y ~ 1 + x + (1 | oid/gid) + (1 | did/gid), data=data)
在这里，我拟合oid嵌套在gid和did嵌套在gid内的随机拦截。这很好用。但是，我希望拟合一个模型，其中截距的方差随着gid的变化而变化。 nlme包能够做到这一点。但是，它不清楚如何。我能做的最好就是这样：
model <- lme(y ~ 1 + x, random=list(gid=~1, oid=~1, did=~1), weights=varIdent(form=~1|gid), data = data)
但是这会嵌套在did和oid嵌套在一起的gid。我尝试使用来自类似question的想法，这似乎是一个非常接近的问题，但答案在这个问题中没有得到很好的解释。我希望有人可以解决这个问题。