Chapter Contents
Chapter Contents
Previous
Previous
Next
Next
Forecasting Process Details

Notation for ARIMA Models

A dependent time series that is modeled as a linear combination of its own past values and past values of an error series is known as a (pure) ARIMA model.

Nonseasonal ARIMA Model Notation

The order of an ARIMA 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 (rarely should d > 2 be needed)

q
is the order of the moving-average process

Given a dependent time series {\{Y_{t} : 1 \le t \le n \}},mathematically the ARIMA model is written as

(1-B)^dY_{t}
={\mu}+\frac{{\theta}(B)}{{\phi}(B) }a_{t}

where
t
indexes time

{\mu}
is the mean term

B
is the backshift operator; that is, BXt=Xt-1

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

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

at
is the independent disturbance, also called the random error

For example, the mathematical form of the ARIMA(1,1,2) model is

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

Seasonal ARIMA Model Notation

Seasonal ARIMA models are expressed in factored form by the notation ARIMA(p,d,q)(P,D,Q)s, where

P
is the order of the seasonal autoregressive part

D
is the order of the seasonal differencing (rarely should D > 1 be needed)

Q
is the order of the seasonal moving-average process

s
is the length of the seasonal cycle

Given a dependent time series {\{Y_{t} : 1 \le t \le n \}},mathematically the ARIMA seasonal model is written as

(1-B)^d(1-B^s)^DY_{t}
={\mu}+\frac{{\theta}(B){\theta}_{s}(B^s)}{{\phi}(B) {\phi}_{s}(B^s) }a_{t}

where
{{\phi}_{s}(B^s)}
is the seasonal autoregressive operator, represented as a polynomial in the back shift operator:
{{\phi}_{s}(B^s)=1-{\phi}_{s,1}B^s -{ ... }- {\phi}_{s,P}B^{sP}}
{{\theta}_{s}(B^s)}
is the seasonal moving-average operator, represented as a polynomial in the back shift operator: {{\theta}_{s}(B^s)=1-{\theta}_{s,1}B^s -{ ... }- {\theta}_{s,Q}B^{sQ}}

For example, the mathematical form of the ARIMA(1,0,1)(1,1,2)12 model is

(1-B^{12})Y_{t}
={\mu}+\frac{(1-{\theta}_{1}B)(1-{\theta}_{s,1}B^{12}
 -{\theta}_{s,2}B^{24})}{(1-{\phi} _{1}B)(1-{\phi} _{s,1}B^{12}) }
 a_{t}

Abbreviated Notation for ARIMA Models

If the differencing order, autoregressive order, or moving-average order is zero, the notation is further abbreviated as

I(d)(D)s
integrated model or ARIMA(0,d,0)(0,D,0)

AR(p)(P)s
autoregressive model or ARIMA(p,0,0)(P,0,0)

IAR(p,d)(P,D)s
integrated autoregressive model or ARIMA(p,d,0)(P,D,0)s

MA(q)(Q)s
moving average model or ARIMA(0,0,q)(0,0,Q)s

IMA(d,q)(D,Q)s
integrated moving average model or ARIMA(0,d,q)(0,D,Q)s

ARMA(p,q)(P,Q)s
autoregressive moving-average model or ARIMA(p,0,q)(P,0,Q)s

Notation for Transfer Functions

A transfer function can be used to filter a predictor time series to form a dynamic regression model.

Let Yt be the dependent series and let Xt be the predictor series, and let {{\Psi}(B)} be a linear filter or transfer function for the effect of Xt on Yt. The ARIMA model is then

(1-B)^d(1-B^s)^DY_{t}
= {\mu} + {\Psi}(B)(1-B)^d(1-B^s)^DX_{t}
+ \frac{{\theta}(B){\theta}_{s}(B^s)}{{\phi}(B){\phi}_{s}(B^s)}a_{t}

This model is called a dynamic regression of Yt on Xt.

Nonseasonal Transfer Function Notation



Given the ith predictor time series {\{X_{i,t} : 1 \le t \le n \}},the transfer function is written as [Dif(di)Lag(ki)N(qi)/ D(pi)] where

di
is the simple order of the differencing for the ith predictor time series, (1-B)diXi,t (rarely should di > 2 be needed)

ki
is the pure time delay (lag) for the effect of the ith predictor time series, Xi,tBki = Xi,t-ki

pi
is the simple order of the denominator for the ith predictor time series

qi
is the simple order of the numerator for the ith predictor time series

The mathematical notation used to describe a transfer function is

{\Psi}_{i}(B)=
\frac{{\omega}_{i}(B)}{{\delta}_{i}(B) }
(1-B)^{d_{i}}B^{k_{i}}

where
B
is the backshift operator; that is, BXt=Xt-1

{{\delta}_{i}(B)}
is the denominator polynomial of the transfer function for the ith predictor time series: {{\delta}_{i}(B)=1-{\delta}_{i,1}B-{ ... }- {\delta}_{i,p_{i}}B^{p_{i}}}

{{\omega}_{i}(B)}
is the numerator polynomial of the transfer function for the ith predictor time series: {{\omega}_{i}(B)=1-{\omega}_{i,1}B-{ ... }- {\omega}_{i,q_{i}}B^{q_{i}}}

The numerator factors for a transfer function for a predictor series are like the MA part of the ARMA model for the noise series. The denominator factors for a transfer function for a predictor series are like the AR part of the ARMA model for the noise series. Denominator factors introduce exponentially weighted, infinite distributed lags into the transfer function.

For example, the transfer function for the ith predictor time series with

ki=3
time lag is 3

di=1
simple order of differencing is one

pi=1
simple order of the denominator is one

qi=2
simple order of the numerator is two

would be written as [Dif(1)Lag(3)N(2)/D(1)]. The mathematical notation for the transfer function in this example is

{\Psi}_{i}(B)=
{(1-{\omega}_{i,1}B-{\omega}_{i,2}B^2)
 \over (1-{\delta}_{i,1}B)}
(1-B)B^3

Seasonal Transfer Function Notation



The general transfer function notation for the ith predictor time series Xi,t with seasonal factors is [Dif(di)(Di)s Lag(ki) N(qi)(Qi)s/ D(pi)(Pi)s] where

Di
is the seasonal order of the differencing for the ith predictor time series (rarely should Di > 1 be needed)

Pi
is the seasonal order of the denominator for the ith predictor time series (rarely should Pi > 2 be needed)

Qi
is the seasonal order of the numerator for the ith predictor time series, (rarely should Qi > 2 be needed)

s
is the length of the seasonal cycle

The mathematical notation used to describe a seasonal transfer function is

{\Psi}_{i}(B)=
 \frac{{\omega}_{i}(B){\omega}_{s,i}(B^s)}{{\delta}_{i}(B){\delta}_{s,i}(B^s) }
 (1-B)^{d_{i}}(1-B^s)^{D_{i}}B^{k_{i}}

where

{{\delta}_{s,i}(B^s)}
is the denominator seasonal polynomial of the transfer function for the ith predictor time series:
{{\delta}_{s,i}(B)=1-{\delta}_{s,i,1}B-{ ... }- {\delta}_{s,i,P_{i}}B^{sP_{i}}}

{{\omega}_{s,i}(B^s)}
is the numerator seasonal polynomial of the transfer function for the ith predictor time series:
{{\omega}_{s,i}(B)=1-{\omega}_{s,i,1}B-{ ... }- {\omega}_{s,i,Q_{i}}B^{sQ_{i}}}

For example, the transfer function for the ith predictor time series Xi,t whose seasonal cycle s=12 with

di=2
simple order of differencing is two

Di = 1
seasonal order of differencing is one

qi=2
simple order of the numerator is two

Qi = 1
seasonal order of the numerator is one

would be written as [Dif(2)(1)s N(2)(1)s]. The mathematical notation for the transfer function in this example is

{\Psi}_{i}(B)=
(1-{\omega}_{i,1}B-{\omega}_{i,2}B^2)
(1-{\omega}_{s,i,1}B^{12})
(1-B)^2(1-B^{12})

Note: In this case, [Dif(2)(1)s N(2)(1)s] = [Dif(2)(1)sLag(0)N(2)(1)s/D(0)(0)s].

Chapter Contents
Chapter Contents
Previous
Previous
Next
Next
Top
Top

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