Gráficos para ilustrar los resultados del modelo lineal de efectos mixtos

He estado analizando algunos datos usando el modelado de efectos mixtos lineales en R. Estoy planeando hacer un póster con los resultados y me preguntaba si alguien con experiencia en modelos de efectos mixtos podría sugerir qué gráficos usar para ilustrar los resultados de modelo. Estaba pensando en gráficos residuales, gráfico de valores ajustados frente a valores originales, etc.

Sé que esto dependerá en gran medida de mis datos, pero solo estaba tratando de tener una idea de la mejor manera de ilustrar los resultados de los modelos lineales de efectos mixtos. Estoy usando el paquete nlme en R.

Gracias

Depende de su modelo, pero, en mi experiencia, incluso los colegas, que no tienen una buena comprensión de los modelos de efectos mixtos, realmente les gusta si traza las predicciones con diferentes niveles de agrupación:

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)