R comparePred 比较预测

说明

预测值是在每个对象的 primary 的指定值处获得的。如果 object1 或 object2 具有分组结构(即 getGroups(object) 不是 NULL )，则为每个组获得预测值。当两个对象确定组时，组级别必须相同。如果在预测模型中使用除 primary 之外的其他协变量，则使用它们的分组平均值(数值协变量)或最常见值(分类协变量)来获取预测值。原始观察结果也包含在返回的对象中。

用法

comparePred(object1, object2, primary, minimum, maximum,
    length.out, level, ...)

参数

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

值

具有四列的 DataFrame ，分别表示主协变量的值、组(如果 object 没有分组结构，则所有元素将为 1 )、预测值或观察值以及类型第三列中的值：对象名称用于对预测值进行分类，original 用于对观测值进行分类。返回的对象继承自类 comparePred 和 augPred 。

注意

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

例子

fm1 <- lme(distance ~ age * Sex, data = Orthodont, random = ~ age)
fm2 <- update(fm1, distance ~ age)
comparePred(fm1, fm2, length.out = 2)

作者

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

也可以看看

augPred , getGroups