R 中的 Cochrane-Orcutt 预测答案

【问题标题】：Cochrane-Orcutt Prediction in RR 中的 Cochrane-Orcutt 预测
【发布时间】：2017-11-22 01:19:27
【问题描述】：

我试图在使用 cochrane-orcutt 调整后预测特定值。我的问题是，使用 predict() 或 predict.orcutt()，只会返回使用用于构建模型的原始数据的预测值。下面的例子不是我的数据，但它可以用来复制问题。

# Values used to build model
build_model = mtcars[1:22,]

# Input values for future predictions
predict_values = mtcars[23:32,]

# Linear Model
fit = lm(mpg~cyl, build_model)

# Cochrane.Orcutt correction
coch_class = cochrane.orcutt(fit)

# Use Predict Function 
predict.orcutt(coch_class, newdata = predict_values$cyl)

[,1]
马自达 RX4 20.46226
马自达 RX4 Wag 20.46226
达特桑 710 26.01009
大黄蜂 4 驱动器 20.46226
大黄蜂运动关于 14.91443
勇敢的 20.46226
除尘器 360 14.91443
Merc 240D 26.01009
Merc 230 26.01009
Merc 280 20.46226
Merc 280C 20.46226
Merc 450SE 14.91443
Merc 450SL 14.91443
Merc 450SLC 14.91443
凯迪拉克弗利特伍德 14.91443
林肯大陆 14.91443
克莱斯勒帝国 14.91443
菲亚特 128 26.01009
本田思域 26.01009
丰田卡罗拉 26.01009
丰田 Corona 26.01009
道奇挑战者 14.91443

【问题讨论】：

它也不能作为 data.frame 正常工作（对于这个例子或我的原始数据）。

标签： r statistics

【解决方案1】：

orcutt 库似乎不支持使用新数据进行预测。 predict.orcutt 函数的源仅返回拟合值。它忽略任何参数。

orcutt:::predict.orcutt
# function (object, ...) 
# {
#     return(object$fitted.values)
# }

也许可以联系包的作者，看看他们是否推导出了调整后进行预测的方法。

【讨论】：