我有一个时间序列,我用预测包中的ets()和auto.arima()预测。 现在我想创建预测值的平均值。下面的代码有效,但不是很好。例如,我无法自动绘制类似于ets或auto.arima预测的结果。我需要手动重新创建一个类似的图表。
也许有人有更好的方法?
非常感谢!
## auto.arima and ets forecasts
require(forecast)
myactuals <- c(1:36)
myts <- ts(myactuals, start=2011, frequency=12)
myarima <- auto.arima(myts)
myets <- ets(myts)
myarimaforecast <- forecast(myarima, h=12)
myetsforecast <- forecast(myets, h=12)
## merge the two forecasts
mycomposite <- NULL
mycomposite$x <- myts
mycomposite$mean <- (myarimaforecast$mean + myetsforecast$mean)/2