1.6 将输出用为输入——结果的重用

1.6 将输出用为输入——结果的重用
R的一个非常实用的特点是,分析的输出结果可轻松保存,并作为进一步分析的输入使用。
让我们通过一个R中已经预先安装好的数据集作为示例阐明这一点。如果你无法理解这里涉及的
统计知识,也别担心,我们在这里关注的只是一般原理。
首先,利用汽车数据mtcars执行一次简单线性回归,通过车身重量(wt)预测每加仑行驶
的英里数(mpg)。可以通过以下语句实现:

lm(mpg~wt, data=mtcars)
结果将显示在屏幕上,不会保存任何信息。
下一步,执行回归,区别是在一个对象中保存结果:

lmfit <- lm(mpg~wt, data=mtcars)
以上赋值语句创建了一个名为lmfit的列表对象,其中包含了分析的大量信息(包括预测值、残
差、回归系数等)。虽然屏幕上没有显示任何输出,但分析结果可在稍后被显示和继续使用。
键入summary(lmfit)将显示分析结果的统计概要,plot(lmfit)将生成回归诊断图形,
而语句cook<-cooks.distance(lmfit)将计算影响度量统计量②,plot(cook)对其绘图。要
在新的车身重量数据上对每加仑行驶的英里数进行预测,不妨使用predict(lmfit, mynewdata)。
要了解某个函数的返回值,查阅这个函数在线帮助文档中的“Value”部分即可。本例中应
当查阅help(lm)或?lm中的对应部分。这样就可以知道将某个函数的结果赋值到一个对象时,
保存下来的结果具体是什么。

posted @ 2016-08-28 16:14  捉老鹰的小鸡鸡  阅读(261)  评论(0编辑  收藏  举报