by hwangii » Mon Feb 22, 2010 1:20 pm
here is my .mod file:
periods 1000;
var w n t_n c k y x g d b R r lamda;
varexo e1 e2 e3 e4 e5;
parameters ros rod theta cdy xdy rhog beta_t bdy t_c t_k delta delta_t psi t_nbar z_bar eta beta kappa;
ros=0.7;
rod=0.38;
theta=0.62;
cdy=0.59;
xdy=0.27;
rhog=0.9;
beta_t=0.9;
bdy=2.52;
t_c=0.05;
t_k=0.36;
delta=0.02;
delta_t=0.025;
psi=0.05;
t_nbar=0.28;
z_bar=1.005;
eta=2;
beta=0.8;
kappa=0.6;
model;
w=(ros)*n+(1/(1-t_nbar))*t_n+c+e1;
w=(rod)*(k-n);
y=theta*n+(1-theta)*k+e2;
y=cdy*c+xdy*x+g;
g=rhog*g(-1)+e3;
d=g+(1/beta_t)*b(-1)+(bdy/beta_t)*(R-e2)-t_c*cdy*c-t_k*(1-theta-(delta/delta_t)*xdy)*k-t_k*(1-theta)*r;
psi*d-e4=theta*t_n+theta*t_nbar*(w+n);
(1-psi)*d+e4=b;
k=(1-delta_t)*k(-1)+delta_t*x(-1)-e2;
r=-(theta/(1-theta))*w;
R=(1-(1-(1-t_k)*delta)*beta*(z_bar^(-eta)))*r;
lamda=-eta*c-(1-eta)*kappa*n;
0=lamda(+1)-lamda+R(+1)+e5;
end;
initval;
w=0;
n=0;
t_n=0;
c=0;
k=0;
y=0;
x=0;
g=0;
d=0;
b=0;
R=0;
r=0;
lamda=0;
end;
shocks;
var e1;
stderr 1;
var e2;
stderr 1;
var e3;
stderr 1;
var e4;
stderr 1;
var e5;
stderr 1;
end;
steady;
stoch_simul(dr_algo=0, periods=1000, irf=40);
rplot w n t_n c k y x g d b R r lamda;
and i did not download this from the web--i created with in the m-file editor and save it as .mod file...i don't know a standard way to create .mod file so that's what i did. and the name of this .mod file is pset2.mod.