Multivariate exponential smoothing: A Bayesian forecast approach based on simulation

This paper deals with the prediction of time series with correlated errors at each time point using a Bayesian forecast approach based on the multivariate Holt–Winters model. Assuming that each of the univariate time series comes from the univariate Holt–Winters model, all of them sharing a common structure, the multivariate Holt–Winters model can be formulated as a traditional multivariate regression model. This formulation facilitates obtaining the posterior distribution of the model parameters, which is not analytically tractable: simulation is needed. An acceptance sampling procedure is used in order to obtain a sample from this posterior distribution. Using Monte Carlo integration the predictive distribution is then approached. The forecasting performance of this procedure is illustrated using the hotel occupancy time series data from three provinces in Spain.