dynare运行出现Error: File: WWDD.m Line: 210 Column: 1

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：column dynare Error colum file

mod文件源码放上来看一下

きずな 发表于 2019-4-18 22:15
mod文件源码放上来看一下

var a n c k i r g w y m;
varexo ea ep;
parameters theta enta beta delta alpha1 alpha2 alpha3 sigma rho1 rho2 sda;
parameters ns cs ks as is rs ws gs ps;

alpha1 = 0.4;
alpha2 = 0.45;
alpha3 = 0.15;
beta  = 0.95;
delta = 0.11;
sigma = 1;
rho1 = 0.85;
rho2 = 0.75; %calibrated
theta = 0.84;
enta = 3;
sda = 0.009 %calibrated

as = 1;
ns = 1/3;
rs = 1/beta;
m = ns^(1+enta)/(1/beta+delta);
cs = (m*alpha1/alpha2)^(alpha1/alpha2)*(ns^(1+enta)*alpha3/alpha2)^(alpha3/alpha2)*ns;
ks = cs*m*alpha1/alpha2;
ys = ks*(1/beta+delta)/alpha1;
is=  delta*ks;
ws = alpha2*ys/ns;
gs = alpha3*ys/ps;


exp(n)^(enta)=alpha2*exp(y)/exp(n);

exp(p)=alpha3*exp(y)/exp(g);

1=beta*(exp(c(+1))^(-theta)/exp(c)^(-theta)*alpha1*exp(y(+1)/exp(k)-delta);


exp(r)=alpha1*exp(a)*exp(k(-1))^(alpha1)*exp(n)^(alpha2)*(exp(p)*exp(g))^(alpha3);


exp(w)=alpha2*exp(a)*exp(k(-1))^(alpha1)*exp(n)^(alpha2-1)*(exp(p)*exp(g))^(alpha3);


exp(y)=exp(c)+exp(i)+exp(p)*exp(g);
exp(y)=exp(a)*exp(k(-1))^(alpha1)*exp(n)^(alpha2)*(exp(p)*exp(g))^(alpha3);


exp(k)=exp(i)+(1-delta)*exp(k(-1));


a=rho1*a(-1)+ea;
p=rho2*p(-1)+ep;



end;
initval;
k=log(ks);
y=log(ys);
c=log(cs);
i=log(is);
a=log(as);
r=log(rs);
w=log(ws);
n=log(ns);
g=log(gs);
p=log(ps);
end;

shocks;
var ea = sda^2;
var ep = sda^2;
end;
resid(1);
steady;
check;

stoch_simul(order =1, hp_filter =1600,periods=1000);

きずな 发表于 2019-4-18 22:15
mod文件源码放上来看一下

var a n c k i r g w y m;
varexo ea ep;
parameters theta enta beta delta alpha1 alpha2 alpha3 sigma rho1 rho2 sda;
parameters ns cs ks as is rs ws gs ps;

alpha1 = 0.4;
alpha2 = 0.45;
alpha3 = 0.15;
beta  = 0.95;
delta = 0.11;
sigma = 1;
rho1 = 0.85;
rho2 = 0.75; %calibrated
theta = 0.84;
enta = 3;
sda = 0.009 %calibrated

as = 1;
ns = 1/3;
rs = 1/beta;
m = ns^(1+enta)/(1/beta+delta);
cs = (m*alpha1/alpha2)^(alpha1/alpha2)*(ns^(1+enta)*alpha3/alpha2)^(alpha3/alpha2)*ns;
ks = cs*m*alpha1/alpha2;
ys = ks*(1/beta+delta)/alpha1;
is=  delta*ks;
ws = alpha2*ys/ns;
gs = alpha3*ys/ps;


exp(n)^(enta)=alpha2*exp(y)/exp(n);


exp(p)=alpha3*exp(y)/exp(g);


1=beta*(exp(c(+1))^(-theta)/exp(c)^(-theta)*alpha1*exp(y(+1)/exp(k)-delta);


exp(r)=alpha1*exp(a)*exp(k(-1))^(alpha1)*exp(n)^(alpha2)*(exp(p)*exp(g))^(alpha3);




exp(w)=alpha2*exp(a)*exp(k(-1))^(alpha1)*exp(n)^(alpha2-1)*(exp(p)*exp(g))^(alpha3);


exp(y)=exp(c)+exp(i)+exp(p)*exp(g);
exp(y)=exp(a)*exp(k(-1))^(alpha1)*exp(n)^(alpha2)*(exp(p)*exp(g))^(alpha3);


exp(k)=exp(i)+(1-delta)*exp(k(-1));


a=rho1*a(-1)+ea;
p=rho2*p(-1)+ep;



end;
initval;
k=log(ks);
y=log(ys);
c=log(cs);
i=log(is);
a=log(as);
r=log(rs);
w=log(ws);
n=log(ns);
g=log(gs);
p=log(ps);
end;

shocks;
var ea = sda^2;
var ep = sda^2;
end;
resid(1);
steady;
check;

stoch_simul(order =1, hp_filter =1600,periods=1000);

きずな 发表于 2019-4-18 22:15
mod文件源码放上来看一下

var a n c k i r g w y m;
varexo ea ep;
parameters theta enta beta delta alpha1 alpha2 alpha3 sigma rho1 rho2 sda;
parameters ns cs ks as is rs ws gs ps;

alpha1 = 0.4;
alpha2 = 0.45;
alpha3 = 0.15;
beta  = 0.95;
delta = 0.11;
sigma = 1;
rho1 = 0.85;
rho2 = 0.75; %calibrated
theta = 0.84;
enta = 3;
sda = 0.009 %calibrated

as = 1;
ns = 1/3;
rs = 1/beta;
m = ns^(1+enta)/(1/beta+delta);
cs = (m*alpha1/alpha2)^(alpha1/alpha2)*(ns^(1+enta)*alpha3/alpha2)^(alpha3/alpha2)*ns;
ks = cs*m*alpha1/alpha2;
ys = ks*(1/beta+delta)/alpha1;
is=  delta*ks;
ws = alpha2*ys/ns;
gs = alpha3*ys/ps;


exp(n)^(enta)=alpha2*exp(y)/exp(n);


exp(p)=alpha3*exp(y)/exp(g);


1=beta*(exp(c(+1))^(-theta)/exp(c)^(-theta)*alpha1*exp(y(+1)/exp(k)-delta);


exp(r)=alpha1*exp(a)*exp(k(-1))^(alpha1)*exp(n)^(alpha2)*(exp(p)*exp(g))^(alpha3);




exp(w)=alpha2*exp(a)*exp(k(-1))^(alpha1)*exp(n)^(alpha2-1)*(exp(p)*exp(g))^(alpha3);


exp(y)=exp(c)+exp(i)+exp(p)*exp(g);
exp(y)=exp(a)*exp(k(-1))^(alpha1)*exp(n)^(alpha2)*(exp(p)*exp(g))^(alpha3);


exp(k)=exp(i)+(1-delta)*exp(k(-1));


a=rho1*a(-1)+ea;
p=rho2*p(-1)+ep;



end;
initval;
k=log(ks);
y=log(ys);
c=log(cs);
i=log(is);
a=log(as);
r=log(rs);
w=log(ws);
n=log(ns);
g=log(gs);
p=log(ps);
end;

shocks;
var ea = sda^2;
var ep = sda^2;
end;
resid(1);
steady;
check;


stoch_simul(order =1, hp_filter =1600,periods=1000);

きずな 发表于 2019-4-18 22:15
mod文件源码放上来看一下

1. var a n c k i r g w y m;
   
2. varexo ea ep;
   
3. parameters theta enta beta delta alpha1 alpha2 alpha3 sigma rho1 rho2 sda;
   
4. parameters ns cs ks as is rs ws gs ps;
   
5. 
   
6. alpha1 = 0.4;
   
7. alpha2 = 0.45;
   
8. alpha3 = 0.15;
   
9. beta  = 0.95;
   
10. delta = 0.11;
    
11. sigma = 1;
    
12. rho1 = 0.85;
    
13. rho2 = 0.75; %calibrated
    
14. theta = 0.84;
    
15. enta = 3;
    
16. sda = 0.009 %calibrated
    
17. 
    
18. as = 1;
    
19. ns = 1/3;
    
20. rs = 1/beta;
    
21. m = ns^(1+enta)/(1/beta+delta);
    
22. cs = (m*alpha1/alpha2)^(alpha1/alpha2)*(ns^(1+enta)*alpha3/alpha2)^(alpha3/alpha2)*ns;
    
23. ks = cs*m*alpha1/alpha2;
    
24. ys = ks*(1/beta+delta)/alpha1;
    
25. is=  delta*ks;
    
26. ws = alpha2*ys/ns;
    
27. gs = alpha3*ys/ps;
    
28. 
    
29. 
    
30. exp(n)^(enta)=alpha2*exp(y)/exp(n);
    
31. 
    
32. 
    
33. exp(p)=alpha3*exp(y)/exp(g);
    
34. 
    
35. 
    
36. 1=beta*(exp(c(+1))^(-theta)/exp(c)^(-theta)*alpha1*exp(y(+1)/exp(k)-delta);
    
37. 
    
38. 
    
39. exp(r)=alpha1*exp(a)*exp(k(-1))^(alpha1)*exp(n)^(alpha2)*(exp(p)*exp(g))^(alpha3);
    
40. 
    
41. 
    
42. 
    
43. exp(w)=alpha2*exp(a)*exp(k(-1))^(alpha1)*exp(n)^(alpha2-1)*(exp(p)*exp(g))^(alpha3);
    
44. 
    
45. 
    
46. exp(y)=exp(c)+exp(i)+exp(p)*exp(g);
    
47. exp(y)=exp(a)*exp(k(-1))^(alpha1)*exp(n)^(alpha2)*(exp(p)*exp(g))^(alpha3);
    
48. 
    
49. 
    
50. exp(k)=exp(i)+(1-delta)*exp(k(-1));
    
51. 
    
52. 
    
53. a=rho1*a(-1)+ea;
    
54. p=rho2*p(-1)+ep;
    
55. 
    
56. 
    
57. 
    
58. end;
    
59. initval;
    
60. k=log(ks);
    
61. y=log(ys);
    
62. c=log(cs);
    
63. i=log(is);
    
64. a=log(as);
    
65. r=log(rs);
    
66. w=log(ws);
    
67. n=log(ns);
    
68. g=log(gs);
    
69. p=log(ps);
    
70. end;
    
71. 
    
72. shocks;
    
73. var ea = sda^2;
    
74. var ep = sda^2;
    
75. end;
    
76. resid(1);
    
77. steady;
    
78. check;
    
79. 
    
80. stoch_simul(order =1, hp_filter =1600,periods=1000);

复制代码

きずな 发表于 2019-4-18 22:15
mod文件源码放上来看一下

var a n c k i r g w y m;
varexo ea ep;
parameters theta enta beta delta alpha1 alpha2 alpha3 sigma rho1 rho2 sda;
parameters ns cs ks as is rs ws gs ps;

alpha1 = 0.4;
alpha2 = 0.45;
alpha3 = 0.15;
beta  = 0.95;
delta = 0.11;
sigma = 1;
rho1 = 0.85;
rho2 = 0.75; %calibrated
theta = 0.84;
enta = 3;
sda = 0.009 %calibrated

as = 1;
ns = 1/3;
rs = 1/beta;
m = ns^(1+enta)/(1/beta+delta);
cs = (m*alpha1/alpha2)^(alpha1/alpha2)*(ns^(1+enta)*alpha3/alpha2)^(alpha3/alpha2)*ns;
ks = cs*m*alpha1/alpha2;
ys = ks*(1/beta+delta)/alpha1;
is=  delta*ks;
ws = alpha2*ys/ns;
gs = alpha3*ys/ps;


exp(n)^(enta)=alpha2*exp(y)/exp(n);


exp(p)=alpha3*exp(y)/exp(g);


1=beta*(exp(c(+1))^(-theta)/exp(c)^(-theta)*alpha1*exp(y(+1)/exp(k)-delta);


exp(r)=alpha1*exp(a)*exp(k(-1))^(alpha1)*exp(n)^(alpha2)*(exp(p)*exp(g))^(alpha3);



exp(w)=alpha2*exp(a)*exp(k(-1))^(alpha1)*exp(n)^(alpha2-1)*(exp(p)*exp(g))^(alpha3);


exp(y)=exp(c)+exp(i)+exp(p)*exp(g);
exp(y)=exp(a)*exp(k(-1))^(alpha1)*exp(n)^(alpha2)*(exp(p)*exp(g))^(alpha3);


exp(k)=exp(i)+(1-delta)*exp(k(-1));


a=rho1*a(-1)+ea;
p=rho2*p(-1)+ep;



end;
initval;
k=log(ks);
y=log(ys);
c=log(cs);
i=log(is);
a=log(as);
r=log(rs);
w=log(ws);
n=log(ns);
g=log(gs);
p=log(ps);
end;

shocks;
var ea = sda^2;
var ep = sda^2;
end;
resid(1);
steady;
check;

stoch_simul(order =1, hp_filter =1600,periods=1000);

这是原码希望大神指点，感谢

きずな 发表于 2019-4-18 22:15
mod文件源码放上来看一下

您好直接放代码一直在审核，我把代码文件放下面了希望得到您的指点

YYQX的鹤 发表于 2019-4-20 21:08
您好直接放代码一直在审核，我把代码文件放下面了希望得到您的指点

model模块没有声明，有一些参数没有定义，小错误自己仔细检查一下吧