R augPred 增强预测

说明

预测值是在 primary 的指定值处获得的。如果 object 具有分组结构(即 getGroups(object) 不是 NULL )，则为每个组获取预测值。如果level 具有多个元素，则针对max(level) 分组因子的每一级别获得预测。如果在预测模型中使用除 primary 之外的其他协变量，则使用它们的平均值(数值协变量)或最频繁值(分类协变量)来获取预测值。原始观察结果也包含在返回的对象中。

用法

augPred(object, primary, minimum, maximum, length.out, ...)

## S3 method for class 'lme'
augPred(object, primary = NULL,
        minimum = min(primary), maximum = max(primary),
        length.out = 51, level = Q, ...)

参数

object	使用 predict 方法可以从中提取预测的拟合模型对象。
primary	一个可选的单边公式，指定用于生成增强预测的主要协变量。默认情况下，如果可以从用于生成 object 的数据中提取协变量(使用 getCovariate )，则它将用作 primary 。
minimum	主要协变量的可选下限。默认为 min(primary) 。
maximum	主要协变量的可选上限。默认为 max(primary) 。
length.out	一个可选整数，其中包含用于评估预测的主要协变量值的数量。默认为 51。
level	指定所需预测级别的可选整数向量。级别从最外层到最内层分组递增，级别 0 代表总体(固定效应)预测。默认为最内层。
...	一些泛型方法可能需要额外的参数。

值

具有四列的 DataFrame ，分别表示主协变量的值、组(如果 object 没有分组结构，则所有元素将为 1 )、预测值或观察值以及类型第三列中的值：original 表示观察值，predicted(单个或无分组因子)或 predict.groupVar(多级分组)，其中 groupVar 替换为实际分组变量名称 (fixed用于人口预测)。返回的对象继承自类 "augPred" 。

注意

这个函数是通用的；可以编写方法函数来处理特定类的对象。已经具有此函数方法的类包括： gls 、 lme 和 lmList 。

例子

fm1 <- lme(Orthodont, random = ~1)
augPred(fm1, length.out = 2, level = c(0,1))

作者

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

参考

Pinheiro, J. C. and Bates, D. M. (2000), Mixed-Effects Models in S and S-PLUS, Springer, New York.

也可以看看

plot.augPred , getGroups , predict