我正在学习https://www.cmm.bris.ac.uk/lemma/login/index.php的混合模型，这个问题是关于模块5中的练习 .
该数据来自508所学校的33,988次观察 . 学生数据是第一级，学校是第二级 . 我运行null模型
nullmodel <- lmer(score ~ (1|schoolid), data = mydata, REML = FALSE)
（我找不到对数据集的开放访问权限 . 我想因为这个问题更具概念性，所以无关紧要）
我们通过运行nullmodel的摘要得到的截距估计是否应该等于'得分'的平均值？我得到拦截估计值为30.60，平均值[均值（mydata $得分）]为31.09 . 值很接近但不相同 .
另外，有没有直接的方法来获得学校级残差u0？
1 回答
随机截距效应u0通常以均值0和方差σ_u0^ 2分布 .
σ_u0^ 2是校际差异 .
您可以使用直接方式计算σ_u0^ 2 .
计算每所学校的平均分数
计算这些值的方差
方差接近σ_u0^ 2 .