var y y_e i p_e p b y_g g t i_n;

varexo e_d e_s e_b e_t e_i e_pe e_ye e_yg e_g;

parameters sigma kapa beta fi_p fi_y tet_b tet_g tet_i tet_p tet_y ro_d ro_s ro_b ro_t ro_i;

fi_p=1.5;
fi_y=1.5;
tet_b=0.5;
tet_g=0.5;
tet_p=0.6;
tet_y=2.5;
tet_i=1.1;
beta=0.9;
kapa=2;
sigma=0.3;

ro_i=0.2;
ro_d=0.2;
ro_s=0.2;
ro_b=1;
ro_t=0.2;

model (linear);

y=y_e-sigma*(i-p_e)+e_d;
p=kapa*y+beta*y_e+e_s;
b=b(-1)+b(-1)*(i-p-y_g)+(g-t)+e_b;
i=i_n+fi_p*p+fi_y*y;
t=tet_b*b(-1)+tet_g*g+tet_p*p+tet_y*y+tet_i*i+e_t;

i_n=ro_i*i_n(-1)+e_i;

p_e=ro_d*p_e(-1)+e_pe;
y_e=ro_s*y_e(-1)+e_ye;
y_g=ro_b*y+p+e_yg;
g=ro_t*g(-1)+e_g;

end;

varobs y i p b t i_n p_e y_e g;

initval;
y_g=0.1;
end;

shocks;
var e_d; stderr 0.2;
var e_s; stderr 0.2;
var e_b; stderr 0.2;
var e_t; stderr 0.2;
var e_i; stderr 0.2;
var e_pe; stderr 0.2;
var e_ye; stderr 0.2;
var e_yg; stderr 0.2;
var e_g; stderr 0.2;
end;

check;
steady;

estimated_params;
sigma, uniform_pdf, 1,  0.7599, 0, 2;
kapa, uniform_pdf, 1, 0.7599, 0, 2;
beta, uniform_pdf, 0.5, 0.2886, 0, 1;
fi_p, uniform_pdf, 2.5, 1.2014, 0, 5;
fi_y, uniform_pdf, 2.5, 1.2014, 0, 5;
tet_b, uniform_pdf, 2.5, 1.2014, 0, 5;
tet_g, uniform_pdf, 2.5, 1.2014, 0, 5;
tet_i, uniform_pdf, 2.5, 1.2014, 0, 5;
tet_p, uniform_pdf, 2.5, 1.2014, 0, 5;
tet_y, uniform_pdf, 2.5, 1.2014, 0, 5;
ro_d, uniform_pdf, 0.5, 0.2886, 0, 1;
ro_s, uniform_pdf, 0.5, 0.2886, 0, 1;
ro_b, uniform_pdf, 0.5, 0.2886, 0, 1;
ro_t, uniform_pdf,  0.5, 0.2886, 0, 1;
ro_i,  uniform_pdf, 0.5, 0.2886, 0, 1;

stderr e_d, inv_gamma_pdf, 0.02, inf;
stderr e_s, inv_gamma_pdf, 0.02, inf;
stderr e_b, inv_gamma_pdf, 0.02, inf;
stderr e_t, inv_gamma_pdf, 0.02, inf;
stderr e_i, inv_gamma_pdf, 0.02, inf;
stderr e_pe, inv_gamma_pdf, 0.02, inf;
stderr e_ye, inv_gamma_pdf, 0.02, inf;
stderr e_yg, inv_gamma_pdf, 0.02, inf;
stderr e_g, inv_gamma_pdf, 0.02, inf;
end;

estimation(datafile=dadosdissertestac,nobs=1595,first_obs=1,mh_replic=2000,mh_nblocks=5,mh_drop=0.45,mh_jscale=0.4,mode_compute=6);