求助，dynare一直报错。-经管之家官网！

我在运行dynare的时候一直出现error: dynare (line 229)DYNARE: preprocessing failed的问题，求大佬看下是怎么回事。下面是代码addpath C:\dynare\4.5.6\matlab
dynare reported2.mod
var c i y w l u q r rk k pi ypot e_a e_b e_g e_i e_l e_u piba dc dh di dp dw dy r=obs;
varexo v_a v_b v_g v_i v_l v_u v_pi np nw nq nr;
parameters alpha beta delta cyss gyss lambdaw lambdap theta sigma chi
zetainv psi phi_1 xip xiw gammap gammaw r_i r_pi r_y r_dpi r_dy rho_a rho_b
rho_g rho_i rho_l rho_u rho_pi zst lst pist rst;
alpha=0.37;
beta=0.995;
delta=0.06/4;
cyss=0.6;
gyss=0.2;
lambdap=0.2;
lambdaw=0.2;
model(linear);
rkba = 1/beta-(1-delta);
np = 0;
nq = 0;
nw = 0;
c = 1/(1+theta+beta*(theta^2))*(theta*c(-1)+(1+beta*(theta^2)+beta*theta)*c(+1)
-beta*theta*c(+2)-((1-theta)/sigma)*((1-beta*theta)*(r-pi(+1))-e_b+
(1+beta*theta)*e_b(+1)-beta*theta*e_b(+2)));
q = -(r-pi(+1))+(1/(1-delta+rkba))*((1-delta)*q(+1)+rkba*rk(+1))+nq;
i = (1/(1+beta))*(i(-1)+beta*i(+1)+(1/zetainv)*q+beta*e_i(+1)-e_i);
k = (1-delta)*k(-1)+delta*i-rkba*(u+e_u);
u = psi*(rk-q-e_u)-e_u;
l = -w+rk+u+k(-1);
y = cyss*c+gyss*e_g+delta*i*k;
y = (1+phi_1)*(e_a+alpha*(u+k(-1))+(1-alpha)*l);
w = (1/(1+beta))*(beta*w(+1)+w(-1)+beta*pi(+1)-(1+beta*gammaw)*pi+gammaw*pi(-1)
-(lambdaw*(1-beta*xiw)*(1-xiw)/((lambdaw+(1+lambdaw)*chi)*xiw))*
(w-chi*l-e_l+(beta*theta/(1-beta*theta))*(e_b-e_b(+1))-nw-
(sigma/((1-theta)*(1-beta*theta)))*((1+beta*(theta^2))*c-theta*c(-1)-
beta*theta*c(+1))));
pi = (1/(1+beta*gammap))*(beta*pi(+1)+gammap*pi(-1)+((1-beta*xip)*(1-xip)/xip)
*(w+alpha*(l-u-k(-1))-e_a+np));
r = r_i*r(-1)+(1-r_i)*piba+r_pi*(pi(-1)-piba)+r_y*(y(-1)-ypot(-1))+r_dpi*(
pi-pi(-1))+r_dy*((y-ypot)-(y(-1)-ypot(-1)))+nr;
ypot = -alpha*(1+phi_1)*e_a;
e_a = rho_a*e_a(-1)+v_a;
e_b = rho_b*e_b(-1)+v_b;
e_g = rho_g*e_g(-1)+v_g;
e_i = rho_i*e_i(-1)+v_i;
e_l = rho_l*e_l(-1)+v_l;
e_u = rho_u*e_u(-1)+v_u;
piba = rho_pi*piba(-1)+v_pi;
dc = zst+e_a+c-c(-1);
dh = lst+1;
di = zst+e_a+i-i(-1);
dp = pist+pi;
dw = zst+e_a+w-w(-1);
dy = zst+e_a+y-y(-1);
r_obs = rst+pist+r;
end;
estimated_params;
theta,beta_pdf,0.7,0.15;
sigma,normal_pdf,1,0.375;
chi,normal_pdf,2,0.75;
zetainv,normal_pdf,4,1.5;
mu,normal_pdf,1,1;
phi_1,gamma_pdf,0.075,0.0125;
xip,beta_pdf,0.375,0.1;
xiw,beta_pdf,0.375,0.1;
gammap,beta_pdf,0.5,0.25;
gammaw,beta_pdf,0.5,0.25;
r_i,normal_pdf,1,0.15;
r_pi,normal_pdf,0.5,0.2;
r_y,normal_pdf,0.01,0.1;
r_dpi,normal_pdf,0.1,0.1;
r_dy,normal_pdf,0.1,0.5;
stderr v_a,inv_gamma_pdf,1,inf;
stderr v_pi,inv_gamma_pdf,0.1,inf;
stderr v_b,inv_gamma_pdf,0.1,inf;
stderr v_g,inv_gamma_pdf,0.4,inf;
stderr v_l,inv_gamma_pdf,2,inf;
stderr v_i,inv_gamma_pdf,1,inf;
stderr v_u,inv_gamma_pdf,0.5,inf;
stderr nr,inv_gamma_pdf,0.1,inf;
stderr nq,inv_gamma_pdf,,,;
stderr np,inv_gamma_pdf,0.2,inf;
stderr nw,inv_gamma_pdf,0.2,inf;
rho_a,beta_pdf,0.85,0.1;
rho_pi,beta_pdf,0.85,0.1;
rho_b,beta_pdf,0.85,0.1;
rho_g,beta_pdf,0.85,0.1;
rho_l,beta_pdf,0.85,0.1;
rho_i,beta_pdf,0.85,0.1;
rho_u,beta_pdf,0.85,0.1;
varobs dc dh di dp dw dy r_obs;
estimation(datafile = DataSet,mode_check,mh_replic = 500000,mh_nblocks=2,
mh_jscale=0.35,compute_mode=4,bayesian_irf,smoother);
stoch_simul(irf = 40, order = 1) y c i w l u pi r;