R predict.lme 来自 lme 对象的预测

说明

i 级别的预测是通过将总体预测(仅基于固定效应估计)和随机效应对小于或等于 i 的分组级别的预测的估计贡献相加而获得的。结果值估计 i 级别的最佳线性无偏预测 (BLUP)。如果 newdata 中存在未包含在原始分组因子中的组值，则对于大于或等于未知组出现的级别的级别，相应的预测将设置为 NA。

用法

## S3 method for class 'lme'
predict(object, newdata, level = Q, asList = FALSE,
        na.action = na.fail, ...)

参数

object	继承自类 "lme" 的对象，表示拟合的线性混合效果模型。
newdata	用于获得预测的可选数据帧。固定效应模型和随机效应模型中使用的所有变量以及分组因子都必须出现在 DataFrame 中。如果缺失，则返回拟合值。
level	一个可选的整数向量，给出用于获得预测的分组级别。级别值从最外层到最内层分组递增，级别零对应于总体预测。默认为最高或最内层的分组。
asList	可选的逻辑值。如果 TRUE 并在 level 中给出单个值，则返回的对象是一个列表，其中预测按组划分；否则，根据 level 的长度，返回值是向量或数据帧。
na.action	一个函数，指示当 newdata 包含 NA 时应该发生什么。默认操作 ( na.fail ) 会导致函数打印一条错误消息，并在存在任何不完整的观察结果时终止。
...	该泛型的某些方法需要额外的参数。此方法中没有使用任何内容。

值

如果在 level 中指定了单个分组级别，则返回值是包含按组分割的预测的列表 ( asList = TRUE ) 或包含预测的向量 ( asList = FALSE )；否则，当在 level 中指定多个分组级别时，返回的对象是一个数据帧，其中的列由不同级别的预测和分组因子给出。

例子

fm1 <- lme(distance ~ age, Orthodont, random = ~ age | Subject)
newOrth <- data.frame(Sex = c("Male","Male","Female","Female","Male","Male"),
                      age = c(15, 20, 10, 12, 2, 4),
                      Subject = c("M01","M01","F30","F30","M04","M04"))
## The 'Orthodont' data has *no* 'F30', so predict  NA  at level 1 :
predict(fm1, newOrth, level = 0:1)

作者

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

也可以看看

lme , fitted.lme