var rho i_d p i_f e zeta_c c w zeta_l l c_h c_f p_h p_f epsilon_d i_s d alpha_s alpha_d i i_t i_l epsilon_l i_b l_d z pi_b ad_i_d ad_i_l ad_z y theta k mc r r_k q f omega lev pi pi_ff rer n pi_f prob v c_e c_e_h c_e_f zeta_i inv ad_inv inv_h inv_f epsilon_var inf_h inf_f g_h g_f g q_x p_x zeta; varexo epsilon_zeta_c epsilon_zeta_l e_d epsilon_alpha_s epsilon_alpha_d epsilon_i_t m e_l epsilon_theta epsilon_omega epsilon_pi_f epsilon_prob epsilon_i e_var e_i epsilon_g epsilon_z; parameters beta b psi alpha etha rho_zeta_c rho_zeta_l phi_i_d rho_epsilon_d rho_alpha_s rho_alpha_d rho_i_t phi_i_l phi_z kappa_d delta_b omega_b gamma rho_epsilon_l rho_theta delta rho_omega omega_barra sigma rho_pi_f rho_prob kappa rho_zeta_i nu_d nu_f rho_epsilon_var rho_i rho_pi rho_y rho_g rho_z choque; beta = 0.998158027126013; b = 0.4986; psi = 1.5; alpha = 0.8; etha = 1.5; rho_zeta_c = 0.5; rho_zeta_l = 0.5; phi_i_d = 20; rho_epsilon_d = 0.5; rho_alpha_s = 0.5; rho_alpha_d = 0.5; rho_i_t = 0.5; phi_i_l = 20; phi_z = 10; kappa_d = 0.1; delta_b = 0.035; omega_b = 0.0422616449087736; gamma = 0.33; rho_epsilon_l = 0.5; rho_theta = 0.5; delta = 0.025; rho_omega = 0.5; omega_barra = 0.38859981897582; sigma = -0.0566; rho_pi_f = 0.5; rho_prob = 0.5; kappa = 23.008; rho_zeta_i = 0.5; nu_d = 118.22; nu_f = 100.043; rho_epsilon_var = 0.5; rho_i = 0.8; rho_pi = 2; rho_y = 0.01; rho_g = 0.5; rho_z = 0.5; choque = 0.00007; model; // 1.HOUSEHOLDS rho = i_d*(p/p(+1)); // 1. Ecuación de Euler i_f*(e(+1)/e) = i_d; // 2. Eficiencia de Mercados rho = (zeta_c/zeta_c(+1))*(1/beta)*(((1/(c-b*c(-1)))-(beta*b)/(c(+1)-b*c))/((1/(c(+1)-b*c))-(beta*b)/(c(+2)-b*c(+1)))); // 3. Definición de Variable Auxiliar w/p = (zeta_l/zeta_c)*(l^(psi-1))/((1/(c-b*c(-1)))-(beta*b)/(c(+1)-b*c)); // 4. Eficiencia del Trabajo c_h = alpha*((p_h/p)^(-etha))*c; // 5. Consumo Doméstico c_f = (1-alpha)*((p_f/p)^(-etha))*c; // 6. Consumo Foráneo p = (alpha*(p_h^(1-etha))+(1-alpha)*(p_f^(1-etha)))^(1/(1-etha)); // 7. Índice de Precios // 2. BANKS ((1+epsilon_d)/epsilon_d)*i_d = i_s-(phi_i_d/epsilon_d)*((i_d/i_d(-1))-1)*(i_d/i_d(-1))+((beta*phi_i_d)/epsilon_d)*((i_d(+1)/i_d)-1)*(i_d(+1)/i_d)*(d(+1)/d); // 1. Tasa de Interés para Depósitos i_s = (1-alpha_s-alpha_d)*i+alpha_d*i_t; // 2. Tasa de Interés para Depósitos Wholesale i_l = (epsilon_l/(epsilon_l-1))*i_b-(phi_i_l/(epsilon_l-1))*((i_l/i_l(-1))-1)*(i_l/i_l(-1))+(beta*phi_i_l/(epsilon_l-1))*((i_l(+1)/i_l)-1)*(i_l(+1)/i_l)*(l_d(+1)/l_d); // 3. Tasa de Interés para Préstamos i_b = i-phi_z*((z/l_d)-kappa_d)*((z/l_d)^2); // 4. Tasa de Interés para Préstamos Wholesale z = (1-delta_b)*z(-1)+omega_b*pi_b(-1)-m; // 5. Capital Bancario pi_b = i_l*l_d-i_d*d-ad_i_d-ad_i_l-ad_z; // 6. Ganancias Bancarias ad_i_d = (phi_i_d/2)*(((i_d/i_d(-1))-1)^2)*d; // 7. Definición de Variable Auxiliar ad_i_l = (phi_i_l/2)*(((i_l/i_l(-1))-1)^2)*l_d; // 8. Definición de Variable Auxiliar ad_z = (phi_z/2)*(((z/l_d)-kappa_d)^2)*z; // 9. Definición de Variable Auxiliar // 3. ENTERPRENEURS y = theta*(l^(1-gamma))*(k(-1)^gamma); // 1. Producción w = (1-gamma)*(y/l)*mc; // 2. Productividad Marginal del Trabajo r = gamma*(y/k(-1))*(mc/p); // 3. Productividad Marginal del Capital r_k = (r(+1)+(1-delta)*q)/q(-1); // 4. Rendimiento Bruto del Capital f = (1+omega)*(omega_barra*(lev^sigma)*(i_l/pi(+1))+(1-omega_barra)*(lev^sigma)*(i_f/pi_ff)*(rer(+1)/rer)); // 5. Costo de Préstamos r_k = f; // 6. Condición de Eficiencia lev = n/(q*k); // 7. Apalancamiento rer = (e*pi_f)/p; // 8. Tipo de Cambio Real pi = p/p(-1); // 9. Inflación Bruta n = prob*v; // 10. Valor Esperado de la Riqueza de Emprendedores v = r_k(-1)*q(-1)*k(-1)-(1+omega(-1))*(omega_barra*(q(-1)*k(-1)-n(-1))*(lev(-1)^sigma)*(i_l(-1)/pi)+(1-omega_barra)*(q(-1)*k(-1)-n(-1))*(lev(-1)^sigma)*(i_f(-1)/pi_ff(-1))*(rer/rer(-1))); // 11. Riqueza de Emprendedores c_e = (1-prob)*v; // 12. Consumo Emprendedores c_e_h = alpha*((p_h/p)^(-etha))*c_e; // 13. Consumo Doméstico Emprendedores c_e_f = (1-alpha)*((p_f/p)^(-etha))*c_e; // 14. Consumo Foráneo Emprendedores mc = (w*l+r*k(-1))/y; // 15. Costo Marginal Emprendedores pi_ff = pi_f/pi_f(-1); // 16. Inflación Foránea // 4. CAPITAL PRODUCERS q = ((zeta_i-kappa*((inv/k)-delta))^-1); // 1. Costo del Capital k = (1-delta)*k(-1)+zeta_i*inv-ad_inv; // 2. Evolución del Capital ad_inv = (kappa/2)*(((inv/k)-delta)^2)*k; // 3. Definición de Variable Auxiliar inv_h = alpha*((p_h/p)^(-etha))*inv; // 4. Inversión Doméstica inv_f = (1-alpha)*((p_f/p)^(-etha))*inv; // 5. Inversión Foránea // 5. RETAILERS p_h = (epsilon_var/(epsilon_var-1))*mc-(nu_d/(epsilon_var-1))*p*inf_h*(inf_h-1)+(nu_d/(epsilon_var-1))*p*rho*(y(+1)/y)*inf_h(+1)*(inf_h(+1)-1); // 1. Precio Retail Doméstico p_f = (epsilon_var/(epsilon_var-1))*mc-(nu_f/(epsilon_var-1))*p*inf_f*(inf_f-1)+(nu_f/(epsilon_var-1))*p*rho*(y(+1)/y)*inf_f(+1)*(inf_f(+1)-1); // 2. Precio Retail Foráneo inf_h = p_h/p_h(-1); // 3. Inflación Precio Retail Doméstico inf_f = p_f/p_f(-1); // 4. Inflación Precio Retail Foráneo // 6. MONETARY POLICY i = rho_i*i(-1)+(1-rho_i)*1+rho_y*(y-1)+rho_pi*(pi-1)+e_i; // 1. Tasa de Referencia del Banco Central // 7. GOVERNMENT g_h = alpha*((p_h/p)^(-etha))*g; // 1. Gasto Doméstico g_f = (1-alpha)*((p_f/p)^(-etha))*g; // 2. Gasto Foráneo // 8. MARKET CLEARING y = c_h+c_e_h+inv_h+g_h+q_x; // 1. Vaciado de Mercado q_x = (p_x/pi_f)^(-zeta); // 2. Definición de Variable Auxiliar p_x = p_h/e; // 3. Ley de un Solo Precio d = l_d; // 4. Vaciado del Mercado Bancario l_d = omega_barra*(q*k-n); // 5. Préstamos Locales // 9. EXOGENOUS PROCESSES zeta_c = (1-rho_zeta_c)+rho_zeta_c*zeta_c(-1)+epsilon_zeta_c; // 1. Choque a Preferencias de Consumo zeta_l = (1-rho_zeta_l)*(w*(zeta_c/(l^(psi-1)))*((1/(c*(1-b)))-(beta*b)/(c*(1-b))))+rho_zeta_l*zeta_l(-1)+epsilon_zeta_l; // 2. Choque a Preferencias de Ocio epsilon_d = (1-rho_epsilon_d)*31.12+rho_epsilon_d*epsilon_d(-1)+e_d; // 3. Choque a Elasticidad de Depósitos alpha_s = (1-rho_alpha_s)+rho_alpha_s*alpha_s(-1)+epsilon_alpha_s; // 4. Choque a CRR alpha_d = (1-rho_alpha_d)*0.584875199219751+rho_alpha_d*alpha_d(-1)+epsilon_alpha_d; // 5. Choque a SLR i_t = (1-rho_i_t)*1.05819933333333+rho_i_t*i_t(-1)+epsilon_i_t; // 6. Choque a Cetes epsilon_l = (1-rho_epsilon_l)*(i_l/(i_l-i_b))+rho_epsilon_l*epsilon_l(-1)+e_l; // 7. Choque a Elasticidad de Préstamos theta = (1-rho_theta)+rho_theta*theta(-1)+epsilon_theta; // 8. Choque a la Productividad omega = (1-rho_omega)*((r_k/f)-1)+rho_omega*omega(-1)+epsilon_omega; // 9. Choque a Costo de Fondos pi_f = (1-rho_pi_f)+rho_pi_f*pi_f(-1)+epsilon_pi_f; // 10. Índice Foráneo de Precios prob = (1-rho_prob)*(n/v)+rho_prob*prob(-1)+epsilon_prob; // 11. Población de Emprendedores zeta_i = (1-rho_zeta_i)+rho_zeta_i*zeta_i(-1)+epsilon_i; // 12. Choque a la Productividad de la Inversión epsilon_var = (1-rho_epsilon_var)*(1.09/0.09)+rho_epsilon_var*epsilon_var(-1)+e_var; // 13. Choque a Preferencia por Variedad g = (1-rho_g)*0.345367122373301+rho_g*g(-1)+epsilon_g; // 14. Choque a Gasto del Gobierno zeta = (1-rho_z)*4.89+rho_z*zeta(-1)+epsilon_z; // 15. Elasticidad Precio de las Exportaciones end; initval; rho = 1.00184537199915; i_d = 1.00184537199915; p = 1; i_f = 1.00184537199915; e = 0.5; zeta_c = 1; w = 1.05863389054603; l = 0.580634064875372; p_h = 1; p_f = 1; alpha_s = 0; alpha_d = 0.584875199219751; i = 1; i_t = 1.05819933333333; i_s = (1-alpha_s-alpha_d)*i+alpha_d*i_t; epsilon_d = i_d/(i_s-i_d); d = 0.585898038987953; i_l = 1.08466278262184; i_b = 1; epsilon_l = i_l/(i_l-i_b); l_d = 0.585898038987953; z = 0.0585898038987953; pi_b = 0.0485225584778911; ad_i_d = 0; ad_i_l = 0; ad_z = 0; y = 1; theta = 1; k = 3.01543135317009; mc = 0.91743119266055; r = 0.1004009901468; r_k = 1.0754009901468; q = 1; f = 1.07540160896276; omega = (r_k/f)-1; lev = 0.5; pi = 1; pi_ff = 1; rer = 0.5; n = 1.50771567658504; pi_f = 1; v = 1.62139799846089; prob = n/v; c_e = 0.113682321875847; c_e_h = 0.0909458575006777; c_e_f = 0.0227364643751694; zeta_i = 1; inv = 0.0753857838292522; ad_inv = 0; inv_h = 0.0603086270634017; inv_f = 0.0150771567658504; epsilon_var = 1.09/0.09; inf_h = 1; inf_f = 1; g_h = 0.276293697898641; g_f = 0.0690734244746602; g = 0.345367122373301; q_x = 0.0337258823907634; p_x = 2; zeta = 4.89; epsilon_zeta_c = 0; epsilon_zeta_l = 0; e_d = 0; epsilon_alpha_s = -0.5; epsilon_alpha_d = 0; epsilon_i_t = 0; m = 0; e_l = 0; epsilon_theta = 0; epsilon_omega = 0; epsilon_pi_f = 0; epsilon_prob = 0; epsilon_i = 0; e_var = 0; e_i = 0; epsilon_g = 0; epsilon_z = 0; c_h = 1-c_e_h-inv_h-g_h-q_x; c = c_h/alpha; c_f = (1-alpha)*c; zeta_l = w*(zeta_c/(l^(psi-1)))*((1/(c*(1-b)))-(beta*b)/(c*(1-b))); end; steady; check; shocks; var m = choque^2; end; stoch_simul;