04.01.2013 Views

Springer - Read

Springer - Read

Springer - Read

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Figure 1-22<br />

Strike data smoothed<br />

by elimination of high<br />

frequencies with f 0.4.<br />

1.5 Estimation and Elimination of Trend and Seasonal Components 29<br />

(thousands)<br />

3.5 4.0 4.5 5.0 5.5 6.0<br />

1950 1960 1970 1980<br />

a0,a1, and a2 to minimize the sum of squares, n<br />

t1 (xt − mt) 2 (see Example 1.3.4).<br />

The method of least squares estimation can also be used to estimate higher-order<br />

polynomial trends in the same way. The Regression option of ITSM allows least<br />

squares fitting of polynomial trends of order up to 10 (together with up to four harmonic<br />

terms; see Example 1.3.6). It also allows generalized least squares estimation<br />

(see Section 6.6), in which correlation between the residuals is taken into account.<br />

Method 2: Trend Elimination by Differencing<br />

Instead of attempting to remove the noise by smoothing as in Method 1, we now<br />

attempt to eliminate the trend term by differencing. We define the lag-1 difference<br />

operator ∇ by<br />

∇Xt Xt − Xt−1 (1 − B)Xt, (1.5.9)<br />

where B is the backward shift operator,<br />

BXt Xt−1. (1.5.10)<br />

Powers of the operators B and ∇ are defined in the obvious way, i.e., B j (Xt) Xt−j<br />

and ∇ j (Xt) ∇(∇ j−1 (Xt)), j ≥ 1, with ∇ 0 (Xt) Xt. Polynomials in B and ∇ are<br />

manipulated in precisely the same way as polynomial functions of real variables. For<br />

example,<br />

∇ 2 Xt ∇(∇(Xt)) (1 − B)(1 − B)Xt (1 − 2B + B 2 )Xt<br />

Xt − 2Xt−1 + Xt−2.

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

Saved successfully!

Ooh no, something went wrong!