Bayesian Analysis of Vector ARMA Models using Gibbs Sampling

We present a methodology for estimation, prediction, and model assessment of vector autoregressive moving-average (VARMA) models in the Bayesian framework using Markov chain Monte Carlo algorithms. The sampling-based Bayesian framework for inference allows for the incorporation of parameter restrictions, such as stationarity restrictions or zero constraints, through appropriate prior specifications. It also facilitates extensive posterior and predictive analyses through the use of numerical summary statistics and graphical displays, such as box plots and density plots for estimated parameters. We present a method for computationally feasible evaluation of the joint posterior density of the model parameters using the exact likelihood function, and discuss the use of backcasting to approximate the exact likelihood function in certain cases. We also show how to incorporate indicator variables as additional parameters for use in coefficient selection. The sampling is facilitated through a Metropolis–Hastings algorithm. Graphical techniques based on predictive distributions are used for informal model assessment. The methods are illustrated using two data sets from business and economics. The first example consists of quarterly fixed investment, disposable income, and consumption rates for West Germany, which are known to have correlation and feedback relationships between series. The second example consists of monthly revenue data from seven different geographic areas of IBM. The revenue data exhibit seasonality, strong inter-regional dependence, and feedback relationships between certain regions.© 1997 John Wiley & Sons, Ltd.     

Robust seasonal adjustment by Bayesian modelling

Akaike's BAYSEA approach to seasonal decomposition is designed to capture the respective merits of several pre-existing adjustment techniques. BAYSEA is computationally efficient, requires only weak assumptions about the data-generating process, and is based on solid inferential (namely, Bayesian) foundations. We present a model similar to that used in BAYSEA, but based on a double exponential rather than a Gaussian error model. The resulting procedure has the advantages of Akaike's method, but in addition is resistant to outliers. The optimal decomposition is obtained rapidly using a sparse linear programming code. Confidence bands and predictive intervals can be obtained using Gibbs sampling.     

Variance estimation for multivariate dynamic linear models

The problem of estimating unknown observational variances in multivariate dynamic linear models is considered. Conjugate procedures are possible for univariate models and also for special very restrictive common components models but they are not generally applicable. However, for clarity of operation and in order to avoid numerical integration, it is desirable to have conjugacy or approximate conjugacy. Such an approximate procedure is proposed based upon a simple analytic approximation. It is exact for the sub-class of conjugate models and improves on a previous procedure based upon the Robust filter.     

ARIMA model building and the time series analysis approach to forecasting

This paper reviews the approach to forecasting based on the construction of ARIMA time series models. Recent developments in this area are surveyed, and the approach is related to other forecasting methodologies.     

Short‐term forecasting of industrial electricity consumption in Brazil

This paper presents short‐term forecasting methods applied to electricity consumption in Brazil. The focus is on comparing the results obtained after using two distinct approaches: dynamic non‐linear models and econometric models. The first method, that we propose, is based on structural statistical models for multiple time series analysis and forecasting. It involves non‐observable components of locally linear trends for each individual series and a shared multiplicative seasonal component described by dynamic harmonics. The second method, adopted by the electricity power utilities in Brazil, consists of extrapolation of the past data and is based on statistical relations of simple or multiple regression type. To illustrate the proposed methodology, a numerical application is considered with real data. The data represents the monthly industrial electricity consumption in Brazil from the three main power utilities: Eletropaulo, Cemig and Light, situated at the major energy‐consuming states, Sao Paulo, Rio de Janeiro and Minas Gerais, respectively, in the Brazilian Southeast region. The chosen time period, January 1990 to September 1994, corresponds to an economically unstable period just before the beginning of the Brazilian Privatization Program. Implementation of the algorithms considered in this work was made via the statistical software S‐PLUS. Copyright © 1999 John Wiley & Sons, Ltd.     

Time-series analysis supported by power transformations

This paper presents some procedures aimed at helping an applied time-series analyst in the use of power transformations. Two methods are proposed for selecting a variance-stabilizing transformation and another for bias-reduction of the forecast in the original scale. Since these methods are essentially model-independent, they can be employed with practically any type of time-series model. Some comparisons are made with other methods currently available and it is shown that those proposed here are either easier to apply or are more general, with a performance similar to or better than other competing procedures.     

Highest-density forecast regions for nonlinear and non-normal time series models

Forecast regions are a common way to summarize forecast accuracy. They usually consist of an interval symmetric about the forecast mean. However, symmetric intervals may not be appropriate forecast regions when the forecast density is not symmetric and unimodal. With many modern time series models, such as those which are non-linear or have non-normal errors, the forecast densities are often asymmetric or multimodal. The problem of obtaining forecast regions in such cases is considered and it is proposed that highest-density forecast regions be used. A graphical method for presenting the results is discussed.     

Temporal Aggregation in Dynamic Linear Models

One important aspect concerning the analysis and forecasting of time series that is sometimes neglected is the relationship between a model and the sampling interval, in particular, when the observation is cumulative over the sampling period. This paper intends to study the temporal aggregation in Bayesian dynamic linear models (DLM). Suppose that a time series Y_t is observed at time units t and the observations of the process are aggregated over r units of time, defining a new time series Z_k=Σ^r_i=1Y_rk+i. The relevant factors explaining the variation of Z_k can, and in general will, be different, depending on how the sampling interval r is chosen. It is shown that if Y_t follows certain dynamic linear models, then the aggregated series can also be described by possibly different DLM. In the examples, the industrial production of Brazil is analysed under various aggregation periods and the results are compared. © 1997 John Wiley & Sons, Ltd.     

Assessing the historical accuracy of regional economic forecasts

This study addresses for the first time systematic evaluation of a widely used class of forecasts, regional economic forecasts. Ex ante regional structural equation model forecasts are analysed for 19 metropolitan areas. One- to ten-quarter-ahead forecasts are considered and the seven-year sample spans a complete business cycle. Counter to previous speculation in the literature, (1) dependency on macroeconomic forecasting model inputs does not substantially erode accuracy relative to univariate extrapolative methodologies and (2) stochastic time series models do not on average, yield more accurate regional economic predictions than structural models. Similar to findings in other studies, clear preferences among extrapolative methodologies do not emerge. Most general conclusions, however, are subject to caveats based on step-length effects and region-specific effects.     

Multiple time-series modelling: Another look at the canadian money and income data

The practice of modelling the components of a vector time series to arrive at a joint model for the vector is considered. It is shown that in some cases this is not unreasonable. A vector ARMA model is used to model the Canadian money and income data. We also use these data to discuss the issue of differencing a multiple time series. Finally, models based on first and second differences are compared using forecasts.