J'ai remarqué que les BLUP issus des estimations de la composante de variance REML entraîneraient systématiquement une tendance à la hausse dans un tracé résiduel conditionnel par rapport au prédictif conditionnel. Les prédictions d'effets aléatoires EMS ne représentent toujours que la moyenne de chaque niveau d'effet aléatoire, les graphes apparaissent donc toujours comme prévu. Je me suis un peu demandé pourquoi REML faisait cela.
Voici un exemple où je génère de manière aléatoire 100 groupes de taille 2 avec une composante de variance de groupe. Le modèle est juste l'effet de groupe aléatoire, pas d'effets fixes. Si je sauvegarde les prévisions et les résidus conditionnels, je trace le graphique suivant.
En faisant la même chose avec EMS au lieu de REML, je reçois cette intrigue:
Ce graphique équivaut à traiter Groupe comme un effet fixe plutôt que aléatoire, et correspond exactement à ce à quoi il devrait ressembler.
Donc, il semble que REML BLUP surestime les moyennes du groupe inférieur et sous-estime les moyennes du groupe qui sont plus élevées. Une explication à ce comportement?
La table de données avec les scripts est jointe. Toute idée serait grandement apprécié!
