楼主: 孤星者
338 7

dynare运行问题There are 9 eigenvalue(s) larger than 1 in modulus for 8 forward- [推广有奖]

  • 0关注
  • 0粉丝

高中生

50%

还不是VIP/贵宾

-

威望
0
论坛币
5 个
通用积分
0.8033
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
139 点
帖子
7
精华
0
在线时间
52 小时
注册时间
2018-5-11
最后登录
2024-4-5

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
var Ni C wi Nj wj rK RD pi m K I di Yi Ai Ki mci piistar piii xi1 xi2 pii dj Yj Aj Kj mcj pijstar pijj xj1 xj2 pij Y gm;
varexo ei ej em;
parameters etai sigma etaj beta delta alpha theta varepsilon k phi rhom rhoi rhoj;
parameters dis djs pis piistars piiis piis pijstars pijjs pijs Ais Ajs rKs RDs mcis mcjs gms wjs wis Kis Kjs Nis Njs Cs Yis Yjs Ks Is Ys xi1s xi2s xj1s xj2s ms;

etai = 1;
sigma = 2;
etaj= 1.803228369;
beta = 0.99;
delta = 0.025;
alpha = 1/3;
theta = 0.75;
varepsilon = 10;
k = 0.8;
phi = 2;
rhom = 0.5;
rhoi = 0.9583;
rhoj = 0.7348;
dis =1;
djs =1;
pis =1;
piistars = 1;
piiis =1;
piis =1;
pijstars =1;
pijjs =1;
pijs =1;
Ais = 1;
Ajs = 0.5;
rKs = 1/beta-1+delta;
RDs = pis/beta;
mcis = (varepsilon-1)/varepsilon;
mcjs = (varepsilon-1)/varepsilon;
gms = 0;
wjs = Ajs*(Ajs*mcjs*alpha/rKs)^(alpha/(1-alpha))*mcjs*(1-alpha) ;
wis = Ais*(Ais*mcis*alpha/rKs)^(alpha/(1-alpha))*mcis*(1-alpha);
Kis = (wis*(Ais*alpha*mcis/rKs)^(1/(1-alpha))/(Ais/k/(Ais*alpha*mcis/rKs)-delta/k)^sigma)^(1/(sigma+1));
Kjs = Kis/k*(1-k);
Nis = Kis/(Ais*alpha*mcis/rKs)^(1/(1-alpha));
Njs = wis*Nis*(1-k)/wjs/k;
Cs = (wis/Nis)^(1/sigma);
Yis = Ais*Kis^alpha*Nis^(1-alpha);
Yjs = Ajs*Kjs^alpha*Njs^(1-alpha);
Ks = Kis+Kjs;
Is = delta*Ks;
Ys = Is+Cs;
xi1s = Cs^(-sigma)*mcis*Yis/(1-beta*theta);
xi2s = Cs^(-sigma)*Yis/(1-beta*theta);
xj1s = Cs^(-sigma)*mcjs*Yjs/(1-beta*theta);
xj2s = Cs^(-sigma)*Yjs/(1-beta*theta);
ms = Cs^sigma/(1-1/RDs);


model;
%家庭
%(1) 新兴产业劳动供给
Ni^(etai) = C^(-sigma)*wi;

%(2) 传统产业劳动供给
Nj^(etaj) = C^(-sigma)*wj;

%(3) 资本存量一阶条件
beta*(rK+1-delta) = (C/C(-1))^sigma;

%(4) 债券的一阶条件
beta*RD/pi(+1)*C(+1)^(-sigma)=C^(-sigma);

%(5) 货币需求方程
m^(-1)=C^(-sigma)*(1-1/RD);

%(6) 资本积累
K = (1-delta)*K(-1)+I;

% 中间品厂商i
%(7) 中间品厂商新兴产业生产函数
di*Yi = Ai*(Ki)^(alpha)*Ni^(1-alpha);

%(8) 边际成本方程
rK=mci*alpha*Yi/Ki;

%(9) 资本收益方程
wi=mci*(1-alpha)*Yi/Ni;

%(10) 价格离散核递归方程
di=(1-theta)*(piistar)^(-varepsilon)*piii^(varepsilon)+theta*(piii)^(varepsilon)*di(-1);

%(11) CPI通胀递归方程
piii^(1-varepsilon)=(1-theta)*(piistar)^(1-varepsilon)+theta;

%(12) 最优价格方程
piistar=varepsilon/(varepsilon-1)*piii*xi1/xi2;

%(13) 最优价格辅助方程
xi1=C^(-sigma)*mci*Yi+beta*theta*xi1(+1)*(piii(+1))^(varepsilon);

%(14) 最优价格辅助方程
xi2=C^(-sigma)*pii/pi*Yi+beta*theta*xi2(+1)*(piii(+1))^(varepsilon-1);

% 中间品厂商j
%(15) 中间品厂商传统产业生产函数
dj*Yj = Aj*(Kj)^(alpha)*Nj^(1-alpha);

%(16) 边际成本方程
rK=mcj*alpha*Yj/Kj;

%(17) 资本收益方程
wj=mcj*(1-alpha)*Yj/Nj;

%(18) 价格离散核递归方程
dj=(1-theta)*(pijstar)^(-varepsilon)*pijj^(varepsilon)+theta*(pijj)^(varepsilon)*dj(-1);

%(19) CPI通胀递归方程
pijj^(1-varepsilon)=(1-theta)*(pijstar)^(1-varepsilon)+theta;

%(20) 最优价格方程
pijstar=varepsilon/(varepsilon-1)*pijj*xj1/xj2;

%(21) 最优价格辅助方程
xj1=C^(-sigma)*mcj*Yj+beta*theta*xj1(+1)*(pijj(+1))^(varepsilon);

%(22) 最优价格辅助方程
xj2=C^(-sigma)*pij/pi*Yj+beta*theta*xj2(+1)*(pijj(+1))^(varepsilon-1);

%最终品厂商部门
%(23) 生产函数
Y = (k^(1/phi)*Yi^((phi-1)/phi)+(1-k)^(1/phi)*Yj^((phi-1)/phi))^(phi/(phi-1));

%(24) 新兴行业生产量
Yi = k*Y*(pi/pii)^phi;

%(25) 传统行业生产量
Yj = (1-k)*Y*(pi/pij)^phi;

%央行部门
%(26) 传统行业生产量
gm = -log(pi)+rhom*gm(-1)+rhom*log(pi(-1))+em;

%(27) 货币供给
gm = log(m)-log(m(-1));

%市场出清
%(28)
Y = C+I;

%(29)
K(-1) = Ki+Kj;

%(30)
pii/pii(-1) = piii/pi(-1);

%(31)
pij/pij(-1) = pijj/pi(-1);

%(32)
log(Ai)-log(Ais) = rhoi*(log(Ai(-1))-log(Ais))+ei;

%(33)
log(Aj)-log(Ajs) = rhoj*(log(Aj(-1))-log(Ajs))+ej;

end;

initval;               
Ni = Nis;
C = Cs;
wi = wis;
Nj = Njs;
wj = wjs;
rK = rKs;
RD = RDs;
pi = pis;
m = ms;
K = Ks;
I = Is;
di =dis;
Yi = Yis;
Ai = Ais;
Ki = Kis;
mci = mcis;
piistar =piistars;
piii = piiis;
xi1 = xi1s;
xi2 = xi2s;
pii =piis;
dj =djs;
Yj = Yjs;
Aj = Ajs;
Kj = Kjs;
mcj = mcjs;
pijstar = pijstars;
pijj = pijjs;
xj1 = xj1s;
xj2 = xj2s;
pij = pijs;
Y = Ys;
gm = gms;
end;

shocks;
var ei = 0.01^2;
var ej = 0.01^2;
var em = 0.01^2;
end;

resid(1);
steady;
check;
stoch_simul(order=1);


二维码

扫码加我 拉你入群

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

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

关键词:eigenvalue forward modulus dynare value DSGE;Dynare; dynare求助

沙发
孤星者 学生认证  发表于 2023-11-19 18:15:02 |只看作者 |坛友微信交流群
Starting Dynare (version 5.4).
Calling Dynare with arguments: none
Starting preprocessing of the model file ...
Found 33 equation(s).
Evaluating expressions...done
Computing static model derivatives (order 1).
Computing dynamic model derivatives (order 1).
Processing outputs ...
done
Preprocessing completed.





Residuals of the static equations:

Equation number 1 : 0 : 1
Equation number 2 : 0 : 2
Equation number 3 : 0 : 3
Equation number 4 : 0 : 4
Equation number 5 : 0 : 5
Equation number 6 : 0 : K
Equation number 7 : 0 : 7
Equation number 8 : 0 : rK
Equation number 9 : 0 : wi
Equation number 10 : 0 : di
Equation number 11 : 0 : 11
Equation number 12 : 0 : piistar
Equation number 13 : 0 : xi1
Equation number 14 : 0 : xi2
Equation number 15 : 0 : 15
Equation number 16 : 0 : 16
Equation number 17 : 0 : wj
Equation number 18 : 0 : dj
Equation number 19 : 0 : 19
Equation number 20 : 0 : pijstar
Equation number 21 : 0 : xj1
Equation number 22 : 0 : xj2
Equation number 23 : 0 : Y
Equation number 24 : 0 : Yi
Equation number 25 : 0 : Yj
Equation number 26 : 0 : gm
Equation number 27 : 0 : 27
Equation number 28 : 0 : 28
Equation number 29 : 0 : 29
Equation number 30 : 0 : 30
Equation number 31 : 0 : 31
Equation number 32 : 0 : 32
Equation number 33 : 0 : 33



STEADY-STATE RESULTS:

Ni                       0.596678
C                        1.71457
wi                       1.75409
Nj                       0.421915
wj                       0.620165
rK                       0.035101
RD                       1.0101
pi                       1
m                        293.976
K                        18.636
I                        0.4659
di                       1
Yi                       1.74438
Ai                       1
Ki                       14.9088
mci                      0.9
piistar                  1
piii                     1
xi1                      2.07393
xi2                      2.30437
pii                      1
dj                       1
Yj                       0.436095
Aj                       0.5
Kj                       3.7272
mcj                      0.9
pijstar                  1
pijj                     1
xj1                      0.518482
xj2                      0.576091
pij                      1
Y                        2.18047
gm                       0

EIGENVALUES:
         Modulus             Real        Imaginary

       1.296e-14       -1.296e-14                0
       2.321e-09        2.321e-09                0
       2.321e-09       -2.321e-09                0
             0.5              0.5                0
          0.6413           0.6413                0
          0.7348           0.7348                0
            0.75             0.75                0
            0.75             0.75                0
          0.9583           0.9583                0
          0.9616           0.9616                0
            1.01             1.01                0
           1.057            1.057                0
           1.347            1.347                0
           1.347            1.347                0
           1.574            1.574                0
           2.788            2.788                0
       1.192e+32        1.192e+32                0
       1.535e+32        1.535e+32                0
       8.359e+32       -8.359e+32                0


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

The rank condition ISN'T verified!

错误使用 print_info
Blanchard & Kahn conditions are not satisfied: no stable
equilibrium.
出错 stoch_simul (第 120 行)
    print_info(info, options_.noprint, options_);
出错 first.driver (第 597 行)
[info, oo_, options_, M_] = stoch_simul(M_, options_, oo_, var_list_);
出错 dynare (第 278 行)
    evalin('base',[fname '.driver']);

使用道具

藤椅
孤星者 学生认证  发表于 2023-11-19 18:15:30 |只看作者 |坛友微信交流群
出现如下问题,请问有大佬知晓问题吗?

使用道具

板凳
孤星者 学生认证  发表于 2023-11-19 18:15:33 |只看作者 |坛友微信交流群
出现如下问题,请问有大佬知晓问题吗?

使用道具

报纸
oliyiyi 发表于 2023-12-20 09:21:03 |只看作者 |坛友微信交流群
模型中的前瞻变量(forward-looking variables)太多,超过了稳定的特征根(stable eigenvalues)的个数。这就导致了阶跃响应(impulse response)不确定。

使用道具

地板
WELCOMETOHELP 学生认证  发表于 2023-12-29 14:27:10 |只看作者 |坛友微信交流群
孤星者 发表于 2023-11-19 18:15
Starting Dynare (version 5.4).
Calling Dynare with arguments: none
Starting preprocessing of the  ...
楼主,你好,请问你是如何解决的?

使用道具

7
WELCOMETOHELP 学生认证  发表于 2023-12-29 14:27:55 |只看作者 |坛友微信交流群
oliyiyi 发表于 2023-12-20 09:21
模型中的前瞻变量(forward-looking variables)太多,超过了稳定的特征根(stable eigenvalues)的个数。这就导 ...
请问该如何解决,删减变量吗?

使用道具

8
ggfree 学生认证  发表于 2024-3-10 15:33:34 |只看作者 |坛友微信交流群
楼主,您好,请问解决了吗?

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加JingGuanBbs
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-27 14:21