请教DSGE的程序问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

请教各位老师，dynare运行结果怎么都出不来呢？麻烦各位老师帮着看一下，谢谢各位老师了！
蓝字为我自己编的一个程序，红字为程序运行报错结果。


??? Error using ==> print_info at 51
One of the eigenvalues is close to 0/0 (the absolute value of numerator and denominator is smaller than
1e-6)


Error in ==> check at 51
    print_info(info, options_.noprint);


Error in ==> dierci at 215
check;


Error in ==> dynare at 120
evalin('base',fname) ;


var c h y i q r dp n m k w;varexo er;parameters beta1 beta2 sigma d1 phi x nu miu sigma_k theta rho_r theta_dp theta_y theta_q phi_r j dpbar mbar rbar cbar ybar nbar kbar hbar qbar;beta1=0.9;beta2=0.95;sigma=0.65;d1=0.7;phi=0.7;x=0.3;nu=0.3;miu=0.3;sigma_k=0.6;theta=0.75;rho_r=0.75; theta_dp=1.5;theta_y=0.6;theta_q=0.35; phi_r=0.7;j=0.65;dpbar=1;rbar=1/beta2;cbar=(1-beta1)^(1/sigma)/(1-j);mbar=(cbar-j*cbar)^sigma/(1-beta1);ybar=(1-beta2)*(cbar-j*cbar)^sigma/nu/x/(1-beta1+beta1*d1-beta2*d1);nbar=((1-nu-miu)*x*y/(cbar-j*cbar)^sigma)^(1/(1-phi));kbar=miu*x*ybar/(1-beta2*(1-sigma_k));hbar=(ybar/kbar^miu/nbar^(1-nu-miu))^(1/nu);qbar=nu*x*ybar/(1-beta2)/hbar;model(linear);w=phi*n+sigma*(c-j*c(-1));-h/hbar=qbar/(cbar-j*cbar)^sigma*(q-sigma*(c-j*c(-1)))-beta1*qbar/(cbar-j*cbar)^sigma*(q(+1)-sigma*(c(+1)-j*c))+beta1*d1*qbar/(cbar-j*cbar)^sigma*(q(+1)-sigma*(c(+1)-j*c))-d1*qbar*dpbar/rbar/(cbar-j*cbar)^sigma*(q(+1)+dp(+1)-r-sigma*(c(+1)-j*c));-m/mbar=(-sigma*(c-j*c(-1)))/(cbar-j*cbar)^sigma+beta1/(cbar-j*cbar)^sigma/dpbar*(-sigma*(c(+1)-j*c)-dp(+1));-c=r-c(+1)-dp(+1);miu*x*ybar/kbar*(y-c-k)=beta2*(1-sigma_k)*c(+1)-c;nu*x*ybar/hbar*(y-c-h)=qbar*(q-c)-beta2*qbar*(q(+1)-c(+1));(1-nu-miu)*x*(y-n)=w;y=0.0968*c+0.1452*i;r=rho_r*r(-1)+(1-rho_r)*(theta_dp*dp+theta_y*y+theta_q*q);dp=beta1*dp(+1)+(1-beta1*theta)*(1-theta)/theta*y;r=r(-1)*phi_r+er;end;initval;c=0;h=0;n=0;m=0;q=0;r=0;i=0;dp=0;w=0;k=0;y=0;er=0;end;steady;check;shocks;var er;stderr 0.1;end;estimated_params;phi_r,beta_pdf,0.7,0.1;stderr er,inv_gamma_pdf,0.1,0.1;end;varobs dp;estimation(datafile=cpi,prefilter=1,nobs=1,first_obs=197,mode_compute=6,mode_check,moments_varendo,order=1);stoch_simul(periods=2000,irf=100);

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：DSGE Eigenvalues Estimation eigenvalue parameters 程序

不知道怎么回事，程序都连在一起了，现在再发一遍程序，麻烦大家帮忙看看，谢谢各位老师
var c h y i q r dp n m k w;
varexo er;
parameters beta1 beta2 sigma d1 phi x nu miu sigma_k theta rho_r theta_dp theta_y theta_q phi_r j dpbar mbar rbar cbar ybar nbar kbar hbar qbar;
beta1=0.9;
beta2=0.95;
sigma=0.65;
d1=0.7;
phi=0.7;
x=0.3;
nu=0.3;
miu=0.3;
sigma_k=0.6;
theta=0.75;
rho_r=0.75;
theta_dp=1.5;
theta_y=0.6;
theta_q=0.35;
phi_r=0.7;
j=0.65;
dpbar=1;
rbar=1/beta2;
cbar=(1-beta1)^(1/sigma)/(1-j);
mbar=(cbar-j*cbar)^sigma/(1-beta1);
ybar=(1-beta2)*(cbar-j*cbar)^sigma/nu/x/(1-beta1+beta1*d1-beta2*d1);
nbar=((1-nu-miu)*x*y/(cbar-j*cbar)^sigma)^(1/(1-phi));
kbar=miu*x*ybar/(1-beta2*(1-sigma_k));
hbar=(ybar/kbar^miu/nbar^(1-nu-miu))^(1/nu);
qbar=nu*x*ybar/(1-beta2)/hbar;
model(linear);
w=phi*n+sigma*(c-j*c(-1));
-h/hbar=qbar/(cbar-j*cbar)^sigma*(q-sigma*(c-j*c(-1)))-beta1*qbar/(cbar-j*cbar)^sigma*(q(+1)-sigma*(c(+1)-j*c))+beta1*d1*qbar/(cbar-j*cbar)^sigma*(q(+1)-sigma*(c(+1)-j*c))-d1*qbar*dpbar/rbar/(cbar-j*cbar)^sigma*(q(+1)+dp(+1)-r-sigma*(c(+1)-j*c));
-m/mbar=(-sigma*(c-j*c(-1)))/(cbar-j*cbar)^sigma+beta1/(cbar-j*cbar)^sigma/dpbar*(-sigma*(c(+1)-j*c)-dp(+1));
-c=r-c(+1)-dp(+1);
miu*x*ybar/kbar*(y-c-k)=beta2*(1-sigma_k)*c(+1)-c;
nu*x*ybar/hbar*(y-c-h)=qbar*(q-c)-beta2*qbar*(q(+1)-c(+1));
(1-nu-miu)*x*(y-n)=w;
y=0.0968*c+0.1452*i;
r=rho_r*r(-1)+(1-rho_r)*(theta_dp*dp+theta_y*y+theta_q*q);
dp=beta1*dp(+1)+(1-beta1*theta)*(1-theta)/theta*y;
r=r(-1)*phi_r+er;
end;
initval;
c=0;
h=0;
n=0;
m=0;
q=0;
r=0;
i=0;
dp=0;
w=0;
k=0;
y=0;
er=0;
end;
steady;
check;
shocks;
var er;stderr 0.1;
end;
estimated_params;
phi_r,beta_pdf,0.7,0.1;
stderr er,inv_gamma_pdf,0.1,0.1;
end;
varobs dp;
estimation(datafile=cpi,prefilter=1,nobs=1,first_obs=197,mode_compute=6,mode_check,moments_varendo,order=1);
stoch_simul(periods=2000,irf=100);