Укажите несколько вложенных случайных эффектов в Im e с гетероскедастической дисперсией по группе
Я хочу, чтобы модель случайных эффектов соответствовала двум отдельным вложенным случайным эффектам. Я могу легко сделать это, используя lmer Пакет в R. Вот как:
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 вложенный вместе Я попытался использовать идею из подобного вопроса, который кажется близкой проблемой, но ответ на этот вопрос не был хорошо объяснен. Я надеюсь, что кто-то может понять это.