使用预测程序包从预测中获取预测点估计值和时间间隔

Tomas Alonso Rehor 发表于 Dev

托马斯·阿隆索·雷霍（Tomas Alonso Rehor）

我需要一种从预测中打印值的方法。

我需要打印深蓝色的线值，如果可能，还需要打印下图中灰色区域的值。

打印该值或打印2019年预测值的代码是什么？

library(forecast)

timese <- ts(WWWusage, start = c(2008, 1), end = c(2016, 1), frequency = 12)

### Structural Time Series Model 
# Trend likelihood    
fit <- StructTS(timese, "trend")

### Make the plot
plot(forecast(fit, level = c(70, 90)), 
     sub = "Confidence Interval 70% ~ 90% or Determined by user", 
     ylab = "Y Axis Variable",
     main = "Forecast Linear Structural Model @ Trend-Wise",
     ylim = c(0, 400))

杰森·艾兹卡恩斯（JasonAizkalns）

只需存储forecast对象并打印即可：

fc <- forecast(fit, level = c(70, 90))
fc
#          Point Forecast    Lo 70    Hi 70    Lo 90    Hi 90
# Feb 2016            234 230.3083 237.6917 228.1411 239.8589
# Mar 2016            240 231.7450 248.2550 226.8991 253.1009
# Apr 2016            246 232.1868 259.8132 224.0780 267.9220
# May 2016            252 231.7796 272.2204 219.9095 284.0905
# Jun 2016            258 230.6214 285.3786 214.5493 301.4507
# Jul 2016            264 228.7832 299.2168 208.1097 319.8903
# Aug 2016            270 226.3189 313.6811 200.6767 339.3233
# Sep 2016            276 223.2716 328.7284 192.3183 359.6817
# Oct 2016            282 219.6765 344.3235 183.0905 380.9095
# Nov 2016            288 215.5631 360.4369 173.0402 402.9598

要提取单个行，将其转换为可能会更容易data.frame：

df_fc <- as.data.frame(fc)
df_fc["Jul 2016", ]
#          Point Forecast    Lo 70    Hi 70    Lo 90    Hi 90
# Jul 2016            264 228.7832 299.2168 208.1097 319.8903

本文收集自互联网，转载请注明来源。

如有侵权，请联系[email protected] 删除。