Diagrammes pour illustrer les résultats du modèle linéaire à effets mixtes

J'ai analysé certaines données à l'aide de la modélisation linéaire à effets mixtes dans R. Je prévois de faire une affiche avec les résultats et je me demandais simplement si quelqu'un expérimenté avec les modèles à effets mixtes pourrait suggérer les tracés à utiliser pour illustrer les résultats de la modèle. Je pensais aux graphiques résiduels, au graphique des valeurs ajustées par rapport aux valeurs d'origine, etc.

Je sais que cela dépendra beaucoup de mes données, mais j'essayais juste d'avoir une idée de la meilleure façon d'illustrer les résultats des modèles à effets mixtes linéaires. J'utilise le paquet nlme dans R.

Merci

Cela dépend de votre modèle, mais, selon mon expérience, même mes collègues, qui n'ont pas une bonne compréhension des modèles d'effets mixtes, aiment vraiment si vous tracez les prédictions avec différents niveaux de regroupement:

library(nlme)
fm2 <- lme(distance ~ age + Sex, data = Orthodont, random = ~ 1|Subject)

newdat <- expand.grid(Sex=unique(Orthodont$Sex),
                  age=c(min(Orthodont$age),
                            max(Orthodont$age)))

library(ggplot2)
p <- ggplot(Orthodont, aes(x=age, y=distance, colour=Sex)) +
  geom_point(size=3) +
  geom_line(aes(y=predict(fm2), group=Subject, size="Subjects")) +
  geom_line(data=newdat, aes(y=predict(fm2, level=0, newdata=newdat), size="Population")) +
  scale_size_manual(name="Predictions", values=c("Subjects"=0.5, "Population"=3)) +
  theme_bw(base_size=22) 
print(p)