I have been trying to estimate a model using dynare.

It runs well until computing Bayesian IRFs but shows errors in computing Bayesian variance decomposition.

[code]

Subscript indices must either be real

positive integers or logicals.

Error in kernel_density_estimate (line 82)

xi(indx+[1 2]) = xi(indx+[1 2]) +

[1-temp,temp]’;

Error in posterior_moments (line 84)

[density(:,1),density(:,2)] =

kernel_density_estimate(xx,number_of_grid_points,…

Error in correlation_mc_analysis (line 95)

posterior_moments(tmp,1,mh_conf_sig);

Error in posterior_analysis>job (line 71)

oo_ =

correlation_mc_analysis(SampleSize,‘posterior’,M_.dname,M_.fname,…

Error in posterior_analysis (line 38)

oo_ =

job(type,SampleSize,arg1,arg2,arg3,options_,M_,oo_,nvar,vartan);

Error in compute_moments_varendo (line 83)

oo_ =

posterior_analysis(‘correlation’,var_list_(i,:),var_list_(j,:),h,options_,M_,oo_);

Error in dynare_estimation_1 (line 818)

oo_ =

compute_moments_varendo(‘posterior’,options_,M_,oo_,var_list_);

Error in dynare_estimation (line 89)

dynare_estimation_1(var_list,dname);

Error in est_3edy_post (line 484)

dynare_estimation(var_list_);

Error in dynare (line 180)

evalin(‘base’,fname) ;[/code]

I hope to know what makes these errors.

I attached my mod-file and data.