不知道怎么回事，程序都连在一起了，现在再发一遍程序，麻烦大家帮忙看看，谢谢各位老师
var c h y i q r dp n m k w;
varexo er;
parameters beta1 beta2 sigma d1 phi x nu miu sigma_k theta rho_r theta_dp theta_y theta_q phi_r j dpbar mbar rbar cbar ybar nbar kbar hbar qbar;
beta1=0.9;
beta2=0.95;
sigma=0.65;
d1=0.7;
phi=0.7;
x=0.3;
nu=0.3;
miu=0.3;
sigma_k=0.6;
theta=0.75;
rho_r=0.75;
theta_dp=1.5;
theta_y=0.6;
theta_q=0.35;
phi_r=0.7;
j=0.65;
dpbar=1;
rbar=1/beta2;
cbar=(1-beta1)^(1/sigma)/(1-j);
mbar=(cbar-j*cbar)^sigma/(1-beta1);
ybar=(1-beta2)*(cbar-j*cbar)^sigma/nu/x/(1-beta1+beta1*d1-beta2*d1);
nbar=((1-nu-miu)*x*y/(cbar-j*cbar)^sigma)^(1/(1-phi));
kbar=miu*x*ybar/(1-beta2*(1-sigma_k));
hbar=(ybar/kbar^miu/nbar^(1-nu-miu))^(1/nu);
qbar=nu*x*ybar/(1-beta2)/hbar;
model(linear);
w=phi*n+sigma*(c-j*c(-1));
-h/hbar=qbar/(cbar-j*cbar)^sigma*(q-sigma*(c-j*c(-1)))-beta1*qbar/(cbar-j*cbar)^sigma*(q(+1)-sigma*(c(+1)-j*c))+beta1*d1*qbar/(cbar-j*cbar)^sigma*(q(+1)-sigma*(c(+1)-j*c))-d1*qbar*dpbar/rbar/(cbar-j*cbar)^sigma*(q(+1)+dp(+1)-r-sigma*(c(+1)-j*c));
-m/mbar=(-sigma*(c-j*c(-1)))/(cbar-j*cbar)^sigma+beta1/(cbar-j*cbar)^sigma/dpbar*(-sigma*(c(+1)-j*c)-dp(+1));
-c=r-c(+1)-dp(+1);
miu*x*ybar/kbar*(y-c-k)=beta2*(1-sigma_k)*c(+1)-c;
nu*x*ybar/hbar*(y-c-h)=qbar*(q-c)-beta2*qbar*(q(+1)-c(+1));
(1-nu-miu)*x*(y-n)=w;
y=0.0968*c+0.1452*i;
r=rho_r*r(-1)+(1-rho_r)*(theta_dp*dp+theta_y*y+theta_q*q);
dp=beta1*dp(+1)+(1-beta1*theta)*(1-theta)/theta*y;
r=r(-1)*phi_r+er;
end;
initval;
c=0;
h=0;
n=0;
m=0;
q=0;
r=0;
i=0;
dp=0;
w=0;
k=0;
y=0;
er=0;
end;
steady;
check;
shocks;
var er;stderr 0.1;
end;
estimated_params;
phi_r,beta_pdf,0.7,0.1;
stderr er,inv_gamma_pdf,0.1,0.1;
end;
varobs dp;
estimation(datafile=cpi,prefilter=1,nobs=1,first_obs=197,mode_compute=6,mode_check,moments_varendo,order=1);
stoch_simul(periods=2000,irf=100);