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 8-5<br />

The one-step predictors ˆYt<br />

for the airline passenger<br />

data (solid line) and the<br />

actual data (square boxes).<br />

8.6 State-Space Models with Missing Observations 283<br />

(thousands)<br />

100 200 300 400 500 600<br />

1949 1951 1953 1955 1957 1959 1961<br />

must be nearly constant and equal to ˆX12 2.171. The first three components of the<br />

predictors ˆXt are plotted in Figure 8.4. Notice that the first component varies like a<br />

random walk around a straight line, while the second component is nearly constant as<br />

a result of ˆσ 2 2 ≈ 0. The third component, corresponding to the seasonal component,<br />

exhibits a clear seasonal cycle that repeats roughly the same pattern throughout the<br />

12 years of data. The one-step predictors ˆXt1 + ˆXt3 of Yt are plotted in Figure 8.5<br />

(solid line) together with the actual data (square boxes). For this model the predictors<br />

follow the movement of the data quite well.<br />

8.6 State-Space Models with Missing Observations<br />

State-space representations and the associated Kalman recursions are ideally suited<br />

to the analysis of data with missing values, as was pointed out by Jones (1980) in the<br />

context of maximum likelihood estimation for ARMA processes. In this section we<br />

shall deal with two missing-value problems for state-space models. The first is the<br />

evaluation of the (Gaussian) likelihood based on {Yi1 ,...,Yir }, where i1,i2,...,ir<br />

are positive integers such that 1 ≤ i1 < i2 < ··· < ir ≤ n. (This allows for<br />

observation of the process {Yt} at irregular intervals, or equivalently for the possibility<br />

that (n−r)observations are missing from the sequence {Y1,...,Yn}.) The solution of<br />

this problem will, in particular, enable us to carry out maximum likelihood estimation<br />

for ARMA and ARIMA processes with missing values. The second problem to be<br />

considered is the minimum mean squared error estimation of the missing values<br />

themselves.

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

Saved successfully!

Ooh no, something went wrong!