Page 1 of 1
mode check-shock variance log likelihood kernel is flat
Posted:
Thu Sep 12, 2013 7:08 pm
by bigbig001
Hi,
I am working on DSGE Bayesian estimation. I found out that the green line (log likelihood kernel) in mode check for all the shock variances are always flat (i tried different versions of model specification but I always have this problem). I was wondering if this is a big problem and how to solve it. In MCMC diagnostic graph, these variances tend to converge very slowly compared to other parameters.
I appreciate any help on it.
Re: mode check-shock variance log likelihood kernel is flat
Posted:
Sat Sep 14, 2013 8:18 am
by jpfeifer
For estimation, this is a bad sign. Your variances are not identified from the data. Thus, your prior will be your posterior. Try using the identification command to see if this is really the case.
Re: mode check-shock variance log likelihood kernel is flat
Posted:
Sun Sep 15, 2013 2:45 am
by bigbig001
Thanks a lot, jpfeifer.
I did the identification test but the program stopped in the sensitivity analysis. The message dynare gave me is:
STEADY: numerical initial values or parameters incompatible with the following equations
1 2 3 4 13 14 20 30
Please check for example
i) if all parameters occurring in these equations are defined
ii) that no division by an endogenous variable initialized to 0 occurs
All parameter values in the specified ranges are not acceptable!
I calculated the steady state by hand and set the initial values equal to the steady state. I use resid(1) to check the initial values and the residuals of all equations are 0. I think all the initial values are correct steady states. I am confused about this message.
I was wondering what is wrong with model or estimation.
Thanks a lot
Re: mode check-shock variance log likelihood kernel is flat
Posted:
Sun Sep 15, 2013 8:11 am
by jpfeifer
Please provide the mod-file
Re: mode check-shock variance log likelihood kernel is flat
Posted:
Sun Sep 15, 2013 3:55 pm
by bigbig001
Thanks a lot, jpfeifer
Here is my mod file and data file. My model is similar to Iacoveillo (2010), but much more simplified. And I am using demeaned log difference of GDP,consumption, investment and housing price, which are stationary.
Re: mode check-shock variance log likelihood kernel is flat
Posted:
Mon Sep 16, 2013 12:36 pm
by jpfeifer
This is a quasi-bug in Dynare. Your combination of using a parameter defined within the shocks-block to estimate a variance does not work. Instead of
- Code: Select all
shocks;
var ec;sicg^2;
end;
estimated_params;
sigc, inv_gamma_PDF,0.05,0.1;
end;
use
- Code: Select all
shocks;
var ec;sicg^2;
end;
estimated_params;
stderr ec, inv_gamma_PDF,0.05,0.1;
end;
Re: mode check-shock variance log likelihood kernel is flat
Posted:
Tue Sep 17, 2013 12:00 am
by bigbig001
Thanks a lot, Jpfeifer.
The standard deviation of shocks are identified after I modified my code. I really appreciate your help. I have been struggling on this problem for a month and tried many different model specifications. I was desperate. Thanks very much.
But I still get the same error message when I run the sensitivity analysis:
STEADY: numerical initial values or parameters incompatible with the following equations
1 2 3 4 13 14 20 30
Please check for example
i) if all parameters occurring in these equations are defined
ii) that no division by an endogenous variable initialized to 0 occurs
All parameter values in the specified ranges are not acceptable!
The initial values were the actual steady states which were calculated by hand and the the residues of all equations are 0 with resid(1). I am confused about this message.
Re: mode check-shock variance log likelihood kernel is flat
Posted:
Tue Sep 17, 2013 8:24 am
by jpfeifer
Your computed steady state is for the calibrated model. But it is not correctly updated for the estimated values. You might need a steady state file to do this.
Re: mode check-shock variance log likelihood kernel is flat
Posted:
Sat Sep 21, 2013 2:43 pm
by bigbig001
Thanks a lot, jpfeifer.
I will write a separate steady state file and try again.