Sơ đồ để minh họa kết quả của mô hình hiệu ứng hỗn hợp tuyến tính

Tôi đã phân tích một số dữ liệu bằng cách sử dụng mô hình hiệu ứng hỗn hợp tuyến tính trong R. Tôi dự định tạo một poster với kết quả và tôi chỉ tự hỏi liệu có ai có kinh nghiệm với các mô hình hiệu ứng hỗn hợp có thể đề xuất sử dụng lô nào để minh họa kết quả của mô hình. Tôi đã suy nghĩ về các lô dư, lô của các giá trị được trang bị so với các giá trị ban đầu, v.v.

Tôi biết điều này sẽ phụ thuộc rất nhiều vào dữ liệu của tôi nhưng tôi chỉ cố gắng cảm nhận cách tốt nhất để minh họa kết quả của các mô hình hiệu ứng hỗn hợp tuyến tính. Tôi đang sử dụng gói nlme trong R.

Cảm ơn

Nó phụ thuộc vào mô hình của bạn, nhưng, theo kinh nghiệm của tôi, ngay cả các đồng nghiệp, những người không hiểu rõ về các mô hình hiệu ứng hỗn hợp, thực sự thích nếu bạn vẽ các dự đoán với các mức nhóm khác nhau:

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)