我想用两个独立的嵌套随机效应拟合随机效应模型。lmer我可以使用R 中的包轻松完成此操作。方法如下：
model<-lmer(y ~ 1 + x + (1 | oid/gid) + (1 | did/gid), data=data)
在这里，我为oidnested withingid和didnested within拟合随机截距gid。这很好用。但是，我想拟合一个模型，其中截距的方差随gid随机效应的变化而变化。nlmepackage 能够做到这一点。但是，目前尚不清楚如何。我能做的最好的就是这样：
model <- lme(y ~ 1 + x, random=list(gid=~1, oid=~1, did=~1), weights=varIdent(form=~1|gid), data = data)
但这将did内部嵌套在一起oid。gid我试图从一个类似的问题中使用这个想法，这似乎是一个很接近的问题，但在那个问题中答案并没有得到很好的解释。我希望有人能解决这个问题。