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