02.12.2012 Views

Applications of state space models in finance

Applications of state space models in finance

Applications of state space models in finance

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.

3.4 Maximum likelihood estimation 27<br />

elements <strong>of</strong> a multivariate series. As it is straightforward to deal with this problem for<br />

univariate series, Durb<strong>in</strong> and Koopman (2001, ¢ 5.1) propose to br<strong>in</strong>g <strong>in</strong> the components<br />

<strong>of</strong> a multivariate series <strong>in</strong>to the analysis one at a time. Thus, a multivariate series is<br />

converted <strong>in</strong>to a univariate series. This leads to computational ga<strong>in</strong>s and significantly<br />

simplifies the process <strong>of</strong> diffuse <strong>in</strong>itialization for multivariate series. For a detailed<br />

discussion <strong>of</strong> the univariate treatment <strong>of</strong> multivariate series, the reader is referred to<br />

Durb<strong>in</strong> and Koopman (2001, ¢ 6.4).<br />

3.3.7 The Kalman filter with non-Gaussian errors<br />

The derivation <strong>of</strong> the Kalman filter presented above and the related estimation procedures<br />

are based on the assumption <strong>of</strong> normally distributed disturbances. It was shown<br />

how the conditional distribution <strong>of</strong> the <strong>state</strong> vector ξ t can be calculated recursively given<br />

the <strong>in</strong>formation set at time t, for all t = 1, . . . , T . As these conditional distributions are<br />

also normal, they are fully specified by their first two moments. These can be computed<br />

by the Kalman filter. The conditional mean <strong>of</strong> the <strong>state</strong> vector represents an optimal<br />

estimator <strong>in</strong> the sense that is has m<strong>in</strong>imum MSE matrix.<br />

In case <strong>of</strong> non-Gaussian disturbances, the Kalman filter is no longer guaranteed to<br />

yield the conditional mean <strong>of</strong> the <strong>state</strong> vector. However, it nevertheless represents an optimal<br />

estimator <strong>in</strong> the sense that no other l<strong>in</strong>ear estimator has a smaller MSE. Therefore,<br />

the Kalman filter can still be employed when the normality assumption is dropped. The<br />

estimators obta<strong>in</strong>ed by maximiz<strong>in</strong>g the Gaussian likelihood function with observations<br />

that are not normally distributed are referred to as quasi-maximum likelihood (QML)<br />

estimators; cf. Hamilton (1994a). 6<br />

3.4 Maximum likelihood estimation<br />

So far, the system matrices have been assumed to be known. In the more general<br />

case, they depend at least partially on ψ, the vector <strong>of</strong> unknown parameters. This<br />

section expla<strong>in</strong>s how the vector <strong>of</strong> unknowns can be estimated by means <strong>of</strong> maximum<br />

likelihood (ML). After briefly summariz<strong>in</strong>g the general idea beh<strong>in</strong>d the concept <strong>of</strong> ML,<br />

Subsection 3.4.1 <strong>in</strong>troduces the loglikelihood function for the general Gaussian <strong>state</strong><br />

<strong>space</strong> model. It will be demonstrated how the general model can be reparameterized<br />

<strong>in</strong> order to reduce the dimension <strong>of</strong> the vector <strong>of</strong> parameters by one. Subsections 3.4.2<br />

and 3.4.3 briefly overview the two alternative concepts <strong>of</strong> maximiz<strong>in</strong>g the loglikelihood:<br />

direct numerical maximization and the EM algorithm. F<strong>in</strong>ally, Subsection 3.4.4 shows<br />

how parameter restrictions can be implemented.<br />

3.4.1 The loglikelihood function<br />

In order to estimate a model by ML, it has to be parametric and fully specified through<br />

the jo<strong>in</strong>t probability density function; the parameter values have to conta<strong>in</strong> all the<br />

necessary <strong>in</strong>formation for a simulation <strong>of</strong> the dependent variables. For the T sets <strong>of</strong><br />

6 For an <strong>in</strong>troduction to <strong>state</strong> <strong>space</strong> <strong>models</strong> that take non-Gaussianity explicitly <strong>in</strong>to account,<br />

which is beyond the scope <strong>of</strong> this thesis, see, for example, Durb<strong>in</strong> and Koopman (2001, Part II).

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

Saved successfully!

Ooh no, something went wrong!