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)，我只需对残差做一个简单的散点图

qqnorm(fm1$residuals)

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

Answer 1

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

oo <- ranef(fm1)$Asym

现在，您可以使用QQ阴谋来检查正常情况：

qqnorm(oo)
qqline(oo)

​

​