when I am trying to replicate the Bayesian estimation with the example of RBC_Est which is done in the User Guide version beta (chapter 5). I get the error information as following :
??? Error using ==> chol
Matrix must be positive definite.
Error in ==> metropolis_hastings_initialization at 52
d = chol(vv);
Error in ==> random_walk_metropolis_hastings at 43
ix2, ilogpo2, ModelName, MhDirectoryName, fblck, fline, npar,
nblck, nruns, NewFile, MAX_nruns, d ] = …
Error in ==> dynare_estimation_1 at 940
feval(options_.posterior_sampling_method,‘DsgeLikelihood’,options_.proposal_distribution,xparam1,invhess,bounds,gend,data);
Error in ==> dynare_estimation at 62
dynare_estimation_1(var_list,varargin{:});
Error in ==> RBC_Est at 132
dynare_estimation(var_list_);
Error in ==> dynare at 102
evalin(‘base’,fname) ;
Could someone tell me the reason? (I did not do any change with the code; the softwares which I use are dynare 4.04 and Matlab R2008a ;the attachments are the code of RBC_Est.mod and the data)
Thank you in advance.
david
simuldataRBC.m (141 KB)
RBC_Est.mod (808 Bytes)