DSGE模型讨论之七——最简单的DSGE模型的Dynare模拟和MLE，Bayesian估计

接下来按老师的模型，然后代入可观测数据，我用的是外汇储备和货币供给量的月度数据，经过HP滤波的处理，跳出来说模型 no stable equilibrium，这是怎么回事？没有稳定均衡是不是意味着我的模型建立的不好？这样的情况下，还有没有什么补救措施啊？
var y k l z c i g ex ms r f;
varexo eps_g eps_z eps_f eps_ms eps_ex;
parameters beta alpha delta gam_l gam_c gam_i gam_g gam_ex
           gam_ms j eta_ms  eta_g eta_f k2 h2
           rho_g rho_z rho_f rho_ex rho_fe rho_ef;
predetermined_variables k;         
           
beta=0.99;
alpha=0.9;
delta=0.05;
gam_l=1;
gam_c=0.3778;
gam_i=0.9;
gam_g=0.2;
gam_ex=0.5328;
gam_ms=0.01;
j=0.1;
eta_ms=0.6;
eta_g=0.7;
eta_f=2.2281;
k2=0.3;
h2=0.1;
rho_g=0.67;
rho_z=0.5;
rho_ex=0.2828;
rho_f=2;
rho_fe=-1.5;
rho_ef=1.1;
model(linear);
y=alpha*k+(1-alpha)*l+(log((1-beta+delta*beta)/(alpha*beta))+(1-alpha)*log(alpha/(1-alpha)))*z;
c=alpha*k-(alpha-gam_l)*l+(log((1-beta+delta*beta)/(alpha*beta))+(1-alpha)*log(alpha/(1-alpha)))*z;
c-c(+1)+(1-beta+delta*beta)*(1-alpha)*(1/(1-alpha)*log((1-beta+delta*beta)/(alpha*beta))+log(alpha/(1-alpha)))*z=0;
k(+1)=delta*i+(1-delta)*k;
i=-j*r;
y=gam_c*c+gam_i*i+gam_g*g+gam_ex*ex-gam_ms*(ms-ms(-1));
eta_ms*ms(-1)+eta_g*eps_g+eta_f*eps_f=k2/gam_ms*y-h2*r;
g=rho_g*g(-1)+eps_g;
z=rho_z*z(-1)+eps_z;
f=rho_f*f(-1)+rho_fe*ex(-1)+eps_f;
ex=rho_ex*ex(-1)+rho_ef(-1)*f+eps_ex;
end;
initval;
y=0;
k=0;
l=0;
z=0;
c=0;
i=0;
g=0;
ex=0;
ms=0;
r=0;
f=0;
end;
steady;
check;
shocks;
var eps_g=0.09^2;
var eps_z=0.09^2;
var eps_f=0.09^2;
var eps_ex=0.09^2;
var eps_ms=0.09^2;
end;
stoch_simul(order=1,nograph);

varobs ms f;
estimated_params;
beta, beta_pdf, 0.99, 0.002;
alpha, beta_pdf, 0.9, 0.002;
delta, beta_pdf, 0.05, 0.002;
gam_c, beta_pdf, 0.3778, 0.002;
gam_i, beta_pdf, 0.9, 0.002;
gam_g, beta_pdf, 0.2, 0.002;
gam_ex, beta_pdf, 0.5328, 0.002;
gam_ms, beta_pdf, 0.01, 0.002;
j, beta_pdf, 0.1, 0.002;
eta_ms, normal_pdf, 0.6, 0.002;
eta_g, normal_pdf, 0.7, 0.002;
eta_f, normal_pdf, 2.2281, 0.002;
k2, normal_pdf, 0.3, 0.002;
h2, gamma_pdf, 0.1,0.05;
rho_g, normal_pdf, 0.67, 0.002;
rho_z, normal_pdf, 0.5, 0.002;
rho_f, normal_pdf, 2, 0.002;
rho_ex, normal_pdf, 0.2828, 0.002;
rho_fe, normal_pdf, -1.5, 0.002;
rho_ef, normal_pdf, 1.1, 0.002;
end;
estimation(datafile=fms) f ms;

fms (2).xls (38.76 KB)

2012-5-24 22:52:45 上传

你好，我好好的看了下你的这个文件，心中困惑大了，你的这个mod里shocks;
var e; stderr 0.1;到底表示什么含义呢，按照你的文件来说是指第一期e的值突然由0变到0.1，以后各期e的值仍然为0，我以前的理解stderr 0.1为产生的关于e的随机数标准差增加为0.1，然后产生各期的e值，再根据各期e值来算各期的x和y.晕了。麻烦解答

hxw_0551 发表于 2012-5-25 20:05
你好，我好好的看了下你的这个文件，心中困惑大了，你的这个mod里shocks;
var e; stderr 0.1;到底表示什么 ...

那句命令就是制造一个标差为.1的白噪音。

老师，您有时间的时候，再帮我看看哈~到底问题在哪里，您给我指点一下，我自己会再好好琢磨修改。是不是在对数线性化的时候出了问题，导致模型有问题啊，我是参照toolkit的论文来进行线性化的，应该出不了问题啊。。exp（aXt+bYt）约等于aXt+bYt+1  然后就化为线性的了。simul那步之前都挺好的，varobs之后，时间序列再代入模型就出现问题了呢期待您的指导啊,您一直以来的帮助对我学习DSGE也是一种莫大的鼓励，再次向您致敬！我也一定要向您学习，成为一名优秀的经济科学家

tekuai5602 发表于 2012-5-28 10:19
老师，您有时间的时候，再帮我看看哈~到底问题在哪里，您给我指点一下，我自己会再好好琢磨修改。是不是在对 ...

早就说了，你这个模型中有几个方程是解不出来的，你所说的论文我也看了，根本和你的程序对不上。
变量C、MS都有问题。
C(t)=EtC(t+1)+....这个差分方程你自已动手解一下看看能不能解出解析解出来？
类似的还有MS。
再说了，经过滤波后的数据何必要求相邻两期的差分，经滤波后的数据应该是平稳的。

tekuai5602 发表于 2012-5-24 22:57
接下来按老师的模型，然后代入可观测数据，我用的是外汇储备和货币供给量的月度数据，经过HP滤波的处理，跳 ...

我之前改出来的那个mod文件，不能直接用来参数估计。我只是要告诉你一个意思，BK条件不满足的时候需要重新calibrate。你在写模型的时候，把k（如果这是capital的话）当成forward-looking variable，这即违反经济学原理有违反动态系统稳定性的。

我只是直接把它挑出来，用predetermined_variables k;    指明了这个backwards-looking。同时改了一些其他变量。你需要自己重新calibrate。这个非常有讲究，你自己试个上百遍就知道这件事多么tricky了。

mj2012 发表于 2012-5-28 11:47
早就说了，你这个模型中有几个方程是解不出来的，你所说的论文我也看了，根本和你的程序对不上。
变量C、 ...

那是不是就是我在线性化的时候出了差错，不然这个模型不应该对不上的吧。那我再试着重新编辑下模型，看看问题在哪里。。这个Et=Et+1我也是从欧拉公式线性化得到的。也许我哪里方法用的不对吧。我解了好几遍，最后得到的结果好像都和这个类似。我一会再重新推演下模型。谢谢老师指点。还有这个ms，那个万老师原文里是指y+ms-ms(-1)=c+i+g+ex作为行为人的效用最大化约束方程，我将它线性化后，不是ms-ms（-1）这项是不能避免的啊，老师说的ms是指这里么？老师看过那篇论文，那应该怎么解决？我的数据还需要HP去除趋势项么？

rastila 发表于 2012-5-28 15:22
我之前改出来的那个mod文件，不能直接用来参数估计。我只是要告诉你一个意思，BK条件不满足的时候需要重新 ...

恩，原来是这样啊，那我再反复的校正模型和估计参数试试，也许是我线性化的时候就出了什么问题。k应该是由看k(-1)与新增加投资决定的啊，我是这么定义的，这个怎么会变成forward-looking了？

tekuai5602 发表于 2012-5-28 17:05
恩，原来是这样啊，那我再反复的校正模型和估计参数试试，也许是我线性化的时候就出了什么问题。k应该是由 ...

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

上面equation就是你mod文件里面的。这个和线性化结果一般无太多关系，主要是parameterisation的问题。

非常感谢