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.

5.1 Preliminary Estimation 149<br />

Example 5.1.4 The lake data<br />

estimated white noise variance and larger Gaussian likelihood. From the ratio of the<br />

estimated coefficient to (1.96× standard error) displayed by ITSM, we obtain the<br />

95% confidence bounds for φ: 0.4371 ± 0.4371/2.1668 (0.2354, 0.6388).<br />

This series {Yt,t 1,...,98} has already been studied in Example 1.3.5. In this<br />

example we shall consider the problem of fitting an AR process directly to the data<br />

without first removing any trend component. A graph of the data was displayed in<br />

Figure 1.9. The sample ACF and PACF are shown in Figures 5.3 and 5.4, respectively.<br />

The sample PACF shown in Figure 5.4 strongly suggests fitting an AR(2) model<br />

to the mean-corrected data Xt Yt − 9.0041. After clicking on the blue preliminary<br />

estimation button of ITSM select Yes to subtract the sample mean from {Yt}. Then<br />

specify 2 for the AR order, 0 for the MA order, and Burg for estimation. Click OK<br />

to obtain the model<br />

Xt − 1.0449Xt−1 + 0.2456Xt−2 Zt, {Zt} ∼WN(0, 0.4706),<br />

with AICC value 213.55 and 95% confidence bounds<br />

φ1 :1.0449 ± 1.0449/5.5295 (0.8559, 1.2339),<br />

φ2 : −0.2456 ± 0.2456/1.2997 (−0.4346, −0.0566).<br />

Selecting the Yule–Walker method for estimation, we obtain the model<br />

ACF<br />

-0.2 0.0 0.2 0.4 0.6 0.8 1.0<br />

Figure 5-3<br />

The sample ACF of the lake<br />

data in Example 5.1.4. Lag<br />

Xt − 1.0538Xt−1 + 0.2668Xt−2 Zt, {Zt} ∼WN(0, 0.4920),<br />

0 10 20 30 40

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

Saved successfully!

Ooh no, something went wrong!