The ARIMA Procedure

General Notation for ARIMA Models

The order of an ARIMA (autoregressive integrated moving-average) model is usually denoted by the notation ARIMA(p,d,q ), where

p

is the order of the autoregressive part

d

is the order of the differencing

q

is the order of the moving-average process

If no differencing is done (d = 0), the models are usually referred to as ARMA(p, q) models. The final model in the preceding example is an ARIMA(1,1,1) model since the IDENTIFY statement specified d = 1, and the final ESTIMATE statement specified p = 1 and q = 1.

Notation for Pure ARIMA Models

Mathematically the pure ARIMA model is written as

\[  W_{t}={\mu }+\frac{{\theta }({B})}{{\phi }({B})}a_{t}  \]

where

t

indexes time

${W_{t}}$

is the response series ${Y_{t}}$ or a difference of the response series

${\mu }$

is the mean term

$B$

is the backshift operator; that is, ${{B}X_{t}=X_{t-1}}$

${{\phi }({B})}$

is the autoregressive operator, represented as a polynomial in the backshift operator: ${{\phi }({B})=1-{\phi }_{1}{B}-{\ldots }-{\phi }_{p}{B}^{p}}$

${{\theta }({B})}$

is the moving-average operator, represented as a polynomial in the backshift operator: ${{\theta }({B})=1-{\theta }_{1}{B}-{\ldots }-{\theta }_{q}{B}^{q}}$

${a_{t}}$

is the independent disturbance, also called the random error

The series ${W_{t}}$ is computed by the IDENTIFY statement and is the series processed by the ESTIMATE statement. Thus, ${W_{t}}$ is either the response series Y$_{t}$ or a difference of ${Y_{t}}$ specified by the differencing operators in the IDENTIFY statement.

For simple (nonseasonal) differencing, ${W_{t}=(1-{B})^{d}Y_{t}}$. For seasonal differencing ${W_{t}=(1-{B})^{d}(1-{B}^{s})^{D}Y_{t}}$, where d is the degree of nonseasonal differencing, D is the degree of seasonal differencing, and s is the length of the seasonal cycle.

For example, the mathematical form of the ARIMA(1,1,1) model estimated in the preceding example is

\[  (1-{B})Y_{t}={\mu }+\frac{(1-{\theta }_{1}{B})}{(1-{\phi }_{1}{B})}a_{t}  \]

Model Constant Term

The ARIMA model can also be written as

\[  {\phi }({B})(W_{t}-{\mu })={\theta }({B})a_{t}  \]

or

\[  {\phi }({B})W_{t}=const+{\theta }({B})a_{t}  \]

where

\[  const = {\phi }({B}){\mu } = {\mu }-{\phi }_{1}{\mu }-{\phi }_{2}{\mu }-{\ldots }-{\phi }_{p}{\mu }  \]

Thus, when an autoregressive operator and a mean term are both included in the model, the constant term for the model can be represented as ${{\phi }({B}){\mu }}$. This value is printed with the label "Constant Estimate" in the ESTIMATE statement output.

Notation for Transfer Function Models

The general ARIMA model with input series, also called the ARIMAX model, is written as

\[  W_{t}={\mu } +\sum _{i}^{}{\frac{{\omega }_{i}({B})}{{\delta }_{i}({B}) }{B}^{k_{i}}X_{i,t}} +\frac{{\theta }({B})}{{\phi }({B})}a_{t}  \]

where

${X_{i,t}}$

is the ith input time series or a difference of the ith input series at time t

${k_{i}}$

is the pure time delay for the effect of the ith input series

${{\omega }_{i}({B})}$

is the numerator polynomial of the transfer function for the ith input series

${{\delta }_{i}({B})}$

is the denominator polynomial of the transfer function for the ith input series.

The model can also be written more compactly as

\[  W_{t}={\mu }+\sum _{i}^{}{{\Psi }_{i}({B})X_{i,t}}+n_{t}  \]

where

${{\Psi }_{\mi{i}}({B})}$

is the transfer function for the ith input series modeled as a ratio of the ${\omega }$ and ${\delta }$ polynomials: ${{\Psi }_{i}({B})=({\omega }_{i}({B})/{\delta }_{i}({B})){B}^{k_{i}}}$

${n_{t}}$

is the noise series: ${n_{t}=({\theta }({B})/{\phi }({B}))a_{t}}$

This model expresses the response series as a combination of past values of the random shocks and past values of other input series. The response series is also called the dependent series or output series. An input time series is also referred to as an independent series or a predictor series. Response variable, dependent variable, independent variable, or predictor variable are other terms often used.

Notation for Factored Models

ARIMA models are sometimes expressed in a factored form. This means that the ${\phi }$, ${\theta }$, ${\omega }$, or ${\delta }$ polynomials are expressed as products of simpler polynomials. For example, you could express the pure ARIMA model as

\[  W_{t}={\mu } +\frac{{\theta }_{1}({B}){\theta }_{2}({B})}{{\phi }_{1}({B}){\phi }_{2}({B})}a_{t}  \]

where ${{\phi }_{1}({B}){\phi }_{2}({B})={\phi }({B})}$ and ${{\theta }_{1}({B}){\theta }_{2}({B})={\theta }({B})}$.

When an ARIMA model is expressed in factored form, the order of the model is usually expressed by using a factored notation also. The order of an ARIMA model expressed as the product of two factors is denoted as ARIMA(p,d,q)${\times }$(P,D,Q).

Notation for Seasonal Models

ARIMA models for time series with regular seasonal fluctuations often use differencing operators and autoregressive and moving-average parameters at lags that are multiples of the length of the seasonal cycle. When all the terms in an ARIMA model factor refer to lags that are a multiple of a constant s, the constant is factored out and suffixed to the ARIMA(p,d,q ) notation.

Thus, the general notation for the order of a seasonal ARIMA model with both seasonal and nonseasonal factors is ARIMA(p,d,q)${\times }$(P,D,Q)$_{\mi{s}}$. The term (p,d,q) gives the order of the nonseasonal part of the ARIMA model; the term (P,D,Q)$_{\mi{s}}$ gives the order of the seasonal part. The value of s is the number of observations in a seasonal cycle: 12 for monthly series, 4 for quarterly series, 7 for daily series with day-of-week effects, and so forth.

For example, the notation ARIMA(0,1,2)${\times }$(0,1,1)$_{12}$ describes a seasonal ARIMA model for monthly data with the following mathematical form:

\[  (1-{B})(1-{B}^{12})Y_{t}={\mu } +(1-{\theta }_{1,1}{B}-{\theta }_{1,2}{B}^{2}) (1-{\theta }_{2,1}{B}^{12})a_{t}  \]