问传递预报()和ar()中的参数时R中的预测错误

EN

这是forecast.ar()中的一个问题(实际上不是一个bug)。所有forecast.xxx()函数都试图存储用于估计时间序列模型的数据，因为它是绘图和精确计算所需的。但是，ar()不返回数据，因此forecast.ar()试图在调用环境或父环境中查找数据。当您直接调用forecast(ar(...))时，函数成功地找到了数据，但是arf()将对ar()的调用放在了一个更深的层次，这使得forecast很难知道正在使用什么数据。

我可以修改函数，使其更难查找数据(也可以查看祖父母环境)，但是构造仍然会失败，因为predict.ar() ( stats包的一部分)会导致类似的错误。如果ar()返回数据会好得多，但是ar()是stats包的一部分，我无法控制它。

有几种可能的解决办法。

1. 您可以将ar替换为Arima： arf <-函数(dat，.)预报(Arima(dat，order=c(1,0,0))，.) 如果数据是平稳的，这应该返回相同的模型(尽管参数估计会略有不同)。对于问题中的例子，它不会返回相同的答案，因为时间序列是非平稳的。
2. 如果您愿意使用比AR(1)更通用的ARIMA模型，您可以使用auto.arima()。 arf <-函数(dat，.)预报(auto.arima(dat，.)
3. (根据@agstudy的建议)。解决方法是确保数据存储在ar对象中： arf <- function( dat，.){ object <- ar(dat，order.max=1，method="ols")对象$x <- dat预测(object，.)}

问题是S3方法对ar类的预测存在缺陷。predict.ar尝试使用ar对象计算新数据参数。显示getS3method('predict','ar')的第一行

function (object, newdata, n.ahead = 1L, se.fit = TRUE, ...) 
{
    if (n.ahead < 1L) 
        stop("'n.ahead' must be at least 1")
    if (missing(newdata)) {
        newdata <- eval.parent(parse(text = object$series))
        if (!is.null(nas <- object$call$na.action)) 
            newdata <- eval.parent(call(nas, newdata))
    }

  .....
}

有关/窃听的线路是：

 newdata <- eval.parent(parse(text = object$series))

但是对象$series没有正确的表达式/字符。因为它是隐藏在新级别的包装函数arf。这里有一个解决办法，就是为这个术语设置正确的表达式：

arf <- function(dat, ...) 
  {
    object <- ar(dat, order.max=1, method="ols")
    object$series <- as.character(as.expression(as.list(match.call())$dat))
    forecast(object,...)
}
arf( ts(1:100, start=c(2000,1), frequency=4)

请注意；此解决方案还适用于：

aa <- ts(1:100, start=c(2000,1), frequency=4)
arf(aa)