Графики для иллюстрации результатов модели линейного смешанного эффекта

Я анализировал некоторые данные с использованием линейного моделирования смешанных эффектов в R. Я планирую сделать плакат с результатами, и мне было просто интересно, может ли кто-нибудь, имеющий опыт работы с моделями смешанных эффектов, предложить, какие графики использовать для иллюстрации результатов модель. Я думал об остаточных графиках, графике подгоночных значений против исходных значений и т. Д.

Я знаю, что это будет очень сильно зависеть от моих данных, но я просто пытался понять, как лучше всего проиллюстрировать результаты линейных моделей со смешанным эффектом. Я использую пакет nlme в R.

благодаря

Это зависит от вашей модели, но, по моему опыту, даже коллегам, которые плохо разбираются в моделях смешанных эффектов, очень нравится, когда вы строите прогнозы с разными уровнями группировки:

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)