Question

nlme

我正在使用R中的qqnorm(fm1$residuals)包来拟合非线性混合效应模型。我想执行模型诊断并检查1）错误是正态分布的假设以及2）随机效应是否正常分散式。

对于1），我可以只做残差的简单散点图

scinot = re.compile('[-+]?[\d]+\.?[\d]*[Ee](?:[-+]?[\d]+)?')

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

Answer 1

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

oo <- ranef(fm1)$Asym

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

qqnorm(oo)
qqline(oo)