将季度预测r转换为年度预测

时间:2017-12-10 20:33:15

标签: r forecasting forecastr

我正在使用library(forecast)根据季度数据创建季度预测。我的预测保存在forecast对象中。我试图找到一种优雅而简单的方法,将它们转换为相应的年度值,将所有4个季度的平均值罚款。

forecast包中是否有允许执行此操作的选项?或者,如果有办法改变季度值?

我尝试将forecast对象转换为timeseries对象,并将其与我用于创建预测的原始时间序列合并但不起作用。

以下是一个例子:

library('quantmod')
library('xts')
library('forecast')
library(zoo)
library(tis)
library(astsa)
library(xts)
GDP_SAAR<-getSymbols('A191RL1Q225SBEA',src='FRED', auto.assign=FALSE)
GDP_SAAR   <- ts(GDP_SAAR , start=c(1947, 2), end=c(2017, 3), frequency=4) 
fit <- auto.arima(GDP_SAAR)
x<-forecast(fit,h=7)

1 个答案:

答案 0 :(得分:1)

我不知道您是否仍然感兴趣,但这是一个涉及12个步骤(季度)的解决方案,以包含至少3年的汇总预测:

GDP_SAAR<-getSymbols('A191RL1Q225SBEA',src='FRED', auto.assign=FALSE)
fit <- auto.arima(GDP_SAAR)
x<-forecast(fit,h=12)
DateTime=as.Date(seq(from=tail(index(GDP_SAAR),1),by="quarter",length.out = length(x$mean)+1))
myresult=tibble(Forecasts=x$mean,DateTime=DateTime[-1])%>%mutate(DateTime=year(DateTime))%>%
  group_by(DateTime)%>%summarize(Forecasts=sum(Forecasts,na.rm=T))
> myresult
# A tibble: 3 x 2
  DateTime Forecasts
     <int>     <dbl>
1     2018      10.5
2     2019      10.6
3     2020      10.6

希望这有帮助