R：使用 nlme 检查模型假设答案

【问题标题】：R: checking model assumptions with nlmeR：使用 nlme 检查模型假设
【发布时间】：2023-04-07 23:50:01
【问题描述】：

library(nlme)
fm1 <- nlme(height ~ SSasymp(age, Asym, R0, lrc),
            data = Loblolly,
            fixed = Asym + R0 + lrc ~ 1,
            random = Asym ~ 1,
            start = c(Asym = 103, R0 = -8.5, lrc = -3.3))

我正在使用 R 中的 nlme 包拟合非线性混合效应模型。我想执行模型诊断并检查假设 1) 错误是正态分布的，并且 2) 随机效应是正态分布的。

对于 1)，我可以做一个简单的残差散点图

qqnorm(fm1$residuals)

但是如何检查随机效应是否也呈正态分布？

【问题讨论】：

标签： r diagnostics nlme

【解决方案1】：

您可以使用通用函数ranef（或random.effects）提取随机效果。特别是你想要的效果是

oo <- ranef(fm1)$Asym

现在您可以使用 QQ 图来检查正态性：

qqnorm(oo)
qqline(oo)

【讨论】：