求助！！dynare报错！！！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

// 内生变量 31
var c k d i ni nj wi wj rd rk ki kj gs yi a kg gi ri li yj rj lj f l t e g y n lambda mu;
// 外生变量
varexo eps_a eps_ig eps_f eps_e ;
// 参数
parameters beta thet1 thet2 alpha1 alpha2 delta_k delta_g tau s rho_a rho_ig rho_f rho_e
A LI_L LJ_L RK RD MU RI RJ F YI_Y YJ_Y GS_G GI_G C_Y LI_YI LJ_YJ LI_LJ YI_YJ YJ_YI KI_KJ KI_K KJ_K KI_YI KJ_YJ WINI_YI WJNJ_YJ K_Y I_Y G_Y T_Y LI_Y E_Y NI_NJ NI_N NJ_N;
beta=0.99;thet1=0.3;thet2=0.2;alpha1=0.45;alpha2=0.2;delta_k=0.04;delta_g=0.04;tau=0.09;s=0.15;
rho_a=0.75;rho_ig=0.35;rho_f=0.9;rho_e=0.9;
A=1;
LI_L=0.09;
LJ_L=1-LI_L;
RK=1/beta+delta_k-1;
RD=1/beta;
MU=(s-RD)/(s-1);
RI=MU/1.5;
RJ=MU;
F=0.5*RI;
GS_G=0.2;
GI_G=0.8;
C_Y=0.384;

LI_YI=(1-tau)/RI;
LJ_YJ=(1-tau)/RJ;
LI_LJ=LI_L/(1-LI_L);
YI_YJ=LI_LJ*(RI/RJ);
YJ_YI=1/YI_YJ;
YI_Y=YI_YJ/(YI_YJ+1);
YJ_Y=1-YI_Y;

KI_KJ=YI_YJ*(RJ/RI);
KI_K=KI_KJ/(KI_KJ+1);
KJ_K=1-KI_K;
KI_YI=(1-tau)*alpha1/(RI*(RK+delta_k));
KJ_YJ=(1-tau)*alpha1/(RJ*(RK+delta_k));

WINI_YI=(1-tau)*(1-alpha)/RI;
WJNJ_YJ=(1-tau)*(1-alpha)/RJ;
K_Y=KJ_YJ*((1+KI_KJ)/(1+YI_YJ));
I_Y=delta_k*K_Y;
G_Y=1-C_Y-I_Y;
T_Y=tau;
LI_Y=LI_YI/(1+YJ_YI);
E_Y=G_Y+F*LI_Y-T_Y;

NI_NJ=((YI_YJ)*(KI_KJ^alpha1))^(1/(1-alpha1));
NI_N=NI_NJ/(NI_NJ+1);
NJ_N=1-NI_N;

// 26+5 Equations of the model
predetermined_variables k ;
model(linear);
//家庭
k(+1)=delta_k*i+(1-delta_k)*k;
lambda+c=0;
wi=c+thet1*ni;
wj=c+thet2*nj;
c(+1)=c+beta*(1/beta-1+delta_k)*rk(+1);
c(+1)=c+rd;
//企业
kg(+1)=delta_g*gi+(1-delta_g)*kg;
yi=a+alpha1*ki+alpha2*kg+(1-alpha1)*ni;
LI_YI*li=WINI_YI*(wi+ni)+RK*KI_YI*(rk+ki)+delta_k*KI_YI*ki;
yi=ki+ni+wi;
yi+ki+ri+rk*RK/(RK+delta_k);
yj=a+alpha1*kj+alpha2*kg+(1-alpha1)*nj;
LJ_YJ*lj=WJNJ_YJ*(wj+nj)+RK*KJ_YJ*(rk+kj)+delta_k*KJ_YJ*kj;
yj=kj+nj+wj;
yj+kj+rj+rk*RK/(RK+delta_k);
y=YI_Y*yi+YJ_Y*yj;
l=LI_L*li+LJ_L*lj;
//银行
l=d;
MU*mu=RI*ri+F*f;
mu=rj;
mu=rd;
//政府
T_Y*t+E_Y*e=G_Y*g+F*LI_Y*(f+li);
t=y;
g=GS_G*gs+GI_G*gi;
//市场出清
y=C_Y*c+I_Y*i+G_Y*g;
k=KI_K*ki+KJ_K*kj;
n=NI_N*ni+NJ_N*nj;
//冲击
a=rho_a*a(-1)+eps_a;
gi=rho_ig*gi(-1)+eps_ig;
f=rho_f*f(-1)+eps_f;
e=rho_e*e(-1)+eps_e;

end;

steady;
check;

shocks;
var eps_a;  stderr 0.01;
var eps_ig;  stderr 0.01;
var eps_f;  stderr 0.01;
var eps_e;  stderr 0.01;
end;


stoch_simul(irf = 40) c wi wj n ni nj y yi yj i;




Using 64-bit preprocessor
Starting Dynare (version 4.6.4).
Calling Dynare with arguments: none
Starting preprocessing of the model file ...
Found 31 equation(s).
Evaluating expressions...done
Computing static model derivatives (order 1).
Computing dynamic model derivatives (order 2).
Processing outputs ...
done
Preprocessing completed.

错误使用 alpha
输出参数太多。

出错 Untitled3.driver (line 524)
M_.params(38) = (1-M_.params(8))*(1-alpha)/M_.params(20);

出错 dynare (line 293)
evalin('base',[fname '.driver']) ;

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：dynare ARE derivatives expressions Processing

你这个也太复杂了

对的参数名的问题，dynare无法区分你的alpha1写到yj=a+alpha1*kj+alpha2*kg+(1-alpha1)*nj;这里头比如是啥意思，还以为是数字1，无法识别

110199837 发表于 2021-10-21 10:15
你这个也太复杂了

不知道是哪里出了问题

110199837 发表于 2021-10-21 10:17
对的参数名的问题，dynare无法区分你的alpha1写到yj=a+alpha1*kj+alpha2*kg+(1-alpha1)*nj;这里头比如是啥 ...

确实是参数变量的问题，经过修改之后的代码如下：
但是依然出现了问题，麻烦帮忙看一下，万分感谢！
// 内生变量 31
var c k d i ni nj wi wj rd rk ki kj gs yi a kg gi ri li yj rj lj f l t e g y n lambda mu;
// 外生变量
varexo eps_a eps_ig eps_f eps_e ;
// 参数
parameters beta thet1 thet2 alpha1 alpha2 delta_k delta_g tau s rho_a rho_ig rho_f rho_e
A LI_L LJ_L RK RD MU RI RJ F YI_Y YJ_Y GS_G GI_G C_Y LI_YI LJ_YJ LI_LJ YI_YJ YJ_YI KI_KJ KI_K KJ_K KI_YI KJ_YJ WINI_YI WJNJ_YJ K_Y I_Y G_Y T_Y LI_Y E_Y NI_NJ NI_N NJ_N;
beta=0.99;thet1=0.3;thet2=0.2;alpha1=0.45;alpha2=0.2;delta_k=0.04;delta_g=0.04;tau=0.09;s=0.15;
rho_a=0.75;rho_ig=0.35;rho_f=0.9;rho_e=0.9;
A=1;
LI_L=0.09;
LJ_L=1-LI_L;
RK=1/beta+delta_k-1;
RD=1/beta;
MU=(s-RD)/(s-1);
RI=MU/1.5;
RJ=MU;
F=0.5*RI;
GS_G=0.2;
GI_G=0.8;
C_Y=0.384;

LI_YI=(1-tau)/RI;
LJ_YJ=(1-tau)/RJ;
LI_LJ=LI_L/(1-LI_L);
YI_YJ=LI_LJ*(RI/RJ);
YJ_YI=1/YI_YJ;
YI_Y=YI_YJ/(YI_YJ+1);
YJ_Y=1-YI_Y;

KI_KJ=YI_YJ*(RJ/RI);
KI_K=KI_KJ/(KI_KJ+1);
KJ_K=1-KI_K;
KI_YI=(1-tau)*alpha1/(RI*(RK+delta_k));
KJ_YJ=(1-tau)*alpha1/(RJ*(RK+delta_k));

WINI_YI=(1-tau)*(1-alpha1)/RI;
WJNJ_YJ=(1-tau)*(1-alpha1)/RJ;
K_Y=KJ_YJ*((1+KI_KJ)/(1+YI_YJ));
I_Y=delta_k*K_Y;
G_Y=1-C_Y-I_Y;
T_Y=tau;
LI_Y=LI_YI/(1+YJ_YI);
E_Y=G_Y+F*LI_Y-T_Y;

NI_NJ=((YI_YJ)*(KI_KJ^alpha1))^(1/(1-alpha1));
NI_N=NI_NJ/(NI_NJ+1);
NJ_N=1-NI_N;

// 27+4 Equations of the model
predetermined_variables k ;
model(linear);
//家庭
k(+1)=delta_k*i+(1-delta_k)*k;
lambda+c=0;
wi=c+thet1*ni;
wj=c+thet2*nj;
c(+1)=c+beta*(1/beta-1+delta_k)*rk(+1);
c(+1)=c+rd;
//企业
kg(+1)=delta_g*gi+(1-delta_g)*kg;
yi=a+alpha1*ki+alpha2*kg+(1-alpha1)*ni;
LI_YI*li=WINI_YI*(wi+ni)+RK*KI_YI*(rk+ki)+delta_k*KI_YI*ki;
yi=ki+ni+wi;
yi+ki+ri+rk*RK/(RK+delta_k);
yj=a+alpha1*kj+alpha2*kg+(1-alpha1)*nj;
LJ_YJ*lj=WJNJ_YJ*(wj+nj)+RK*KJ_YJ*(rk+kj)+delta_k*KJ_YJ*kj;
yj=kj+nj+wj;
yj+kj+rj+rk*RK/(RK+delta_k);
y=YI_Y*yi+YJ_Y*yj;
l=LI_L*li+LJ_L*lj;
//银行
l=d;
MU*mu=RI*ri+F*f;
mu=rj;
mu=rd;
//政府
T_Y*t+E_Y*e=G_Y*g+F*LI_Y*(f+li);
t=y;
g=GS_G*gs+GI_G*gi;
//市场出清
y=C_Y*c+I_Y*i+G_Y*g;
k=KI_K*ki+KJ_K*kj;
n=NI_N*ni+NJ_N*nj;
//冲击
a=rho_a*a(-1)+eps_a;
gi=rho_ig*gi(-1)+eps_ig;
f=rho_f*f(-1)+eps_f;
e=rho_e*e(-1)+eps_e;

end;

steady;
check;

shocks;
var eps_a;  stderr 0.01;
var eps_ig;  stderr 0.01;
var eps_f;  stderr 0.01;
var eps_e;  stderr 0.01;
end;
stoch_simul(irf = 40) c wi wj n ni nj y yi yj i;

错误：
STEADY-STATE RESULTS:

c                       0
k                       0
d                       0
i                       0
ni                      0
nj                      0
wi                      0
wj                      0
rd                      0
rk                      0
ki                      0
kj                      0
gs                      0
yi                      0
a                       0
kg                      0
gi                      0
ri                      0
li                      0
yj                      0
rj                      0
lj                      0
f                       0
l                       0
t                       0
e                       0
g                       0
y                       0
n                       0
lambda                  0
mu                      0

EIGENVALUES:
         Modulus             Real        Imaginary

            0.35             0.35                0
            0.75             0.75                0
          0.8648           0.8648                0
             0.9              0.9                0
             0.9              0.9                0
            0.96             0.96                0
           1.182            1.182                0
           10.04           -10.04                0


There are 2 eigenvalue(s) larger than 1 in modulus
for 3 forward-looking variable(s)

The rank condition ISN'T verified!

错误使用 print_info (line 32)
Blanchard & Kahn conditions are not satisfied: indeterminacy.

出错 stoch_simul (line 107)
    print_info(info, options_.noprint, options_);

出错 Untitled.driver (line 558)
[info, oo_, options_, M_] = stoch_simul(M_, options_, oo_, var_list_);

出错 dynare (line 293)
evalin('base',[fname '.driver']) ;

wwwwwwwwxy 发表于 2021-10-22 10:51
确实是参数变量的问题，经过修改之后的代码如下：
但是依然出现了问题，麻烦帮忙看一下，万分感谢！
// ...

层主，怎么解决这个问题的，我现在也遇到了