我一直在使用 R 中的线性混合效应建模分析一些数据。我打算用结果制作一张海报，我只是想知道是否有经验过混合效应模型的人可以建议使用哪些图来说明结果模型。我在考虑残差图、拟合值与原始值的图等。
我知道这在很大程度上取决于我的数据，但我只是想了解说明线性混合效应模型结果的最佳方法。我在 R 中使用 nlme 包。
谢谢
图解说明线性混合效应模型的结果
机器算法验证
r
数据可视化
混合模式
2022-03-09 17:24:57
1个回答
这取决于您的模型，但是根据我的经验，即使是对混合效应模型没有很好理解的同事，也非常喜欢您绘制具有不同分组级别的预测：
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)
