我注意到REML方差分量估计的BLUP会系统地导致条件残差与条件预测图的上升趋势。 EMS随机效应预测始终只是每个随机效应级别的平均值,因此绘图总是按预期显示。 我有点担心为什么REML会这样做。
这是一个示例,我随机生成100组大小为2的组方差组件。 该模型只是群体随机效应,没有固定效应。 如果我保存条件预测和残差,我会做出如下图。
用EMS而不是REML做同样的事情,我得到了这个情节:
这个图等同于将Group视为固定效果而不是随机,并且正是我所期望的那样。
因此,看起来REML BLUPs高估了较低的群体意味着低估群体意味着更高的群体意味着。 对此行为有何解释?
附带脚本的数据表。 任何见解将不胜感激!