我注意到REML方差分量估計的BLUP會系統地導致條件殘差與條件預測圖的上升趨勢。 EMS隨機效應預測始終只是每個隨機效應級別的平均值,因此繪圖總是按預期顯示。 我有點擔心為什麼REML會這樣做。
這是一個示例,我隨機生成100組大小為2的組方差組件。 該模型只是群體隨機效應,沒有固定效應。 如果我保存條件預測和殘差,我會做出如下圖。
用EMS而不是REML做同樣的事情,我得到了這個情節:
這個圖等同於將Group視為固定效果而不是隨機,並且正是我所期望的那樣。
因此,看起來REML BLUPs高估了較低的群體意味著低估群體意味著更高的群體意味著。 對此行為有何解釋?
附帶腳本的數據表。 任何見解將不勝感激!