我试图估计 R 中的纵向模型，其中有几个相互关联的随机截距，并且数据是非嵌套的。例如，考虑一个简单的纵向模型，其中学生的考试成绩在三个随机截距、一个随机学校效应、一个随机教师效应和一个随机教师-学校匹配效应上回归。数据是非嵌套的，因为教师在纵向数据中经常更换学校。我主要对随机效应参数本身感兴趣，因此我想指定所有三个随机效应并允许它们相互关联。没有随机斜率，只有随机截距。
我试过了：
model1 <- lmer(test_score ~ FEs + (1|schoolid) + (1|teacherid) + (1|matchid), data)
但这似乎迫使随机效应彼此不相关。有没有办法，使用lmer或任何其他 R 包来估计一个类似的模型，其中允许随机截距相互关联（使用非嵌套数据）？
更具体地说，上述模型强制协方差矩阵是对角的。我想估计的是一个具有无限制协方差矩阵的随机效应模型，例如这个链接