全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 MATLAB等数学软件专版
4906 3
2014-11-28
There are 9 eigenvalue(s) larger than 1 in modulus
for 7 forward-looking variable(s)

The rank condition ISN'T verified!


You did not declare endogenous variables after the estimation/calib_smoother command.

Loading 52 observations from datass.xls


ESTIMATION_CHECKS: There was an error in computing the likelihood for initial parameter values.
ESTIMATION_CHECKS: You should try using the calibrated version of the model as starting values. To do
ESTIMATION_CHECKS: this, add an empty estimated_params_init-block with use_calibration option immediately before the estimation
ESTIMATION_CHECKS: command (and after the estimated_params-block so that it does not get overwritten):


Error using initial_estimation_checks (line 34)
initial_estimation_checks:: Estimation can't take place because there are less declared
shocks than observed variables!

Error in initial_estimation_checks (line 34)
    error(['initial_estimation_checks:: Estimation can''t take place because there are
    less declared shocks than observed variables!'])

Error in dynare_estimation_1 (line 179)
    oo_ =
    initial_estimation_checks(objective_function,xparam1,dataset_,M_,estim_params_,options_,bayestopt_,oo_);

Error in dynare_estimation (line 89)
    dynare_estimation_1(var_list,dname);

Error in miao2 (line 426)
dynare_estimation(var_list_);

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

这是我的代码!
var y pi m r i oh c q of h w rk k u x a z p gama_w gama_p v;
varexo e_z e_x e_a e_p e_gama_w e_gama_p e_v;
parameters sigma_m sigma_c kesi beta sigma fy xi_w sigma_L alpha_1 alpha_2
alpha_3 delta xi_p omiga_r omiga_pi omiga_y omiga_u ep_c ep_i ep_o ro_a ro_x
ro_z ro_p ro_gama_p ro_gama_w ro_v sigma_a sigma_x sigma_gama_w sigma_p sigma_z
sigma_gama_p sigma_v ;

sigma_m=1.4;
sigma_c=2;
kesi=0.98;
beta=0.992;
sigma=0.03;
fy=4;
xi_w=0.5;
sigma_L=1;
alpha_1=0.6;
alpha_2=0.21;
alpha_3=0.19;
delta=0.5;
xi_p=0.5;
omiga_u=1;
omiga_pi=2;
omiga_y=2;
omiga_r=1;
ep_c=0.3649;
ep_i=0.4729;
ep_o=0.1622;
ro_z=0.5;
ro_v=0.5;
ro_a=0.5;
ro_x=0.5;
ro_p=0.5;
ro_gama_w=0.5;
ro_gama_p=0.5;
sigma_a=0.1;
sigma_x=0.1;
sigma_gama_w=0.1;
sigma_p=0.1;
sigma_z=0.1;
sigma_gama_p=0.1;
sigma_v=0.1;

model(linear);
sigma_m*m=(sigma_c-sigma_c*kesi)*oh+(sigma_c*kesi-kesi+1)*c-beta/(1-beta)*r+beta/(1-beta)*pi(+1);
p=(1-kesi)*c-kesi*oh(-1);
r=pi(+1)+a-a(+1)+(sigma_c*kesi-sigma_c)*(oh-oh(-1))+(kesi-sigma_c*kesi-1)*(c-c(-1));
q=beta*(1-sigma)*q(+1)-r+pi(+1)+(1-beta*(1-sigma))*rk(+1);
i=beta/(1+beta)*i(+1)+1/(1+beta)*i+fy/(1+beta)*q+1/(1+beta)*(beta*x(+1)-x);
k=(1-sigma)*k(-1)+sigma*i;
w=k(-1)+rk-h;
rk=of+p-k(-1);
w=1/(1+beta)*(pi(-1)+w(-1))-pi+beta/(1+beta)*(pi(+1)+w(+1))+(1-beta*xi_w)*(1-xi_w)*(-w+sigma_L*h+(sigma_c-sigma_c*kesi)*oh+(sigma_c*kesi-kesi+1)*c+gama_w);
y=z+alpha_1*h+alpha_2*k(-1)+alpha_3*of;
pi=(1-delta)*pi(-1)+beta*(pi(+1)-(1-delta)*pi)+(1-beta*xi_p)*(1-xi_p)/xi_p*(-z+alpha_1*w+alpha_2*rk+alpha_3*p+gama_p);
u=m-m(-1)+pi;
r=omiga_r*r(-1)+(1-omiga_r)*(omiga_pi*pi+omiga_y*y+omiga_u*u)+v;
y=ep_c*c+ep_i*i+ep_o*(oh+of+p);
a=sigma_a*e_a+ro_a*a(-1);
x=sigma_x*e_x+ro_x*x(-1);
z=sigma_z*e_z+ro_z*z(-1);
p=sigma_p*e_p+ro_p*p(-1);
gama_p=sigma_gama_p*e_gama_p+ro_gama_p*gama_p(-1);
gama_w=sigma_gama_w*e_gama_w+ro_gama_w*gama_w(-1);
v=sigma_v*e_v+ro_v*v(-1);
end;

initval;
y=0.0;
pi=0.0;
m=0.0;
r=0.0;
i=0.0;
oh=0.0;
c=0.0;
q=0.0;
of=0.0;
h=0.0;
w=0.0;
rk=0.0;
k=0.0;
u=0.0;
x=0.0;
a=0.0;
z=0.0;
p=0.0;
v=0.0;
gama_w=0.0;
gama_p=0.0;
end;

steady;
check;

shocks;
var e_p;
stderr 1;
var e_z;
stderr 1;
var e_a;
stderr 1;
var e_x;
stderr 1;
var e_gama_p;
stderr 1;
Var e_gama_w;
stderr 1;
var e_v;
stderr 1;
end;

estimated_params;
sigma_c,normal_pdf,2,0.25;
kesi,normal_pdf,0.98,0.01;
fy,normal_pdf,4,1.5;
xi_p,beta_pdf,0.5,0.1;
xi_w,beta_pdf,0.5,0.1;
delta,beta_pdf,0.5,0.015;
omiga_pi,normal_pdf,-0.5,0.1;
omiga_y,normal_pdf,-0.5,0.1;
omiga_u,normal_pdf,-0.5,0.1;
ro_a,beta_pdf,0.5,0.2;
ro_x,beta_pdf,0.5,0.2;
ro_z,beta_pdf,0.5,0.2;
ro_p,beta_pdf,0.5,0.2;
ro_gama_p,beta_pdf,0.5,0.2;
ro_gama_w,beta_pdf,0.5,0.2;
ro_v,beta_pdf,0.5,0.2;
sigma_a,inv_gamma_pdf,0.1,2;
sigma_x,inv_gamma_pdf,0.1,2;
sigma_z,inv_gamma_pdf,0.1,2;
sigma_p,inv_gamma_pdf,0.1,2;
sigma_v,inv_gamma_pdf,0.1,2;
sigma_gama_p,inv_gamma_pdf,0.1,2;
sigma_gama_w,inv_gamma_pdf,0.1,2;
end;
varobs y pi m r i oh c q of h w ;
estimation(datafile=datass,mh_replic=10000);
stoch_simul(periods=2100);


附件列表

miao2.zip

大小:1.15 KB

 马上下载

本附件包括:

  • miao2.mod

二维码

扫码加我 拉你入群

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

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

全部回复
2014-11-29 07:52:13
第一个可能的问题,你的参数校准可能有问题。
第二个问题是,你的观测方程呢?好像没有看到啊。如果是这样的话,你的数据怎么和你的DSGE模型结合呢?
二维码

扫码加我 拉你入群

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

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

2014-12-1 09:19:45
richardgu26 发表于 2014-11-29 07:52
第一个可能的问题,你的参数校准可能有问题。
第二个问题是,你的观测方程呢?好像没有看到啊。如果是这样 ...
参数校准需要怎么修改呢?因为我也是才学这个,不是很懂,这些都是模仿别人论文编的,真心请您指导,万分感谢!
二维码

扫码加我 拉你入群

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

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

2014-12-2 15:49:19
真心求指教,因为毕业需要用这个模型,请懂的不吝赐教,如果能够搞出来一定重谢!报酬可以私信商量!
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群