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 R2008 ;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)