Diagramme zur Veranschaulichung der Ergebnisse des linearen Mischeffektmodells

Ich habe einige Daten mithilfe der linearen Mischeffektmodellierung in R analysiert. Ich habe vor, ein Poster mit den Ergebnissen zu erstellen, und ich habe mich nur gefragt, ob jemand, der Erfahrung mit Mischeffektmodellen hat, vorschlagen könnte, welche Diagramme zur Veranschaulichung der Ergebnisse verwendet werden sollen Modell. Ich dachte über Residuendiagramme, Diagramm angepasster Werte gegen Originalwerte usw. nach.

Ich weiß, dass dies sehr stark von meinen Daten abhängt, aber ich habe nur versucht, ein Gefühl dafür zu bekommen, wie die Ergebnisse linearer Mischeffektmodelle am besten dargestellt werden können. Ich verwende das nlme-Paket in R.

Vielen Dank

Es hängt von Ihrem Modell ab, aber meiner Erfahrung nach mögen es sogar Kollegen, die kein gutes Verständnis für Modelle mit gemischten Effekten haben, wirklich, wenn Sie die Vorhersagen mit verschiedenen Gruppierungsebenen zeichnen:

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)