//Simple MA(1) model

var y x;
varexo eps; 

parameters rho1; 
rho1  = 0.70;

model;
#rho2=1.5-rho1;
x      = eps;
y      = -rho1*x(-1) -rho2*x(-1) + rho2*rho1*x(-2) +eps;
end;

shocks;
var eps; stderr 1;
end;


steady;
check;

estimated_params;

//stderr eps,     INV_GAMMA_PDF,1,200;      //technology
stderr eps,     UNIFORM_PDF,,,0.001,5;      //technology

rho1,  ,            UNIFORM_PDF,,,     0,4;

end;


//stoch_simul(irf=10) x y;
varobs y;


estimation(datafile=ydataF,
                mode_compute=4,prefilter=0,mh_replic=5000 ,mh_nblocks=1,mh_jscale=0.50,mh_drop=0.2,
                mcmc_jumping_covariance =identity_matrix);