REMLの分散成分の推定値からのBLUPは、条件付き残差対条件付き予測プロットで系統的に上昇傾向になることに気付きました。 EMSのランダム効果予測は常に各ランダム効果レベルの平均値であるため、プロットは常に予想通りに表示されます。 私はなぜREMLがこれをするのかについて少し頭をかいています。
これは、グループ分散コンポーネントを使用してサイズ2のグループを100個ランダムに生成する例です。 モデルは単なるグループランダム効果であり、固定効果ではありません。 条件付き予測と残差を保存すると、次のようになります。
REMLの代わりにEMSでも同じことをすると、このプロットが得られます。
このプロットは、Groupをランダム効果ではなく固定効果として扱うことと同じです。まさにそのように見えるはずです。
そのため、REML BLUPは低いグループ平均を過大評価し、高いグループ平均を過小評価しているようです。 この動作についての説明はありますか？
スクリプト付きデータテーブルが添付されています。 任意の洞察力は大歓迎です！
