全部版块 我的主页
论坛 经济学论坛 三区 宏观经济学
2727 5
2021-10-20
// 内生变量 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']) ;

二维码

扫码加我 拉你入群

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

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

全部回复
2021-10-21 10:15:24
你这个也太复杂了
二维码

扫码加我 拉你入群

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

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

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

扫码加我 拉你入群

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

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

2021-10-21 10:29:43
110199837 发表于 2021-10-21 10:15
你这个也太复杂了
不知道是哪里出了问题
二维码

扫码加我 拉你入群

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

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

2021-10-22 10:51:33
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']) ;
二维码

扫码加我 拉你入群

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

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

2022-11-9 22:57:13
wwwwwwwwxy 发表于 2021-10-22 10:51
确实是参数变量的问题,经过修改之后的代码如下:
但是依然出现了问题,麻烦帮忙看一下,万分感谢!
// ...
层主,怎么解决这个问题的,我现在也遇到了
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群