Hi all,

I want to estimate a two-country DSGE model in Dynare but even with mh_replic=10 it takes forever (as in: since yesterday afternoon and still running…)

I use steady_state_model so that can’t be the issue. Do you have other ideas what could be the (or an obvious) issue? Any help is greatly appreciated.

Best,

Maren

nestedCES.mod (22.1 KB)

dataforestimation.mat (7.4 KB)

The steady state file is missing. Did mode-finding conclude? And which Dynare version are you using?

Hi,

I use the steady_state_model block in the .mod-file instead of a separate steady-state-file. Is that inefficient?

I’m using Dynare 4.4.3.

What exactly do you mean by ‘did mode-finding conclude’?

The estimation was ‘climbing the hill’ already but it just doesn’t seem to end

Hi,

I use the latest stable version now with mode_compute=1 but get the error message:

Error using chol

Matrix must be positive definite.

Error in posterior_sampler_initialization (line 84)

d = chol(vv);

Error in posterior_sampler (line 59)

[ ix2, ilogpo2, ModelName, MetropolisFolder, fblck, fline, npar, nblck, nruns,

NewFile, MAX_nruns, d, bayestopt_] = …

Error in dynare_estimation_1 (line 443)

posterior_sampler(objective_function,posterior_sampler_options.proposal_distribution,xparam1,posterior_sampler_options,bounds,dataset_,dataset_info,options_,M_,estim_params_,bayestopt_,oo

Error in dynare_estimation (line 105)

dynare_estimation_1(var_list,dname);

Error in nestedCES (line 795)

oo_recursive_=dynare_estimation(var_list_);

Error in dynare (line 223)

evalin(‘base’,fname) ;

The mode check plots also scream trouble -.-

Any further hints what could be the issue here?

ToT_N_steady_state.m (642 Bytes)

C_steady_state.m (351 Bytes)

C_star_steady_state.m (507 Bytes)

Thanks a lot Johannes,

yeah that was a stupid mistake with the observation equations. I removed the seasonalities in the inflation rates and let the estimation run. The key parameters of interest for me are alpha and gamma.

Dynare is not able to run the estimation with mode_compute=1 (hessian is not positive) so I used mode_compute=6 to check whether this one works at least. The procedure found the mode and gave me the attached results (i included the plots in the .zip file). I’m pretty puzzled by how tight the posterior distributions are. What could be an obvious issue here? These results are clearly fishy.

Also: I am not so sure how to treat R_obs in the data. We have the downward trend due to the crisis in 2008. Would I have to remove this trend?

Sorry for any inconveniences!

<a class=“attachment” href="/uploads/default/original/2X/a/a89028ca6583865452358ab9c9estimation plots.zip (115.9 KB)

a998b354edb762.m">C_star_steady_state.m (507 Bytes)

ToT_Nstar_steady_state.m (578 Bytes)

C_steady_state.m (351 Bytes)

nestedCES.mod (22.1 KB)

dataforestimation.mat (7.4 KB)