R：使用 nlme 检查模型假设

Question

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)

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

Answer 1

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

oo <- ranef(fm1)$Asym

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

qqnorm(oo)
qqline(oo)

R：使用 nlme 检查模型假设

R: checking model assumptions with nlme

r

diagnostics

nlme