李老师，我主要是不太理解这个报错的原因在哪里。code如下
var lambda c tauc l taul w r q tauk rk v i y k g b z ua uz uk ul uc ug;
varexo e_a e_z e_k e_l e_c e_g;
parameters gamma psi kappa beta delta_0 delta_1 delta_2 tau alpha phiz phik phil gammaz gammak gammal gammag rhoz rhok rhol rhoa rhog rhoc;
parameters lambdas cs taucs ls tauls ws rs qs tauks rks vs is ys ks gs bs zs uas uzs uks uls ucs;
parameters sg sb ;
%参数赋值
gamma=1.75;
psi= 1;
kappa= 2;
beta=0.99;
delta0 =0.025;
delta1=0;
delta2=0.7;
tau=2;
alpha=0.33;
%phiz=0.1; phik=0.1; phil=0.1; gammaz=0.1; gammak=0.1; gammal=0.1; gammag=0.1; rhoz=0.9 ;rhok=0.9 ;rhol=0.9; rhoa=0.9; rhog=0.9;rhoc=0.9;
sg=0.1 ;
sb=0.3 ;
%稳态值
vs=1;
qs=1;
ks=(rks/alpha)^((1-alpha)/(kappa+alpha)*(kappa+gamma))*((1-tauls)*(1-alpha))^(1/(kappa+gamma))*((1-sg)*rks/alpha-delta0)^(-gamma/(kappa+alpha));
lambdas=cs^(-gamma);
ys=rks*ks/alpha;
cs=(1-gs/ys)*ys-ks*delta_0;
rs=1/beta;
rks=(1/beta-1+delta_0)/(1-tauks);
ls=(ys/ks^alpha)^(1/(1-alpha));
ws=(1-alpha)*ys/ls;
is=delta0*ks;
gs=sg*ys;
bs=sb*ys;
tauls=0.3 ;
tauks=0.2 ;
taucs=0;
zs=tauks*alpha*ys+tauls*(1-alpha)*ys-ys*(sg+(1-beta)*sb/beta);
uzs=1;
uks=1;
uls=1;
ucs=1;

model;
exp(lambda)=1/exp(c*gamma)*(tauc+1);
exp(lambda)=psi*exp(l*kappa)/(1-taul)*exp(w);
exp(lambda)/exp(lambda(+1))=beta*r;
q=lambda(+1)/lambda*beta*(1-tauk(+1)*rk(+1)*v(+1)+1-(delta_0+delta_1*(v(+1)-1)+1/2*(v(+1)-1)^2)*q(+1));
(1-tauk)*rk=q*(delta_1+delta_2*(v-1));
1=q*(1-tau*(exp(i)/exp(i(-1))-1)*(3/2*exp(i)/exp(i(-1))-1)+beta*q(+1)*lambda(+1)/lambda*tau*(exp(i(+1))/exp(i)-1)*(exp(i(+1))/exp(i))^2);
exp(w)=(1-alpha)*exp(y)/exp(l);
rk*v=alpha*exp(y)/exp(k(-1));
exp(y)=exp(ua)*v*exp(k*alpha)*exp((1-alpha)*l);
exp(k)=(1-(delta_0+delta_1*(v-1)+delta_2*1/2*(v-1)^2))*k(-1)+(1-tau/2*(exp(i)/exp(i(-1))-1)^2)*i;
exp(y)=exp(c)+exp(i)+exp(g);
exp(b)+tauk*alpha*exp(y)+taul*(1-alpha)*exp(y)+tauc*exp(c)=r(-1)*b(-1)+g+z;
%财政规则
exp(z)=-phiz*exp(y)-gammaz*exp(b(-1))/exp(y(-1))+exp(uz);
tauk=phik*exp(y)+gammak*exp(b(-1))/exp(y(-1))+exp(uk);
taul=phil*exp(y)+gammal*exp(b(-1))/exp(y(-1))+exp(ul);
tauc=exp(uc);
exp(g)=-gammag*exp(b(-1))/exp(y(-1))+exp(ug);
%冲击
ua=rhoa*ua(-1)+ e_a;
uz=rhoz*uz(-1)+ e_z;
uk=rhok*uk(-1)+ e_k;
ul=rhol*ul(-1)+ e_l;
uc=rhoc*uc(-1)+ e_c;
ug=rhog*ug(-1)+ e_g;

end;

initval;
lambda=log(lambdas);
c=log(cs) ;
tauc=taucs ;
l=log(ls);
taul=tauls;
w=log(ws) ;
r=rs;
q=qs ;
tauc=taucs;
tauk=tauks;
rk=rks ;
v=vs ;
i=log(is);
y=log(ys) ;
k=log(ks) ;
g=log(gs) ;
b=log(bs) ;
z=log(zs);
ua=log(uas);
uz=log(uzs);
uk=log(uks);
ul=log(uls);
uc=log(ucs);
end;

resid;
steady;

shock;
var e_a ;
stderr 1;
var e_z ;
stderr 1;
var e_k ;
stderr 1;
var e_l ;
stderr 1;
var e_c;
stderr 1;
var e_g ;
stderr 1;
end;
stoch_simul(order=1) y k i c g b ;

你好， gandz123456，非常高兴你能参与到DSGE课程的学习中来。很抱歉，最近比较忙，刚刚回复你。
对于你提出的这个错误信息，其实是你的语法错了，外生shock的声明模块，开头应该用shocks, 而不是shock。即添加s即可。

另外，dynare仍然报错，请你检查一下 下面这些参数是否都有 initial values (delta_0, delta_1, delta_2, phiz, phik, phil, gammaz, gammak, gammal, gammag, rhoz, rhok, rhol, rhoa, rhog, rhoc, lambdas, cs, ls,ws, rks, is, ys, ks, gs, bs, zs, uas)

if errors persist, please post your error info and code again. I will help you to check it.

我明白啦，谢谢，李老师。

求助李老师，我在做贝叶斯估计的时候遇到了问题，就是图只显示先验概率，而且dynare对于结果认为很有可能是错误的，而且貌似变成了极大似然估计了。我对于贝叶斯的代码不是太熟悉，希望老师帮我看看究竟是哪里的问题。万分感谢！！！
代码和数据均在附件内。