STIMATION_CHECKS: There was an error in computing the likelihood for initial parameter values.
If this is not a problem with the setting of options (check the error message below),you should try using the calibrated version of the model as starting values. To do
this, add an empty estimated_params_init-block with use_calibration option immediately before the estimation
ESTIMATION_CHECKS: command (and after the estimated_params-block so that it does not get overwritten):


我在估计的时候老是呈现这个问题，不知道怎么解决？请老师指点下。
var y i k c r n m pi w r_k q z b a mu_p mu_w theta ;//  endogenous17

varexo e_z e_b e_a e_theta e_mu_w  e_mu_p ; //外生冲击6

parameters  beta sigma gamma phi delta alpha psi psi_u mu
            epsilon_p epsilon_w rho_z  rho_a  rho_b rho_mu_p
            rho_mu_w rho_theta sigma_z sigma_b sigma_a sigma_mu_p
            sigma_mu_w  sigma_theta; //23
beta =0.9701;
sigma =0.36;
gamma =0.6;
phi =0.42;
delta =0.025;
alpha =0.41;
psi =0.15;
psi_u = 0.17;
mu =1.2;

epsilon_p=0.75;
epsilon_w =0.75;
rho_z  =0.85;
rho_a =0.85;
rho_b =0.85;
rho_mu_p=0.85;            
rho_mu_w =0.85;
rho_theta =0.85;

sigma_z =0.04;
sigma_b =0.1;
sigma_a =0.04;
sigma_mu_p=0.05;            
sigma_mu_w  =0.25;
sigma_theta =0.014;

model(linear);
   -(1/gamma)*m + c/sigma  = beta/(1-beta)*r;  //1
    n = k(-1)-w +(1+psi_u)*r_k;   //2
    c =c(+1) +sigma*(-r+pi(+1)+(1-rho_b)*b);//3
    i = i(-1)/(1+beta)+beta*i(+1)/(1+beta)+psi*q/(1+beta)-(1-beta*rho_z)*z/(1+beta);//4
    q = beta *(1-delta)*q(+1) - r +pi(+1)+(1-beta*(1-delta))*r_k(+1);//5
    w = (1/(1+beta))*(pi(-1)+w(-1)+beta*(w(+1)+pi(+1))-(1-beta*epsilon_w)*(1-epsilon_w)*(w-mu_w-n/phi-c/sigma)/epsilon_w)-pi;//6
   pi = (1/(1+beta))*(pi(-1)+beta*pi(+1)+(1-beta*epsilon_p)*(1-epsilon_p)*(alpha*r_k +(1-alpha)*w -a+mu_p)/epsilon_p);//7
    y = c*(1-alpha*beta*delta)/(mu*(1-beta*(1-delta)))+alpha*beta*delta*i/(mu*(1-beta*(1-delta)));//8
    y = a +alpha*k(-1)+alpha*psi_u*r_k+(1-alpha)*n;//9
    k =(1-delta)*k(-1)+delta*i;//10
theta = m -m(-1)+pi;  //11
    z = rho_z*z(-1)+sigma_z*e_z;//12
    b = rho_b*b(-1)+sigma_b*e_b;  //13
    a = rho_a*a(-1)+sigma_a*e_a;//14
mu_p = rho_mu_p*mu_p(-1)+sigma_mu_p*e_mu_p;  //15
theta = rho_theta*theta(-1)+sigma_theta*e_theta;  //16
mu_w = rho_mu_w*mu_w(-1)+sigma_mu_w*e_mu_w;   //17
end;
initval;
y = 0;
i=0;
k =0;
c =0;
r =0;
n =0;
m =0;
pi =0;
w =0;
r_k =0;
q =0;
z =0;
b =0;
a =0;
mu_p =0;
mu_w =0;
theta=0 ;
end;

shocks;
var e_z   ; stderr 1;
var e_a   ; stderr 1;
var e_b   ; stderr 1;
var e_theta ; stderr 1;
var e_mu_p ; stderr 1;
var e_mu_w ; stderr 1;
end;
steady;
check;

estimated_params;
sigma_a, inv_gamma_pdf,0.4,30;
sigma_b, inv_gamma_pdf,0.2,30;
sigma_z,    inv_gamma_pdf,0.1,30;
sigma_theta, inv_gamma_pdf,0.1,30;
sigma_mu_p, inv_gamma_pdf,0.15,30;
sigma_mu_w, inv_gamma_pdf,0.25,30;
psi,normal_pdf,0.15,0.075;
psi_u,normal_pdf,0.17,0.075;
sigma,normal_pdf,0.36,0.03;
gamma,normal_pdf,0.6,0.2;
rho_a,beta_pdf,0.85,0.1;
rho_z,beta_pdf,0.85,0.1;
rho_b,beta_pdf,0.85,0.1;
rho_theta,beta_pdf,0.85,0.1;
rho_mu_p,beta_pdf,0.85,0.1;
rho_mu_w,beta_pdf,0.85,0.1;


epsilon_p,beta_pdf,0.75,0.15;
epsilon_w,beta_pdf,0.75,0.15;

end;


varobs  y c pi w r m;
estimation(datafile=chap3data,mh_replic=2000,mh_nblocks=4,mh_jscale=1);
stoch_simul y c;