Ich habe festgestellt, dass BLUPs aus REML-Varianzkomponentenschätzungen systematisch zu einem Aufwärtstrend in einer konditionalen Residual- gegenüber einer prognostizierten Kondition führen. EMS-Zufallseffektvorhersagen sind immer nur der Mittelwert jedes Zufallseffektniveaus, sodass die Diagramme immer wie erwartet erscheinen. Ich kratzte ein wenig darüber, warum REML das tut.
Hier ein Beispiel, bei dem ich zufällig 100 Gruppen der Größe 2 mit einer Gruppenvarianzkomponente generiere. Das Modell ist nur der zufällige Gruppeneffekt, keine festen Effekte. Wenn ich die bedingten Vorhersagen und Residuen speichere, mache ich die folgende Darstellung.
Dasselbe mache ich mit EMS anstelle von REML. Ich bekomme die folgende Darstellung:
Diese Darstellung ist gleichbedeutend mit der Behandlung von Group als festem Effekt statt zufällig, und genau so würde ich es erwarten.
Es sieht also so aus, als würden REML-BLUPs die unteren Gruppenmittelwerte überschätzen und die Gruppenmittelwerte, die höher sind, unterschätzen. Eine Erklärung für dieses Verhalten?
Datentabelle mit Skripten ist beigefügt. Jeder Einblick wäre sehr dankbar!
