楼主: giant
7401 9

[DSGE讨论专题] 用dynare编写程序时遇到个问题,请教 [推广有奖]

  • 1关注
  • 0粉丝

博士生

61%

还不是VIP/贵宾

-

威望
0
论坛币
95174 个
通用积分
5.1621
学术水平
2 点
热心指数
0 点
信用等级
0 点
经验
4606 点
帖子
155
精华
0
在线时间
453 小时
注册时间
2004-8-17
最后登录
2018-12-8

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

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

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
程序报错如下:“Error using print_info (line 40)
Blanchard Kahn conditions are not satisfied: no stable equilibrium
Error in stoch_simul (line 81)
    print_info(info, options_.noprint);
Error in chap5 (line 159)
info = stoch_simul(var_list_);
Error in dynare (line 120)
evalin('base',fname) ;”
源程序如下,请大侠帮忙指点一下是什么地方出问题了?

var rc rd rl m dp c w h k z a g p gg;
varexo ze ae ge gge;
parameters beta delta theta b rhoz rhoa rhog tao rhogg;
beta=0.99;
delta=0.05;
theta=0.36;
b=-2.5805;
rhoz=0.005;
rhoa=0.005;
rhog=0.005;
tao=0.0609;
rhogg=0.005;
predetermined_variables p;
model(linear);
1+z+p(+1)+c(+1)=beta*1.0101*(1+rd+z(+1)+p+c);
1+z+rd(+1)+c(+1)=beta*(1+z(+1)+rd)*(0.0351*(1+rc)+(1-delta));
b*0.7901*1.0101*(1+p+c+rd)=2.3469*(1+z+w)*(tao*1.0101*(1+rd)-(1-tao)*(1+p));
0.7901*(1+c)=3.1586*(1+m(-1))-2.2326*(1+dp)-tao*2.3469*0.9513*(1+w+h);
3.1586*(1+m)+36.1408*(1+k)=(1-tao)*2.3469*0.9513*(1+w+h)+0.0351*36.1408*(1+rc+k(-1))+(1-delta)*36.1408*(1+k(-1))+1.0101*2.2326*(1+rd+dp);
3.5237*(a+theta*k(-1)+(1-theta)*h)=1.0101*2.3469*0.9513*(rl+w+h)+0.0351*36.1408*(rc+k(-1));
rl+w=a+theta*(k(-1)-h);
rc=a+(theta-1)*(k(-1)-h);
2.2326*(rl-rd)+3.1586*g=0;
rl-rd=dp-gg;
z=rhoz*z(-1)+ze;
a=rhoa*a(-1)+ae;
g=rhog*g(-1)+ge;
gg=rhogg*gg(-1)+gge;
end;
check;
steady;
shocks;
var ze;stderr 0.099;
var ae;stderr 0.099;
var ge;stderr 0.099;
end;
stoch_simul(order=1,irf=30);

二维码

扫码加我 拉你入群

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

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

关键词:dynare 编写程序 ARE equilibrium Determined 编写程序

已有 1 人评分经验 论坛币 收起 理由
linmengmiki + 60 + 60 鼓励积极发帖讨论

总评分: 经验 + 60  论坛币 + 60   查看全部评分

沙发
richardgu26 发表于 2013-12-16 08:46:44 |只看作者 |坛友微信交流群
你的变量p是什么,一般来说,k是资本,这个是predetermined_variables。

使用道具

藤椅
giant 发表于 2013-12-31 22:29:10 |只看作者 |坛友微信交流群
p是价格,模型里价格本不是必须的,但修改后似乎价格无法消除,不知道如何处理好?

使用道具

板凳
giant 发表于 2014-1-1 14:22:10 |只看作者 |坛友微信交流群
我修改了模型,取消了价格,程序如下,但仍然提示有问题,请大侠帮忙指点一下:“Error using print_info (line 43)
Blanchard Kahn conditions are not satisfied: indeterminacy

Error in stoch_simul (line 81)
    print_info(info, options_.noprint);

Error in chap50101 (line 155)
info = stoch_simul(var_list_);

Error in dynare (line 120)
evalin('base',fname) ;”

var rc rd rl m dp c w h k z a g gg;
varexo ze ae ge gge;
parameters beta delta theta b rhoz rhoa rhog tao rhogg;
beta=0.99;
delta=0.05;
theta=0.36;
b=-2.5805;
rhoz=0.005;
rhoa=0.005;
rhog=0.005;
tao=0.0609;
rhogg=0.005;

model(linear);
1+z+c(+1)=beta*1.0101*(1+rd+z(+1)+c);
1+z+rd(+1)+c(+1)=beta*(1+z(+1)+rd)*(0.0351*(1+rc)+(1-delta));
b*0.7901*1.0101*(1+c+rd)=2.3469*(1+z+w)*(tao*1.0101*(1+rd)-(1-tao));
0.7901*(1+c)=3.1586*(1+m(-1))-2.2326*(1+dp)-tao*2.3469*0.9513*(1+w+h);
3.1586*(1+m)+36.1408*(1+k)=(1-tao)*2.3469*0.9513*(1+w+h)+0.0351*36.1408*(1+rc+k(-1))+(1-delta)*36.1408*(1+k(-1))+1.0101*2.2326*(1+rd+dp);
3.5237*(a+theta*k(-1)+(1-theta)*h)=1.0101*2.3469*0.9513*(rl+w+h)+0.0351*36.1408*(rc+k(-1));
rl+w=a+theta*(k(-1)-h);
rc=a+(theta-1)*(k(-1)-h);
2.2326*(dp-w-h)+3.1586*g=0;
z=rhoz*z(-1)+ze;
a=rhoa*a(-1)+ae;
g=rhog*g(-1)+ge;
gg=rhogg*gg(-1)+gge;
end;

check;
steady;
shocks;
var ze;stderr 0.099;
var ae;stderr 0.099;
var ge;stderr 0.099;
end;
stoch_simul(order=1,irf=30);

使用道具

报纸
ahnulxy 发表于 2014-2-17 10:17:33 |只看作者 |坛友微信交流群
giant 发表于 2014-1-1 14:22
我修改了模型,取消了价格,程序如下,但仍然提示有问题,请大侠帮忙指点一下:“Error using print_info ( ...
真的不建议模型直接写出linear形式,原因如下
1. 需要手动线性化,耗时费力,容易出错
2. 一旦出错,很难排查
3. 模型一旦线性化做高阶分析,就不可能,尤其做福利分析

所以建议写出非线性形式,自己手动求解steady state,让dynare自己线性化或者高阶近似
琅琅教育 www.longlongedu.com  
人生最大的杯具:白天看中国股市,晚上看中国足球!

使用道具

地板
gssdzc 在职认证  发表于 2014-4-10 18:16:40 |只看作者 |坛友微信交流群
这个建议好!

使用道具

7
信陵公子 发表于 2019-3-3 19:50:06 |只看作者 |坛友微信交流群
ahnulxy 发表于 2014-2-17 10:17
真的不建议模型直接写出linear形式,原因如下
1. 需要手动线性化,耗时费力,容易出错
2. 一旦出错,很 ...
您好,您说的不线性化,如何处理calvo定价的问题呢?

使用道具

8
ahnulxy 发表于 2019-3-4 09:26:28 |只看作者 |坛友微信交流群
写成递归形式啊,请参考 https://item.jd.com/12511678.html, 有非常详细地解释
已有 1 人评分经验 论坛币 收起 理由
linmengmiki + 60 + 60 热心帮助其他会员

总评分: 经验 + 60  论坛币 + 60   查看全部评分

使用道具

9
tianwk 发表于 2019-3-4 16:31:54 |只看作者 |坛友微信交流群
thanks for sharing

使用道具

10
大侠无想 学生认证  发表于 2020-4-12 11:55:28 |只看作者 |坛友微信交流群
代码出现如下错误提示,请问通常是出现了什么问题?

错误使用 print_info (line 32)
The steadystate file did not compute the steady state
出错 stoch_simul (line 103)
    print_info(info, options_.noprint, options_);
出错 HFA2.driver (line 608)
[info, oo_, options_, M_] = stoch_simul(M_, options_, oo_, var_list_);
出错 dynare (line 293)
evalin('base',[fname '.driver']) ;

使用道具

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

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

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

GMT+8, 2024-4-20 02:07