Page 1 of 1

IRF of output to tech. shock in Lubik and Schorfheide (2007)

PostPosted: Fri Jul 03, 2015 6:04 am
by earsmall
Hi, owing to your support, I have made some progress in reproducing Lubik and Schorfheide (2007).

Re: IRF of output to tech. shock in Lubik and Schorfheide (2

PostPosted: Mon Jul 06, 2015 5:37 am
by jpfeifer
The reason is that TFP has a unit root, but the model is entered in detrended form. To get the permanent IRFs to TFP shock, you need to reconstruct the underlying non-stationary series. An example for this can be found in the AguiarGopinath2007.mod on my homepage at https://sites.google.com/site/pfeiferecon/dynare

Re: IRF of output to tech. shock in Lubik and Schorfheide (2

PostPosted: Mon Jul 06, 2015 6:58 am
by earsmall
Thank you very much.
I will try with reference to this mod file !
:lol:

Re: IRF of output to tech. shock in Lubik and Schorfheide (2

PostPosted: Mon Jul 06, 2015 10:54 am
by earsmall
jpfeifer wrote:The reason is that TFP has a unit root, but the model is entered in detrended form. To get the permanent IRFs to TFP shock, you need to reconstruct the underlying non-stationary series. An example for this can be found in the AguiarGopinath2007.mod on my homepage at https://sites.google.com/site/pfeiferecon/dynare



There are three things in my modified code compared to original code (detrended form).
Code: Select all
model(linear);
...
...
z = a - a(-1) ;  //z_t is growth rate of technology (A)

yy = y + a ; // y_t = Y_t/A_t


However, the response from dynare is that "unstable".

I searched a lot to solve this problem and concluded that I should use trend_var command or .
Unfortunately, these do not work as well.

And finally, dynare works when I add diffuse_filter command in estimation()
and I got the IRF that I wanted.

But I am not sure whether I solved well.

Please give me more comments.


-Jerry

Re: IRF of output to tech. shock in Lubik and Schorfheide (2

PostPosted: Mon Jul 06, 2015 3:01 pm
by jpfeifer
You are introducing a unit root in a. To deal with this non-stationarity in estimation, you need the diffuse filter or lik_init=2.