04.01.2013 Views

Springer - Read

Springer - Read

Springer - Read

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

102 Chapter 3 ARMA Models<br />

and<br />

rn → 1.<br />

Algebraic calculation of the coefficients θnj and rn is not feasible except for very simple<br />

models, such as those considered in the following examples. However, numerical<br />

implementation of the recursions is quite straightforward and is used to compute<br />

predictors in the program ITSM.<br />

Example 3.3.1 Prediction of an AR(p) process<br />

Applying (3.3.7) to the ARMA(p, 1) process with θ1 0, we easily find that<br />

ˆXn+1 φ1Xn +···+φpXn+1−p, n ≥ p.<br />

Example 3.3.2 Prediction of an MA(q) process<br />

Applying (3.3.7) to the ARMA(1,q) process with φ1 0gives<br />

ˆXn+1 <br />

min(n,q) <br />

j1<br />

θnj<br />

<br />

Xn+1−j − ˆXn+1−j<br />

<br />

, n ≥ 1,<br />

where the coefficients θnj are found by applying the innovations algorithm to the covariances<br />

κ(i,j) defined in (3.3.3). Since in this case the processes {Xt} and {σ −1 Wt}<br />

are identical, these covariances are simply<br />

κ(i,j) σ −2 γX(i − j) <br />

q−|i−j| <br />

Example 3.3.3 Prediction of an ARMA(1,1) process<br />

If<br />

r0<br />

θrθr+|i−j|.<br />

Xt − φXt−1 Zt + θZt−1, {Zt} ∼WN 0,σ 2 ,<br />

and |φ| < 1, then equations (3.3.7) reduce to the single equation<br />

ˆXn+1 φXn + θn1(Xn − ˆXn), n ≥ 1.<br />

To compute θn1 we first use Example 3.2.1 to find that γX(0)σ 2 1 + 2θφ + θ 2 / 1−<br />

φ2 . Substituting in (3.3.3) then gives, for i, j ≥ 1,<br />

⎧ 2<br />

1 + 2θφ + θ / 1 − φ 2 , i j 1,<br />

⎪⎨ 1 + θ<br />

κ(i,j) <br />

⎪⎩<br />

2 , i j ≥ 2,<br />

θ, |i − j| 1,i ≥ 1,<br />

0, otherwise.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!