R residuals.lme 提取lme残差

说明

通过从响应向量中减去该级别的拟合级别(并除以估计的组内标准误差，如果是 type="pearson" )即可获得级别 i 的残差。 i 级别的拟合值是通过将总体拟合值(仅基于固定效应估计)和随机效应对小于或等于 i 的分组级别的拟合值的估计贡献相加而获得的。

用法

## S3 method for class 'lme'
residuals(object, level = Q,
          type = c("response", "pearson", "normalized"), asList = FALSE, ...)

参数

object	继承自类 "lme" 的对象，表示拟合的线性混合效果模型。
level	一个可选的整数向量，给出用于从 object 中提取残差的分组级别。级别值从最外层到最内层分组增加，级别零对应于总体残差。默认为最高或最内层的分组。
type	一个可选字符串，指定要使用的残差类型。如果 "response" ，默认情况下，使用 “raw” 残差(观察到的拟合)；否则，如果 "pearson" ，则使用标准化残差(原始残差除以相应的标准误差)；否则，如果 "normalized" ，则使用归一化残差(标准化残差预乘以估计误差相关矩阵的逆平方根因子)。使用参数的部分匹配，因此只需提供第一个字符。
asList	可选的逻辑值。如果 TRUE 并在 level 中给出单个值，则返回的对象是一个列表，其中残差按组划分；否则，根据 level 的长度，返回值是向量或数据帧。默认为 FALSE 。
...	该泛型的某些方法需要额外的参数。此方法中没有使用任何内容。

值

如果在 level 中指定了单级分组，则返回值是按组分割残差的列表 ( asList = TRUE ) 或具有残差的向量 ( asList = FALSE )；否则，当 level 中指定多个分组级别时，返回的对象是一个数据帧，其中的列由不同级别的残差和分组因子给出。对于矢量或数据帧结果，应用 naresid 方法。

例子

fm1 <- lme(distance ~ age + Sex, data = Orthodont, random = ~ 1)
head(residuals(fm1, level = 0:1))
summary(residuals(fm1) /
        residuals(fm1, type = "p")) # constant scaling factor 1.432

作者

José Pinheiro and Douglas Bates bates@stat.wisc.edu

也可以看看

lme , fitted.lme