Chapter Contents |
Previous |
Next |

Time Series Analysis and Control Examples |

Selected routines from the TIMSAC package from ISM were converted by SAS Institute staff into SAS/IML routines under an agreement between SAS Institute and ISM. Credit for authorship of these TIMSAC SAS/IML routines goes to ISM, which has agreed to make them available to SAS users without charge.

There are four packages of TIMSAC programs. See the section "ISM TIMSAC Packages" for more information on the TIMSAC package produced by ISM. Since these SAS/IML time series analysis subroutines are adapted from the corresponding FORTRAN subroutines in the TIMSAC package produced by ISM, they are collectively referred to as "the TIMSAC subroutines" in this chapter.

The subroutines analyze and forecast univariate and multivariate time series data. The nonstationary time series and seasonal adjustment models can also be analyzed by using the Interactive Matrix Language TIMSAC subroutines. These subroutines contain the Bayesian modeling of seasonal adjustment and changing spectrum estimation.

Discrete time series modeling has been widely used to analyze dynamic systems in economics, engineering, and statistics. The Box-Jenkins and Box-Tiao approaches are classical examples of unified time series analysis through identification, estimation, and forecasting (or control). The ARIMA procedure in the SAS/ETS product uses these approaches. Bayesian methods are being increasingly applied despite the controversial issues involved in choosing a prior distribution.

The fundamental idea of the Bayesian method is that
uncertainties can be explained by probabilities.
If there is a class model consisting of sets of
member models (), you can describe the uncertainty
of using a prior distribution of .
The member model is directly related to model parameters.
Let the prior probability density function be .When you observe the data **y** that is generated from
the model , the data distribution is described
as given the unknown with a
prior probability density , where the
function is the usual likelihood function.
Then the posterior distribution is the updated
prior distribution given the sample information.
The posterior probability density function
is proportional to *observed likelihood
function × prior density function*.

The IML TIMSAC subroutines contain various time series analysis and Bayesian models. Most of the subroutines are based on the minimum Akaike Information Criterion (AIC) or on the minimum ABIC method to determine the best model among alternative models. The TSBAYSEA subroutine is a typical example of Bayesian modeling. The following subroutines are supported:

**TSBAYSEA**- Bayesian seasonal adjustment modeling
**TSDECOMP**- time series decomposition analysis
**TSMLOCAR**- locally stationary univariate AR model fitting
**TSMLOMAR**- locally stationary multivariate AR model fitting
**TSMULMAR**- multivariate AR model fitting
**TSPERARS**- periodic AR model fitting
**TSPRED**- ARMA model forecasting and forecast error variance
**TSROOT**- polynomial roots or ARMA coefficients computation
**TSTVCAR**- time-varying coefficient AR model estimation
**TSUNIMAR**- univariate AR model fitting

The TSUNIMAR and TSMULMAR subroutines estimate the autoregressive models and select the appropriate AR order automatically by using the minimum AIC method. The TSMLOCAR and TSMLOMAR subroutines analyze the nonstationary time series data. The Bayesian time-varying AR coefficient model (TSTVCAR) offers another nonstationary time series analysis method. The state space and Kalman filter method is systematically applied to the smoothness priors models (TSDECOMP and TSTVCAR), which have stochastically perturbed difference equation constraints. The TSBAYSEA subroutine provides a way of handling Bayesian seasonal adjustment, and it can be an alternative to the SAS/ETS X-11 procedure. The TSBAYSEA subroutine employs the smoothness priors idea through constrained least squares estimation, while the TSDECOMP and TSTVCAR subroutines estimate the smoothness tradeoff parameters using the state space model and Kalman filter recursive computation. The TSPRED subroutine computes the one-step or multi-step predicted values of the ARMA time series model. In addition, the TSPRED subroutine computes forecast error variances and impulse response functions. The TSROOT subroutine computes the AR and MA coefficients given the characteristic roots of the polynomial equation and the characteristic roots for the AR or MA model.

Chapter Contents |
Previous |
Next |
Top |

Copyright © 1999 by SAS Institute Inc., Cary, NC, USA. All rights reserved.