贝叶斯估计出现的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

在做贝叶斯估计的时候，出现报错“Error using print_info (line 36)
The generalized Schur (QZ) decomposition failed. For more information, see the documentation for
Lapack function dgges: info=9, n=9
Error in check (line 76)
    print_info(info, options.noprint);
Error in chap4esti0112 (line 133)
oo_.dr.eigval = check(M_,options_,oo_);
Error in dynare (line 120)
evalin('base',fname) ;”

请教版上的大侠这是什么原因？
源程序如下：
var rc rd rl m dp c w h k z a g p;
varexo ze ae ge;
parameters beta delta theta b rhoz rhoa rhog alph;
predetermined_variables p;
model(linear);
w(+1)=w+beta*0.0351*rd;
p(+1)+c=z(+1)+p+w;
z+w-c-rd=0;
0.9004*c=1.8271*m(-1)-0.9268*dp;
1.8272*m+12.3965*k=2.3469*0.3263*(w(-1)+h)+0.0351*12.3965*(rc+k(-1))+(1-delta)*12.3965*k(-1)+1.0101*0.9268*(rd+dp);
1.2086*(a+theta*k(-1)+(1-theta)*h)=1.0101*2.3469*0.3263*(rl+w+h)+0.0351*12.3965*(rc+k(-1));
rl+w=a+theta*(k(-1)-h);
rc=a+(theta-1)*(k(-1)-h);
0.9268*(rl-rd)+1.8272*g=0;
0.9268*(dp-rl+rd)=2.3469*0.3263*(w+h)+alph*1.2086*(a+theta*k(-1)+(1-theta)*h);
z=rhoz*z(-1)+ze;
a=rhoa*a(-1)+ae;
g=rhog*g(-1)+ge;
end;
shocks;
var ze;stderr 0.099;
var ae;stderr 0.099;
var ge;stderr 0.099;
end;
varobs w dp;
check;
steady;
estimated_params;
beta,beta_pdf,0.99,0.002;
delta,beta_pdf,0.025,0.003;
theta,beta_pdf,0.36,0.002;
b,normal_pdf,-2.5805,0.002;
rhoz,beta_pdf,0.005,0.099;
rhoa,beta_pdf,0.005,0.099;
rhog,beta_pdf,0.005,0.099;
alph,beta_pdf,0.1319,0.002;
end;
estimation(datafile=estimation2,mh_replic=2000,mh_nblocks=2,mh_drop=0.45,mh_jscale=0.8);

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：贝叶斯估计 贝叶斯 composition information Generalized function failed 源程序 Error

大侠们啊快出来指点指点啊

放假了……

很有可能是prior的设置有问题

rastila 发表于 2014-1-17 17:08
很有可能是prior的设置有问题

那么一般先验分布该怎么设好呢？请传授下技巧吧，我用的参数都比价哦传统，如贴现率，折旧率，劳动资本比，这些都用什么分布比较好呢？

giant 发表于 2014-1-19 21:51
那么一般先验分布该怎么设好呢？请传授下技巧吧，我用的参数都比价哦传统，如贴现率，折旧率，劳动资本比 ...

可以参考下苗建军教授的著作，他比较喜欢用dynare，有介绍的。