The following function produces in- and out-of-sample MAPA forecasts, for multiple steps ahead. This is the recommended function to use in forecasting with MAPA.
mapafor(y, mapafit, fh = -1, ifh = 1,
comb = c("w.mean","w.median","mean","median","wght"),
outplot = c(0,1), hybrid = c(TRUE, FALSE),
conf.lvl = NULL, xreg=NULL)
In-sample forecasts.
Out-of-sample forecasts.
Prediction intervals for given confidence levels.
In-sample MSE error.
In-sample MAE error.
In sample observations of a time series (vector).
Fitted MAPA model (from mapaest).
Forecast horizon. Default = ppy.
In-sample forecast horizon. Default = 0.
Combination operator. This can be: "mean"; "median"; "wght" - where each aggregation level is weighted inversly to aggregation; "w.mean" - level and trend components are averaged, but seasonal and xreg follow the wght combination; "w.median" - as w.mean, but with median. It is suggested that for data with high sampling frequency to use one of the "w.mean" and "w.median".
Provide output plot. 0 = no; 1 = yes. Default is 1.
Provide hybrid forecasts, as in Kourentzes et al. paper. If minimumAL > 1 then the minimumAL ETS forecasts are used. Default is TRUE.
Vector of confidence level for prediction intervals. Values must be (0,1). If conf.lvl == NULL then no intervals are calculated. For example to get the intervals for 80% and 95% use conf.lvl=c(0.8,0.95).
Vector or matrix of exogenous variables to be included in the MAPA. If matrix then rows are observations and columns are variables. Must be at least as long as in-sample plus fh. Additional observations are unused.
Nikolaos Kourentzes, nikolaos@kourentzes.com; Fotios Petropoulos.
Kourentzes N., Petropoulos F., Trapero J.R. (2014) Improving forecasting by estimating time series structural components across multiple frequencies. International Journal of Forecasting, 30(2), 291--302.
Kourentzes N., Petropoulos F. (2015) Forecasting with multivariate temporal aggregation: The case of promotional modelling. International Journal of Production Economics.
You can find more information about MAPA at Nikos' blog.
mapa
, mapaest
, mapacalc
.
mapafit <- mapaest(admissions,outplot=0)
out <- mapafor(admissions,mapafit)
Run the code above in your browser using DataLab