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.

20 3 L<strong>in</strong>ear Gaussian <strong>state</strong> <strong>space</strong> <strong>models</strong> and the Kalman filter<br />

where a1 and P 1 are assumed to be known for now. The <strong>in</strong>itialization <strong>of</strong> the <strong>state</strong> vector<br />

will be discussed <strong>in</strong> ¢ 3.3.6.<br />

The matrices T t, ct, Rt, Zt, dt, Q t and Ht are called system matrices. Unless<br />

<strong>state</strong>d otherwise, they are assumed to be non-stochastic, i.e. they change over time <strong>in</strong><br />

a predeterm<strong>in</strong>ed way. This leads to a system <strong>in</strong> which y t can be expressed as a l<strong>in</strong>ear<br />

comb<strong>in</strong>ation <strong>of</strong> ξ 1 and past and present <strong>state</strong> and observation disturbances. A <strong>state</strong> <strong>space</strong><br />

model with constant system matrices over time is referred to as be<strong>in</strong>g time-homogeneous.<br />

Usually, at least some <strong>of</strong> the elements <strong>of</strong> the system matrices T t and Rt <strong>in</strong> the <strong>state</strong><br />

equation and Zt, Q t and Ht <strong>in</strong> the observation equation depend on a vector ψ <strong>of</strong><br />

unknown parameters, which are referred to as hyperparameters. While the hyperparameters<br />

control for the stochastic characteristics <strong>of</strong> the model, the parameters related to<br />

the vectors ct and dt have a determ<strong>in</strong>istic impact on the observations and the expected<br />

<strong>state</strong> value. The values <strong>of</strong> the system matrices are typically unknown. Their estimation<br />

on the basis <strong>of</strong> the available observations will be the subject <strong>of</strong> ¢ 3.4. For now, the system<br />

matrices are assumed to be known.<br />

Throughout this thesis, the s<strong>of</strong>tware package SsfPack 2.3 by Koopman et al. (1999)<br />

is employed to apply the Kalman filter to estimate <strong>models</strong> <strong>in</strong> <strong>state</strong> <strong>space</strong> form. SsfPack is<br />

a comprehensive collection <strong>of</strong> C rout<strong>in</strong>es for estimat<strong>in</strong>g <strong>models</strong> <strong>in</strong> <strong>state</strong> <strong>space</strong> form with<br />

general rout<strong>in</strong>es for smooth<strong>in</strong>g, filter<strong>in</strong>g, likelihood evaluation and forecast<strong>in</strong>g l<strong>in</strong>ked<br />

to the comput<strong>in</strong>g environment Ox 3.30 by Doornik (2001). With regard to the use <strong>of</strong><br />

SsfPack, it will prove useful to represent the <strong>state</strong> <strong>space</strong> model <strong>in</strong>troduced by (3.1)–(3.7)<br />

<strong>in</strong> compact form as<br />

�<br />

ξt+1<br />

yt �<br />

= Φtξt + δt + υt, υt ∼ N(0, Ωt), t = 1, . . . , T, (3.8)<br />

with<br />

�<br />

T t<br />

Φt =<br />

Zt<br />

� �<br />

ct<br />

, δt =<br />

dt<br />

� �<br />

Rtηt , υt =<br />

ɛt<br />

� �<br />

Qt 0<br />

, Ωt =<br />

0 Ht<br />

�<br />

, (3.9)<br />

where Φt is <strong>of</strong> dimension (m + N) × m, δt is (m + N) × 1, υt is (m + N) × 1 and Ωt is<br />

(m + N) × (m + N). The <strong>in</strong>itial values are collected <strong>in</strong> the (m + 1) × m matrix Σ which<br />

is def<strong>in</strong>ed as<br />

�<br />

P 1<br />

Σ =<br />

a ′ 1<br />

�<br />

. (3.10)<br />

For a survey on us<strong>in</strong>g SsfPack for <strong>state</strong> <strong>space</strong> model<strong>in</strong>g <strong>in</strong> macroeconomics and f<strong>in</strong>ance,<br />

see Zivot et al. (2002).<br />

3.3 The Kalman filter and smoother<br />

Once a model is put <strong>in</strong>to <strong>state</strong> <strong>space</strong> form, the Kalman filter can be employed to compute<br />

optimal forecasts <strong>of</strong> the mean and covariance matrix <strong>of</strong> the normally distributed <strong>state</strong><br />

vector ξ t+1, based on the available <strong>in</strong>formation through time t. This section outl<strong>in</strong>es<br />

how the Kalman filter can be used for estimation by filter<strong>in</strong>g and smooth<strong>in</strong>g: filter<strong>in</strong>g<br />

bases an <strong>in</strong>ference about the <strong>state</strong> vector only on the <strong>in</strong>formation up to time t; smooth<strong>in</strong>g

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

Saved successfully!

Ooh no, something went wrong!