% % Status : main Dynare file % % Warning : this file is generated automatically by Dynare % from model file (.mod) clear all tic; global M_ oo_ options_ ys0_ ex0_ estimation_info options_ = []; M_.fname = 'mydata1'; % % Some global variables initialization % global_initialization; diary off; diary('mydata1.log'); M_.exo_names = 'ETA_B'; M_.exo_names_tex = 'ETA\_B'; M_.exo_names_long = 'ETA_B'; M_.exo_names = char(M_.exo_names, 'ETA_P'); M_.exo_names_tex = char(M_.exo_names_tex, 'ETA\_P'); M_.exo_names_long = char(M_.exo_names_long, 'ETA_P'); M_.exo_names = char(M_.exo_names, 'ETA_Y'); M_.exo_names_tex = char(M_.exo_names_tex, 'ETA\_Y'); M_.exo_names_long = char(M_.exo_names_long, 'ETA_Y'); M_.exo_names = char(M_.exo_names, 'ETA_X'); M_.exo_names_tex = char(M_.exo_names_tex, 'ETA\_X'); M_.exo_names_long = char(M_.exo_names_long, 'ETA_X'); M_.exo_names = char(M_.exo_names, 'ETA_T'); M_.exo_names_tex = char(M_.exo_names_tex, 'ETA\_T'); M_.exo_names_long = char(M_.exo_names_long, 'ETA_T'); M_.exo_names = char(M_.exo_names, 'ETA_G'); M_.exo_names_tex = char(M_.exo_names_tex, 'ETA\_G'); M_.exo_names_long = char(M_.exo_names_long, 'ETA_G'); M_.exo_names = char(M_.exo_names, 'ETA_R'); M_.exo_names_tex = char(M_.exo_names_tex, 'ETA\_R'); M_.exo_names_long = char(M_.exo_names_long, 'ETA_R'); M_.endo_names = 'dYobs'; M_.endo_names_tex = 'dYobs'; M_.endo_names_long = 'dYobs'; M_.endo_names = char(M_.endo_names, 'dCobs'); M_.endo_names_tex = char(M_.endo_names_tex, 'dCobs'); M_.endo_names_long = char(M_.endo_names_long, 'dCobs'); M_.endo_names = char(M_.endo_names, 'dIobs'); M_.endo_names_tex = char(M_.endo_names_tex, 'dIobs'); M_.endo_names_long = char(M_.endo_names_long, 'dIobs'); M_.endo_names = char(M_.endo_names, 'dTBYobs'); M_.endo_names_tex = char(M_.endo_names_tex, 'dTBYobs'); M_.endo_names_long = char(M_.endo_names_long, 'dTBYobs'); M_.endo_names = char(M_.endo_names, 'PIobs'); M_.endo_names_tex = char(M_.endo_names_tex, 'PIobs'); M_.endo_names_long = char(M_.endo_names_long, 'PIobs'); M_.endo_names = char(M_.endo_names, 'iobs'); M_.endo_names_tex = char(M_.endo_names_tex, 'iobs'); M_.endo_names_long = char(M_.endo_names_long, 'iobs'); M_.endo_names = char(M_.endo_names, 'dSobs'); M_.endo_names_tex = char(M_.endo_names_tex, 'dSobs'); M_.endo_names_long = char(M_.endo_names_long, 'dSobs'); M_.endo_names = char(M_.endo_names, 'Cn'); M_.endo_names_tex = char(M_.endo_names_tex, 'Cn'); M_.endo_names_long = char(M_.endo_names_long, 'Cn'); M_.endo_names = char(M_.endo_names, 'Cm'); M_.endo_names_tex = char(M_.endo_names_tex, 'Cm'); M_.endo_names_long = char(M_.endo_names_long, 'Cm'); M_.endo_names = char(M_.endo_names, 'P'); M_.endo_names_tex = char(M_.endo_names_tex, 'P'); M_.endo_names_long = char(M_.endo_names_long, 'P'); M_.endo_names = char(M_.endo_names, 'PI'); M_.endo_names_tex = char(M_.endo_names_tex, 'PI'); M_.endo_names_long = char(M_.endo_names_long, 'PI'); M_.endo_names = char(M_.endo_names, 'C'); M_.endo_names_tex = char(M_.endo_names_tex, 'C'); M_.endo_names_long = char(M_.endo_names_long, 'C'); M_.endo_names = char(M_.endo_names, 'D'); M_.endo_names_tex = char(M_.endo_names_tex, 'D'); M_.endo_names_long = char(M_.endo_names_long, 'D'); M_.endo_names = char(M_.endo_names, 'H'); M_.endo_names_tex = char(M_.endo_names_tex, 'H'); M_.endo_names_long = char(M_.endo_names_long, 'H'); M_.endo_names = char(M_.endo_names, 'Y'); M_.endo_names_tex = char(M_.endo_names_tex, 'Y'); M_.endo_names_long = char(M_.endo_names_long, 'Y'); M_.endo_names = char(M_.endo_names, 'Rk'); M_.endo_names_tex = char(M_.endo_names_tex, 'Rk'); M_.endo_names_long = char(M_.endo_names_long, 'Rk'); M_.endo_names = char(M_.endo_names, 'W'); M_.endo_names_tex = char(M_.endo_names_tex, 'W'); M_.endo_names_long = char(M_.endo_names_long, 'W'); M_.endo_names = char(M_.endo_names, 'We'); M_.endo_names_tex = char(M_.endo_names_tex, 'We'); M_.endo_names_long = char(M_.endo_names_long, 'We'); M_.endo_names = char(M_.endo_names, 'MC'); M_.endo_names_tex = char(M_.endo_names_tex, 'MC'); M_.endo_names_long = char(M_.endo_names_long, 'MC'); M_.endo_names = char(M_.endo_names, 'Pn'); M_.endo_names_tex = char(M_.endo_names_tex, 'Pn'); M_.endo_names_long = char(M_.endo_names_long, 'Pn'); M_.endo_names = char(M_.endo_names, 'PIn'); M_.endo_names_tex = char(M_.endo_names_tex, 'PIn'); M_.endo_names_long = char(M_.endo_names_long, 'PIn'); M_.endo_names = char(M_.endo_names, 'K'); M_.endo_names_tex = char(M_.endo_names_tex, 'K'); M_.endo_names_long = char(M_.endo_names_long, 'K'); M_.endo_names = char(M_.endo_names, 'Q'); M_.endo_names_tex = char(M_.endo_names_tex, 'Q'); M_.endo_names_long = char(M_.endo_names_long, 'Q'); M_.endo_names = char(M_.endo_names, 'I'); M_.endo_names_tex = char(M_.endo_names_tex, 'I'); M_.endo_names_long = char(M_.endo_names_long, 'I'); M_.endo_names = char(M_.endo_names, 'u'); M_.endo_names_tex = char(M_.endo_names_tex, 'u'); M_.endo_names_long = char(M_.endo_names_long, 'u'); M_.endo_names = char(M_.endo_names, 'De'); M_.endo_names_tex = char(M_.endo_names_tex, 'De'); M_.endo_names_long = char(M_.endo_names_long, 'De'); M_.endo_names = char(M_.endo_names, 'wbar'); M_.endo_names_tex = char(M_.endo_names_tex, 'wbar'); M_.endo_names_long = char(M_.endo_names_long, 'wbar'); M_.endo_names = char(M_.endo_names, 'Ce'); M_.endo_names_tex = char(M_.endo_names_tex, 'Ce'); M_.endo_names_long = char(M_.endo_names_long, 'Ce'); M_.endo_names = char(M_.endo_names, 'NW'); M_.endo_names_tex = char(M_.endo_names_tex, 'NW'); M_.endo_names_long = char(M_.endo_names_long, 'NW'); M_.endo_names = char(M_.endo_names, 'Re'); M_.endo_names_tex = char(M_.endo_names_tex, 'Re'); M_.endo_names_long = char(M_.endo_names_long, 'Re'); M_.endo_names = char(M_.endo_names, 'i'); M_.endo_names_tex = char(M_.endo_names_tex, 'i'); M_.endo_names_long = char(M_.endo_names_long, 'i'); M_.endo_names = char(M_.endo_names, 'X'); M_.endo_names_tex = char(M_.endo_names_tex, 'X'); M_.endo_names_long = char(M_.endo_names_long, 'X'); M_.endo_names = char(M_.endo_names, 'Pm'); M_.endo_names_tex = char(M_.endo_names_tex, 'Pm'); M_.endo_names_long = char(M_.endo_names_long, 'Pm'); M_.endo_names = char(M_.endo_names, 'PIm'); M_.endo_names_tex = char(M_.endo_names_tex, 'PIm'); M_.endo_names_long = char(M_.endo_names_long, 'PIm'); M_.endo_names = char(M_.endo_names, 'MCm'); M_.endo_names_tex = char(M_.endo_names_tex, 'MCm'); M_.endo_names_long = char(M_.endo_names_long, 'MCm'); M_.endo_names = char(M_.endo_names, 'Ym'); M_.endo_names_tex = char(M_.endo_names_tex, 'Ym'); M_.endo_names_long = char(M_.endo_names_long, 'Ym'); M_.endo_names = char(M_.endo_names, 'S'); M_.endo_names_tex = char(M_.endo_names_tex, 'S'); M_.endo_names_long = char(M_.endo_names_long, 'S'); M_.endo_names = char(M_.endo_names, 'E_B'); M_.endo_names_tex = char(M_.endo_names_tex, 'E\_B'); M_.endo_names_long = char(M_.endo_names_long, 'E_B'); M_.endo_names = char(M_.endo_names, 'E_P'); M_.endo_names_tex = char(M_.endo_names_tex, 'E\_P'); M_.endo_names_long = char(M_.endo_names_long, 'E_P'); M_.endo_names = char(M_.endo_names, 'E_Y'); M_.endo_names_tex = char(M_.endo_names_tex, 'E\_Y'); M_.endo_names_long = char(M_.endo_names_long, 'E_Y'); M_.endo_names = char(M_.endo_names, 'E_G'); M_.endo_names_tex = char(M_.endo_names_tex, 'E\_G'); M_.endo_names_long = char(M_.endo_names_long, 'E_G'); M_.endo_names = char(M_.endo_names, 'istar'); M_.endo_names_tex = char(M_.endo_names_tex, 'istar'); M_.endo_names_long = char(M_.endo_names_long, 'istar'); M_.endo_names = char(M_.endo_names, 'Pstar'); M_.endo_names_tex = char(M_.endo_names_tex, 'Pstar'); M_.endo_names_long = char(M_.endo_names_long, 'Pstar'); M_.endo_names = char(M_.endo_names, 'Cstar'); M_.endo_names_tex = char(M_.endo_names_tex, 'Cstar'); M_.endo_names_long = char(M_.endo_names_long, 'Cstar'); M_.endo_names = char(M_.endo_names, 'BMA'); M_.endo_names_tex = char(M_.endo_names_tex, 'BMA'); M_.endo_names_long = char(M_.endo_names_long, 'BMA'); M_.endo_names = char(M_.endo_names, 'IMA'); M_.endo_names_tex = char(M_.endo_names_tex, 'IMA'); M_.endo_names_long = char(M_.endo_names_long, 'IMA'); M_.endo_names = char(M_.endo_names, 'TMA'); M_.endo_names_tex = char(M_.endo_names_tex, 'TMA'); M_.endo_names_long = char(M_.endo_names_long, 'TMA'); M_.param_names = 'alpha'; M_.param_names_tex = 'alpha'; M_.param_names_long = 'alpha'; M_.param_names = char(M_.param_names, 'beta'); M_.param_names_tex = char(M_.param_names_tex, 'beta'); M_.param_names_long = char(M_.param_names_long, 'beta'); M_.param_names = char(M_.param_names, 'delta'); M_.param_names_tex = char(M_.param_names_tex, 'delta'); M_.param_names_long = char(M_.param_names_long, 'delta'); M_.param_names = char(M_.param_names, 'gamma'); M_.param_names_tex = char(M_.param_names_tex, 'gamma'); M_.param_names_long = char(M_.param_names_long, 'gamma'); M_.param_names = char(M_.param_names, 'theta'); M_.param_names_tex = char(M_.param_names_tex, 'theta'); M_.param_names_long = char(M_.param_names_long, 'theta'); M_.param_names = char(M_.param_names, 'omega'); M_.param_names_tex = char(M_.param_names_tex, 'omega'); M_.param_names_long = char(M_.param_names_long, 'omega'); M_.param_names = char(M_.param_names, 'nu'); M_.param_names_tex = char(M_.param_names_tex, 'nu'); M_.param_names_long = char(M_.param_names_long, 'nu'); M_.param_names = char(M_.param_names, 'epsilon_w'); M_.param_names_tex = char(M_.param_names_tex, 'epsilon\_w'); M_.param_names_long = char(M_.param_names_long, 'epsilon_w'); M_.param_names = char(M_.param_names, 'phi'); M_.param_names_tex = char(M_.param_names_tex, 'phi'); M_.param_names_long = char(M_.param_names_long, 'phi'); M_.param_names = char(M_.param_names, 'eta'); M_.param_names_tex = char(M_.param_names_tex, 'eta'); M_.param_names_long = char(M_.param_names_long, 'eta'); M_.param_names = char(M_.param_names, 'hab'); M_.param_names_tex = char(M_.param_names_tex, 'hab'); M_.param_names_long = char(M_.param_names_long, 'hab'); M_.param_names = char(M_.param_names, 'phi_c'); M_.param_names_tex = char(M_.param_names_tex, 'phi\_c'); M_.param_names_long = char(M_.param_names_long, 'phi_c'); M_.param_names = char(M_.param_names, 'psi_D'); M_.param_names_tex = char(M_.param_names_tex, 'psi\_D'); M_.param_names_long = char(M_.param_names_long, 'psi_D'); M_.param_names = char(M_.param_names, 'l_n'); M_.param_names_tex = char(M_.param_names_tex, 'l\_n'); M_.param_names_long = char(M_.param_names_long, 'l_n'); M_.param_names = char(M_.param_names, 'xi_n'); M_.param_names_tex = char(M_.param_names_tex, 'xi\_n'); M_.param_names_long = char(M_.param_names_long, 'xi_n'); M_.param_names = char(M_.param_names, 'l_m'); M_.param_names_tex = char(M_.param_names_tex, 'l\_m'); M_.param_names_long = char(M_.param_names_long, 'l_m'); M_.param_names = char(M_.param_names, 'xi_m'); M_.param_names_tex = char(M_.param_names_tex, 'xi\_m'); M_.param_names_long = char(M_.param_names_long, 'xi_m'); M_.param_names = char(M_.param_names, 'l_w'); M_.param_names_tex = char(M_.param_names_tex, 'l\_w'); M_.param_names_long = char(M_.param_names_long, 'l_w'); M_.param_names = char(M_.param_names, 'xi_w'); M_.param_names_tex = char(M_.param_names_tex, 'xi\_w'); M_.param_names_long = char(M_.param_names_long, 'xi_w'); M_.param_names = char(M_.param_names, 'chi'); M_.param_names_tex = char(M_.param_names_tex, 'chi'); M_.param_names_long = char(M_.param_names_long, 'chi'); M_.param_names = char(M_.param_names, 'zeta_u'); M_.param_names_tex = char(M_.param_names_tex, 'zeta\_u'); M_.param_names_long = char(M_.param_names_long, 'zeta_u'); M_.param_names = char(M_.param_names, 'zeta_e'); M_.param_names_tex = char(M_.param_names_tex, 'zeta\_e'); M_.param_names_long = char(M_.param_names_long, 'zeta_e'); M_.param_names = char(M_.param_names, 'zeta_l'); M_.param_names_tex = char(M_.param_names_tex, 'zeta\_l'); M_.param_names_long = char(M_.param_names_long, 'zeta_l'); M_.param_names = char(M_.param_names, 'zeta_c'); M_.param_names_tex = char(M_.param_names_tex, 'zeta\_c'); M_.param_names_long = char(M_.param_names_long, 'zeta_c'); M_.param_names = char(M_.param_names, 'zeta_p'); M_.param_names_tex = char(M_.param_names_tex, 'zeta\_p'); M_.param_names_long = char(M_.param_names_long, 'zeta_p'); M_.param_names = char(M_.param_names, 'tau'); M_.param_names_tex = char(M_.param_names_tex, 'tau'); M_.param_names_long = char(M_.param_names_long, 'tau'); M_.param_names = char(M_.param_names, 'rho_r'); M_.param_names_tex = char(M_.param_names_tex, 'rho\_r'); M_.param_names_long = char(M_.param_names_long, 'rho_r'); M_.param_names = char(M_.param_names, 'psi_pi'); M_.param_names_tex = char(M_.param_names_tex, 'psi\_pi'); M_.param_names_long = char(M_.param_names_long, 'psi_pi'); M_.param_names = char(M_.param_names, 'psi_y'); M_.param_names_tex = char(M_.param_names_tex, 'psi\_y'); M_.param_names_long = char(M_.param_names_long, 'psi_y'); M_.param_names = char(M_.param_names, 'psi_ex'); M_.param_names_tex = char(M_.param_names_tex, 'psi\_ex'); M_.param_names_long = char(M_.param_names_long, 'psi_ex'); M_.param_names = char(M_.param_names, 'rho_b'); M_.param_names_tex = char(M_.param_names_tex, 'rho\_b'); M_.param_names_long = char(M_.param_names_long, 'rho_b'); M_.param_names = char(M_.param_names, 'rho_istar'); M_.param_names_tex = char(M_.param_names_tex, 'rho\_istar'); M_.param_names_long = char(M_.param_names_long, 'rho_istar'); M_.param_names = char(M_.param_names, 'rho_p'); M_.param_names_tex = char(M_.param_names_tex, 'rho\_p'); M_.param_names_long = char(M_.param_names_long, 'rho_p'); M_.param_names = char(M_.param_names, 'rho_y'); M_.param_names_tex = char(M_.param_names_tex, 'rho\_y'); M_.param_names_long = char(M_.param_names_long, 'rho_y'); M_.param_names = char(M_.param_names, 'rho_x'); M_.param_names_tex = char(M_.param_names_tex, 'rho\_x'); M_.param_names_long = char(M_.param_names_long, 'rho_x'); M_.param_names = char(M_.param_names, 'rho_t'); M_.param_names_tex = char(M_.param_names_tex, 'rho\_t'); M_.param_names_long = char(M_.param_names_long, 'rho_t'); M_.param_names = char(M_.param_names, 'rho_g'); M_.param_names_tex = char(M_.param_names_tex, 'rho\_g'); M_.param_names_long = char(M_.param_names_long, 'rho_g'); M_.param_names = char(M_.param_names, 'mu_istar'); M_.param_names_tex = char(M_.param_names_tex, 'mu\_istar'); M_.param_names_long = char(M_.param_names_long, 'mu_istar'); M_.param_names = char(M_.param_names, 'mu_p'); M_.param_names_tex = char(M_.param_names_tex, 'mu\_p'); M_.param_names_long = char(M_.param_names_long, 'mu_p'); M_.param_names = char(M_.param_names, 'mu_t'); M_.param_names_tex = char(M_.param_names_tex, 'mu\_t'); M_.param_names_long = char(M_.param_names_long, 'mu_t'); M_.param_names = char(M_.param_names, 'theta1'); M_.param_names_tex = char(M_.param_names_tex, 'theta1'); M_.param_names_long = char(M_.param_names_long, 'theta1'); M_.param_names = char(M_.param_names, 'C_Yss'); M_.param_names_tex = char(M_.param_names_tex, 'C\_Yss'); M_.param_names_long = char(M_.param_names_long, 'C_Yss'); M_.param_names = char(M_.param_names, 'Ce_Yss'); M_.param_names_tex = char(M_.param_names_tex, 'Ce\_Yss'); M_.param_names_long = char(M_.param_names_long, 'Ce_Yss'); M_.param_names = char(M_.param_names, 'I_Yss'); M_.param_names_tex = char(M_.param_names_tex, 'I\_Yss'); M_.param_names_long = char(M_.param_names_long, 'I_Yss'); M_.param_names = char(M_.param_names, 'X_Yss'); M_.param_names_tex = char(M_.param_names_tex, 'X\_Yss'); M_.param_names_long = char(M_.param_names_long, 'X_Yss'); M_.param_names = char(M_.param_names, 'Ym_Yss'); M_.param_names_tex = char(M_.param_names_tex, 'Ym\_Yss'); M_.param_names_long = char(M_.param_names_long, 'Ym_Yss'); M_.param_names = char(M_.param_names, 'A_Yss'); M_.param_names_tex = char(M_.param_names_tex, 'A\_Yss'); M_.param_names_long = char(M_.param_names_long, 'A_Yss'); M_.param_names = char(M_.param_names, 'D_Yss'); M_.param_names_tex = char(M_.param_names_tex, 'D\_Yss'); M_.param_names_long = char(M_.param_names_long, 'D_Yss'); M_.param_names = char(M_.param_names, 'De_Yss'); M_.param_names_tex = char(M_.param_names_tex, 'De\_Yss'); M_.param_names_long = char(M_.param_names_long, 'De_Yss'); M_.param_names = char(M_.param_names, 'Ce_NWss'); M_.param_names_tex = char(M_.param_names_tex, 'Ce\_NWss'); M_.param_names_long = char(M_.param_names_long, 'Ce_NWss'); M_.param_names = char(M_.param_names, 'We_NWss'); M_.param_names_tex = char(M_.param_names_tex, 'We\_NWss'); M_.param_names_long = char(M_.param_names_long, 'We_NWss'); M_.param_names = char(M_.param_names, 'Pn_Pss'); M_.param_names_tex = char(M_.param_names_tex, 'Pn\_Pss'); M_.param_names_long = char(M_.param_names_long, 'Pn_Pss'); M_.param_names = char(M_.param_names, 'Pm_Pss'); M_.param_names_tex = char(M_.param_names_tex, 'Pm\_Pss'); M_.param_names_long = char(M_.param_names_long, 'Pm_Pss'); M_.param_names = char(M_.param_names, 'Pstar_Pss'); M_.param_names_tex = char(M_.param_names_tex, 'Pstar\_Pss'); M_.param_names_long = char(M_.param_names_long, 'Pstar_Pss'); M_.param_names = char(M_.param_names, 'Pm_Pstar'); M_.param_names_tex = char(M_.param_names_tex, 'Pm\_Pstar'); M_.param_names_long = char(M_.param_names_long, 'Pm_Pstar'); M_.param_names = char(M_.param_names, 'awbar'); M_.param_names_tex = char(M_.param_names_tex, 'awbar'); M_.param_names_long = char(M_.param_names_long, 'awbar'); M_.param_names = char(M_.param_names, 'bwbar'); M_.param_names_tex = char(M_.param_names_tex, 'bwbar'); M_.param_names_long = char(M_.param_names_long, 'bwbar'); M_.param_names = char(M_.param_names, 'cwbar'); M_.param_names_tex = char(M_.param_names_tex, 'cwbar'); M_.param_names_long = char(M_.param_names_long, 'cwbar'); M_.param_names = char(M_.param_names, 'D_ss'); M_.param_names_tex = char(M_.param_names_tex, 'D\_ss'); M_.param_names_long = char(M_.param_names_long, 'D_ss'); M_.param_names = char(M_.param_names, 'Rk_ss'); M_.param_names_tex = char(M_.param_names_tex, 'Rk\_ss'); M_.param_names_long = char(M_.param_names_long, 'Rk_ss'); M_.param_names = char(M_.param_names, 'Re_ss'); M_.param_names_tex = char(M_.param_names_tex, 'Re\_ss'); M_.param_names_long = char(M_.param_names_long, 'Re_ss'); M_.param_names = char(M_.param_names, 'istar_ss'); M_.param_names_tex = char(M_.param_names_tex, 'istar\_ss'); M_.param_names_long = char(M_.param_names_long, 'istar_ss'); M_.param_names = char(M_.param_names, 'S_ss'); M_.param_names_tex = char(M_.param_names_tex, 'S\_ss'); M_.param_names_long = char(M_.param_names_long, 'S_ss'); M_.param_names = char(M_.param_names, 'i_ss'); M_.param_names_tex = char(M_.param_names_tex, 'i\_ss'); M_.param_names_long = char(M_.param_names_long, 'i_ss'); M_.param_names = char(M_.param_names, 'pi_ss'); M_.param_names_tex = char(M_.param_names_tex, 'pi\_ss'); M_.param_names_long = char(M_.param_names_long, 'pi_ss'); M_.param_names = char(M_.param_names, 'ggtrend'); M_.param_names_tex = char(M_.param_names_tex, 'ggtrend'); M_.param_names_long = char(M_.param_names_long, 'ggtrend'); M_.param_names = char(M_.param_names, 'gtrend'); M_.param_names_tex = char(M_.param_names_tex, 'gtrend'); M_.param_names_long = char(M_.param_names_long, 'gtrend'); M_.param_names = char(M_.param_names, 'betabar'); M_.param_names_tex = char(M_.param_names_tex, 'betabar'); M_.param_names_long = char(M_.param_names_long, 'betabar'); M_.param_names = char(M_.param_names, 'deltabar'); M_.param_names_tex = char(M_.param_names_tex, 'deltabar'); M_.param_names_long = char(M_.param_names_long, 'deltabar'); M_.exo_det_nbr = 0; M_.exo_nbr = 7; M_.endo_nbr = 47; M_.param_nbr = 69; M_.orig_endo_nbr = 47; M_.aux_vars = []; options_.varobs = []; options_.varobs = 'dYobs'; options_.varobs = char(options_.varobs, 'dCobs'); options_.varobs = char(options_.varobs, 'dIobs'); options_.varobs = char(options_.varobs, 'dTBYobs'); options_.varobs = char(options_.varobs, 'PIobs'); options_.varobs = char(options_.varobs, 'iobs'); options_.varobs = char(options_.varobs, 'dSobs'); options_.varobs_id = [ 1 2 3 4 5 6 7 ]; M_.Sigma_e = zeros(7, 7); M_.Correlation_matrix = eye(7, 7); M_.H = 0; M_.Correlation_matrix_ME = 1; options_.linear = 1; options_.block=0; options_.bytecode=0; options_.use_dll=0; erase_compiled_function('mydata1_static'); erase_compiled_function('mydata1_dynamic'); M_.lead_lag_incidence = [ 0 28 0; 0 29 0; 0 30 0; 0 31 0; 0 32 0; 0 33 0; 0 34 0; 0 35 0; 0 36 0; 1 37 0; 2 38 75; 3 39 76; 4 40 0; 0 41 0; 5 42 0; 0 43 77; 6 44 78; 0 45 0; 0 46 0; 7 47 0; 8 48 79; 9 49 0; 10 50 80; 11 51 81; 0 52 0; 12 53 0; 13 54 0; 0 55 0; 0 56 0; 0 57 82; 14 58 0; 0 59 0; 15 60 0; 16 61 83; 0 62 0; 0 63 0; 17 64 84; 18 65 0; 19 66 0; 20 67 0; 21 68 0; 22 69 85; 23 70 0; 24 71 0; 25 72 0; 26 73 0; 27 74 0;]'; M_.nstatic = 18; M_.nfwrd = 2; M_.npred = 18; M_.nboth = 9; M_.nsfwrd = 11; M_.nspred = 27; M_.ndynamic = 29; M_.equations_tags = { }; M_.static_and_dynamic_models_differ = 0; M_.exo_names_orig_ord = [1:7]; M_.maximum_lag = 1; M_.maximum_lead = 1; M_.maximum_endo_lag = 1; M_.maximum_endo_lead = 1; oo_.steady_state = zeros(47, 1); M_.maximum_exo_lag = 0; M_.maximum_exo_lead = 0; oo_.exo_steady_state = zeros(7, 1); M_.params = NaN(69, 1); M_.NNZDerivatives = zeros(3, 1); M_.NNZDerivatives(1) = 215; M_.NNZDerivatives(2) = 0; M_.NNZDerivatives(3) = -1; M_.params( 1 ) = 0.350; alpha = M_.params( 1 ); M_.params( 2 ) = 0.985; beta = M_.params( 2 ); M_.params( 3 ) = 0.025; delta = M_.params( 3 ); M_.params( 5 ) = 0.492; theta = M_.params( 5 ); M_.params( 4 ) = 0.749; gamma = M_.params( 4 ); M_.params( 6 ) = 0.955; omega = M_.params( 6 ); M_.params( 7 ) = 0.985; nu = M_.params( 7 ); M_.params( 8 ) = 1.01; epsilon_w = M_.params( 8 ); M_.params( 9 ) = 1.5014; phi = M_.params( 9 ); M_.params( 10 ) = 3.0619; eta = M_.params( 10 ); M_.params( 11 ) = 0.2847; hab = M_.params( 11 ); M_.params( 12 ) = 1.5030; phi_c = M_.params( 12 ); M_.params( 13 ) = 0.0007; psi_D = M_.params( 13 ); M_.params( 14 ) = 0.3508; l_n = M_.params( 14 ); M_.params( 15 ) = 0.5844; xi_n = M_.params( 15 ); M_.params( 16 ) = 0.4562; l_m = M_.params( 16 ); M_.params( 17 ) = 0.5966; xi_m = M_.params( 17 ); M_.params( 18 ) = 0.4294; l_w = M_.params( 18 ); M_.params( 19 ) = 0.8280; xi_w = M_.params( 19 ); M_.params( 20 ) = 1.8826; chi = M_.params( 20 ); M_.params( 24 ) = 0.4374; zeta_c = M_.params( 24 ); M_.params( 22 ) = 0.4993; zeta_e = M_.params( 22 ); M_.params( 23 ) = 0.6915; zeta_l = M_.params( 23 ); M_.params( 25 ) = 3.7493; zeta_p = M_.params( 25 ); M_.params( 21 ) = 0.1019; zeta_u = M_.params( 21 ); M_.params( 26 ) = 0.4965; tau = M_.params( 26 ); M_.params( 27 ) = 0.8407; rho_r = M_.params( 27 ); M_.params( 28 ) = 1.1434; psi_pi = M_.params( 28 ); M_.params( 29 ) = 0.0002; psi_y = M_.params( 29 ); M_.params( 30 ) = 1.1389; psi_ex = M_.params( 30 ); M_.params( 31 ) = 0.7655; rho_b = M_.params( 31 ); M_.params( 32 ) = 0.6785; rho_istar = M_.params( 32 ); M_.params( 37 ) = 0.9034; rho_g = M_.params( 37 ); M_.params( 33 ) = 0.7965; rho_p = M_.params( 33 ); M_.params( 34 ) = 0.9998; rho_y = M_.params( 34 ); M_.params( 35 ) = 0.8877; rho_x = M_.params( 35 ); M_.params( 36 ) = 0.9800; rho_t = M_.params( 36 ); M_.params( 38 ) = 0.1206; mu_istar = M_.params( 38 ); M_.params( 39 ) = 0.6111; mu_p = M_.params( 39 ); M_.params( 40 ) = 0.9297; mu_t = M_.params( 40 ); M_.params( 41 ) = 0.4957; theta1 = M_.params( 41 ); M_.params( 42 ) = 0.463; C_Yss = M_.params( 42 ); M_.params( 43 ) = 0.035; Ce_Yss = M_.params( 43 ); M_.params( 44 ) = 0.284; I_Yss = M_.params( 44 ); M_.params( 45 ) = 0.991; X_Yss = M_.params( 45 ); M_.params( 46 ) = 0.861; Ym_Yss = M_.params( 46 ); M_.params( 47 ) = 0.798; A_Yss = M_.params( 47 ); M_.params( 48 ) = 0.120; D_Yss = M_.params( 48 ); M_.params( 49 ) = 0.126; De_Yss = M_.params( 49 ); M_.params( 50 ) = 0.016; Ce_NWss = M_.params( 50 ); M_.params( 51 ) = 0.008; We_NWss = M_.params( 51 ); M_.params( 52 ) = 0.911; Pn_Pss = M_.params( 52 ); M_.params( 53 ) = 0.995; Pm_Pss = M_.params( 53 ); M_.params( 54 ) = 1.056; Pstar_Pss = M_.params( 54 ); M_.params( 55 ) = 0.691; Pm_Pstar = M_.params( 55 ); M_.params( 56 ) = 0.434; awbar = M_.params( 56 ); M_.params( 57 ) = 0.756; bwbar = M_.params( 57 ); M_.params( 58 ) = 0.007; cwbar = M_.params( 58 ); M_.params( 59 ) = 0.027; D_ss = M_.params( 59 ); M_.params( 60 ) = 0.041; Rk_ss = M_.params( 60 ); M_.params( 61 ) = 1.034; Re_ss = M_.params( 61 ); M_.params( 62 ) = 0.032; istar_ss = M_.params( 62 ); M_.params( 63 ) = 1; S_ss = M_.params( 63 ); M_.params( 64 ) = 0.043; i_ss = M_.params( 64 ); M_.params( 65 ) = 2.960; pi_ss = M_.params( 65 ); M_.params( 66 ) = 1.0072; ggtrend = M_.params( 66 ); M_.params( 67 ) = (M_.params(66)-1)*100; gtrend = M_.params( 67 ); M_.params( 68 ) = M_.params(2)*M_.params(66)^(-M_.params(9)); betabar = M_.params( 68 ); M_.params( 69 ) = 1-(1-M_.params(3))/M_.params(66); deltabar = M_.params( 69 ); % % SHOCKS instructions % make_ex_; M_.exo_det_length = 0; M_.Sigma_e(1, 1) = (0.0391)^2; M_.Sigma_e(2, 2) = (0.7636)^2; M_.Sigma_e(3, 3) = (0.8482)^2; M_.Sigma_e(4, 4) = (4.9305)^2; M_.Sigma_e(5, 5) = (1.9104)^2; M_.Sigma_e(6, 6) = (1.9063)^2; M_.Sigma_e(7, 7) = (0.1591)^2; M_.sigma_e_is_diagonal = 1; resid(1) steady; oo_.dr.eigval = check(M_,options_,oo_); global estim_params_ estim_params_.var_exo = []; estim_params_.var_endo = []; estim_params_.corrx = []; estim_params_.corrn = []; estim_params_.param_vals = []; estim_params_.var_exo = [estim_params_.var_exo; 1, 0.0391, 0.01, 4, 4, 0.1, 2, NaN, NaN, NaN ]; estim_params_.var_exo = [estim_params_.var_exo; 2, 0.7636, 0.01, 4, 4, 0.65, 2, NaN, NaN, NaN ]; estim_params_.var_exo = [estim_params_.var_exo; 3, 0.8482, 0.01, 4, 4, 0.9, 2, NaN, NaN, NaN ]; estim_params_.var_exo = [estim_params_.var_exo; 4, 4.9305, 0.01, 8, 4, 3.5, 2, NaN, NaN, NaN ]; estim_params_.var_exo = [estim_params_.var_exo; 5, 1.9104, 0.01, 4, 4, 1.35, 2, NaN, NaN, NaN ]; estim_params_.var_exo = [estim_params_.var_exo; 6, 1.9063, 0.01, 4, 4, 1.25, 2, NaN, NaN, NaN ]; estim_params_.var_exo = [estim_params_.var_exo; 7, 1.3268, 0.01, 4, 4, 0.15, 2, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 31, 0.7655, 0.1, 0.9999, 1, 0.85, 0.10, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 32, 0.6785, 0.1, 0.9999, 1, 0.85, 0.10, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 33, 0.7965, 0.1, 0.9999, 1, 0.85, 0.10, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 34, 0.9998, 0.1, 0.9999, 1, 0.85, 0.10, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 35, 0.8877, 0.1, 0.9999, 1, 0.85, 0.10, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 36, 0.9800, 0.1, 0.9999, 1, 0.85, 0.10, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 37, 0.9034, 0.1, 0.9999, 1, 0.85, 0.10, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 38, 0.1206, 0.01, 0.9999, 1, 0.5, 0.2, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 39, 0.6111, 0.01, 0.9999, 1, 0.7, 0.2, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 40, 0.9297, 0.01, 0.9999, 1, 0.7, 0.2, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 9, 2.3014, 0.25, 4, 3, 2.1, 0.375, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 10, 3.0619, 0.5, 5, 3, 2.5, 0.75, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 11, 0.2847, 0.1, 0.95, 1, 0.57, 0.1, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 12, 1.5030, 1.0, 3, 3, 1.50, 0.125, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 13, 0.0007, 0.00001, 2, 3, 0.001, 0.05, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 14, 0.3508, 0.1, 0.99, 1, 0.40, 0.15, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 15, 0.5844, 0.1, 0.99, 1, 0.57, 0.05, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 16, 0.4562, 0.1, 0.99, 1, 0.40, 0.15, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 17, 0.5966, 0.1, 0.99, 1, 0.57, 0.05, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 18, 0.4294, 0.1, 0.99, 1, 0.50, 0.15, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 19, 0.8280, 0.1, 0.99, 1, 0.75, 0.05, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 20, 1.8826, 0.1, 5, 3, 1.48, 1.5, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 24, 0.4374, 0.01, 0.99, 1, 0.45, 0.15, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 22, 0.4993, 0.01, 0.99, 1, 0.50, 0.15, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 23, 0.6915, 0.01, 0.99, 1, 0.65, 0.15, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 25, 3.7493, 0.001, 7.5, 3, 3.75, 0.5, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 21, 0.1019, 0.01, 0.99, 1, 0.5, 0.15, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 26, 0.4965, 0.01, 0.9999, 1, 0.495, 0.15, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 27, 0.8407, 0.3, 0.9999, 1, 0.75, 0.15, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 28, 1.1434, 1.0, 3, 3, 1.5, 0.25, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 29, 0.0002, 0.000001, 0.25, 3, 0.125, 0.05, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 30, 1.1389, 0.01, 2.0, 3, 1.112, 0.05, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 67, 0.720, 0.10, 1.10, 2, 0.784, 0.011, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 65, 2.960, 4.00, 10.0, 2, 8.712, 0.056, NaN, NaN, NaN ]; estim_params_.param_vals = [estim_params_.param_vals; 41, 0.430, 0.001, 10.0, 2, 0.430, 0.1, NaN, NaN, NaN ]; options_.bayesian_irf = 1; options_.lik_init = 1; options_.mh_jscale = 0.2; options_.mh_nblck = 2; options_.mh_replic = 2000; options_.mode_check.status = 1; options_.mode_compute = 5; options_.datafile = 'mysdata1'; options_.optim_opt = '''MaxIter'',1000000,''MaxFunEvals'',500000'; options_.order = 1; var_list_=[]; var_list_ = 'dYobs'; var_list_ = char(var_list_, 'dCobs'); var_list_ = char(var_list_, 'dIobs'); var_list_ = char(var_list_, 'dTBYobs'); var_list_ = char(var_list_, 'PIobs'); var_list_ = char(var_list_, 'iobs'); var_list_ = char(var_list_, 'dSobs'); dynare_estimation(var_list_); options_.ar = 10; options_.irf = 20; options_.conditional_variance_decomposition = [1;2;4;8;16;32;64;100;]; var_list_=[]; var_list_ = 'dYobs'; var_list_ = char(var_list_, 'dCobs'); var_list_ = char(var_list_, 'dIobs'); var_list_ = char(var_list_, 'PIobs'); var_list_ = char(var_list_, 'iobs'); var_list_ = char(var_list_, 'dSobs'); var_list_ = char(var_list_, 'dTBYobs'); info = stoch_simul(var_list_); save('mydata1_results.mat', 'oo_', 'M_', 'options_'); if exist('estim_params_', 'var') == 1 save('mydata1_results.mat', 'estim_params_', '-append'); end if exist('bayestopt_', 'var') == 1 save('mydata1_results.mat', 'bayestopt_', '-append'); end if exist('dataset_', 'var') == 1 save('mydata1_results.mat', 'dataset_', '-append'); end if exist('estimation_info', 'var') == 1 save('mydata1_results.mat', 'estimation_info', '-append'); end disp(['Total computing time : ' dynsec2hms(toc) ]); if ~isempty(lastwarn) disp('Note: warning(s) encountered in MATLAB/Octave code') end diary off