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.

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

1.5.2 Estimation and Elimination of Both Trend and Seasonality<br />

The methods described for the estimation and elimination of trend can be adapted in<br />

a natural way to eliminate both trend and seasonality in the general model, specified<br />

as follows.<br />

Classical Decomposition Model<br />

Xt mt + st + Yt, t 1,...,n, (1.5.11)<br />

where EYt 0, st+d st, and d<br />

j1 sj 0.<br />

We shall illustrate these methods with reference to the accidental deaths data of<br />

Example 1.1.3, for which the period d of the seasonal component is clearly 12.<br />

Method S1: Estimation of Trend and Seasonal Components<br />

The method we are about to describe is used in the Transform>Classical option<br />

of ITSM.<br />

Suppose we have observations {x1,...,xn}. The trend is first estimated by applying<br />

a moving average filter specially chosen to eliminate the seasonal component<br />

and to dampen the noise. If the period d is even, say d 2q, then we use<br />

ˆmt (0.5xt−q + xt−q+1 +···+xt+q−1 + 0.5xt+q)/d, q < t ≤ n − q. (1.5.12)<br />

If the period is odd, say d 2q + 1, then we use the simple moving average (1.5.5).<br />

The second step is to estimate the seasonal component. For each k 1,...,d,we<br />

compute the average wk of the deviations {(xk+jd−ˆmk+jd), q < k+jd ≤ n−q}. Since<br />

these average deviations do not necessarily sum to zero, we estimate the seasonal<br />

component sk as<br />

ˆsk wk − d −1<br />

d<br />

wi, k 1,...,d, (1.5.13)<br />

i1<br />

and ˆsk ˆsk−d,k >d.<br />

The deseasonalized data is then defined to be the original series with the estimated<br />

seasonal component removed, i.e.,<br />

dt xt −ˆst, t 1,...,n. (1.5.14)<br />

Finally, we reestimate the trend from the deseasonalized data {dt} using one of<br />

the methods already described. The program ITSM allows you to fit a least squares<br />

polynomial trend ˆm to the deseasonalized series. In terms of this reestimated trend<br />

and the estimated seasonal component, the estimated noise series is then given by<br />

ˆYt xt −ˆmt −ˆst, t 1,...,n.

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

Saved successfully!

Ooh no, something went wrong!