next up previous
Next: About this document

Econ 551: Extra Credit Problem on Markov Processes and Invariant Densities
Professor John Rust, Department of Economics, Yale University

Let tex2html_wrap_inline67 be a first order Gaussian process defined by the linear AR(1) system:

displaymath69

where tex2html_wrap_inline71 is and IID Gaussian process with mean 0 and variance tex2html_wrap_inline75 (i.e. tex2html_wrap_inline77 and tex2html_wrap_inline79 are mutually independent whenever tex2html_wrap_inline81 and for each t we have tex2html_wrap_inline85 ). Assume that for all t tex2html_wrap_inline89 is independent of tex2html_wrap_inline77 .

1.
Is tex2html_wrap_inline67 a Markov process? If so, write down a formula for its transition density tex2html_wrap_inline95 .

2.
Under what conditions on tex2html_wrap_inline97 is tex2html_wrap_inline67 a stationary, ergodic Markov process?

3.
What is the recurrent class for tex2html_wrap_inline67 ? (The recurrent class of a Markov process tex2html_wrap_inline67 is defined as the smallest closed set C such that tex2html_wrap_inline107 , where i.o. denotes ``infinitely often''.)

4.
Assuming that tex2html_wrap_inline67 is stationary and ergodic, show that its invariant density p(y) is:

displaymath113

Hint: Equation (2) states that the invariant density p(y) is normally distributed with mean tex2html_wrap_inline117 and variance tex2html_wrap_inline119 . Clearly these quantities will make sense only if tex2html_wrap_inline121 , which should give you a pretty good idea about how to answer question 2 above. To see why the invariant density should be a normal distribution, use lag operators and write equation (1) as:

displaymath123

where tex2html_wrap_inline125 is the identity operator and tex2html_wrap_inline127 is the lag operator. If tex2html_wrap_inline121 then the operator tex2html_wrap_inline131 is invertible and we can write

displaymath133

By analogy to geometric series the inverse of tex2html_wrap_inline131 is given by an infinite geometric series:

displaymath137

Since tex2html_wrap_inline71 is an IID mean zero Gaussian process, it is clear from equation (5) that tex2html_wrap_inline67 is also Gaussian process, but with mean tex2html_wrap_inline117 and variance given by:

displaymath145

Since the time index t+1 is arbitrary, it follows from (5) and (6) that tex2html_wrap_inline89 normally distributed with mean tex2html_wrap_inline117 and variance tex2html_wrap_inline119 for all t. This suggests that tex2html_wrap_inline89 is in ``stochastic steady state'' and that the tex2html_wrap_inline159 distribution for tex2html_wrap_inline89 is indeed the invariant distribution of this process. However to be rigorous, we need to verify that p(y) defined in equation (2) satisfies the definition of an invariant density:

displaymath165

This requires some tedious calculations, but you should be able to verify that p(y) in equation (2) does indeed satisfy equation (7) (in doing the integral in (7) you may need to make use of the fact that a moment generating function for tex2html_wrap_inline169 is given by tex2html_wrap_inline171 ).

Another way to guess that equation (2) gives the invariant density for this case is to assume that stationarity holds, in which case

displaymath173

for all s and t. Similarly

displaymath179

for all s and t. Using these results and the fact that tex2html_wrap_inline85 we have:

displaymath187

and

displaymath189

So stochastic stationarity implies that tex2html_wrap_inline89 has mean tex2html_wrap_inline193 and variance tex2html_wrap_inline195 . The fact that tex2html_wrap_inline67 is Gaussian follows from the fact that tex2html_wrap_inline89 is equal to a sum of normally distributed random variables as you can see from the moving average representation for tex2html_wrap_inline89 given in equation (5).




next up previous
Next: About this document

John Rust
Fri Feb 7 16:00:33 CST 1997