//************************************************************************* // Mon Modele DSGE //************************************************************************* //------------------------------------------------------------------------- // Partie 1: Declaration des variables endogènes //------------------------------------------------------------------------- var Ctild_star % 01- La consommation privée à l'echelle de l'union pi_star % 02- Le taux d'inflation de l'union pi_s % 03- L'inflation à la production le bloc de pays s pi_B % 04- L'inflation à la production le bloc de pays B ytild_star % 05- La production à l'échelle de l'union Ctild_s % 06- La consommation privée dans le bloc de pays s Ctild_B % 07 - La consommation privée dans le bloc de pays B tt_tild_s % 08 - Ecart des termes de l'échange du bloc de pays s par rapport à l'union tt_tild_B % 09- Ecart des termes de l'échange du bloc de pays B par rapport à l'union ytild_s % 10- La production dans le bloc de pays s ytild_B % 11- La production dans le bloc de pays B dchapo_g_s % 12- La valeur réelle de la dette publique par habitant du bloc de pays s dchapo_g_B % 13- La valeur réelle de la dette publique par habitant du bloc de pays B rtild_star % 14- Le taux d'intérêt de court terme (instrument de la politique monétaire) gtild_s % 15- Les dépenses publiques du bloc de pays s gtild_B % 16- Les dépenses publiques du bloc pays B gtild_star % 17- Les dépenses publiques à l'échelle de l'union tautild_s % 18- L'impôt sur le revenu du bloc de pays s tautild_B % 19- L'impôt sur le revenu du bloc de pays B tt_tild_i % 20- Ecart des termes de l'échange du pays i par rapport à l'union Ctild_i % 21- La consommation privée du pays i tt_bar_s % 22- Le log-linéaire du premier ordre des termes de l'échange du petit pays s par rapport au régime permanent tt_bar_B % 23- Le log-linéaire du premier ordre des termes de l'échange du grand pays B par rapport au régime permanent Y_s % 24- La production en régime permanent du bloc de pays s D_g_s % 25- La valeur réelle de la dette en régime permanent du bloc de pays s (exprimée e prix à la consommation), en termes par habitant Y_B % 26- La production en régime permanent du bloc de pays B D_g_B % 27- La valeur réelle de la dette en régime permanent du bloc de pays B (exprimée e prix à la consommation), en termes par habitant r_bar_star % 28- Le log-linéaire du premier ordre du taux d'intérêt de court terme par rapport au régime permanent r_star % 29- taux d'intérêt réel en régime permanent (4% sur une base annuelle) C_bar_s % 30- La consommation privée en régime permanent dans le bloc de pays s C_bar_B % 31- La consommation privée en régime permanent dans le bloc de pays B C_bar_star % 32- La consommation privée en régime permanent à l'echelle de l'union g_bar_s % 33- Le log-linéaire du premier ordre des dépenses publiques du bloc de petit pays s par rapport au régime permanent g_bar_B % 34- Le log-linéaire du premier ordre des dépenses publiques de bloc du pays B par rapport au régime permanent g_bar_star % 35- Les dépenses publiques en régime permanent à l'echelle de l'union ybar_s % 36- Le log-linéaire du premier ordre de la production du bloc de petit pays s par rapport au régime permanent ybar_B % 37- Le log-linéaire du premier ordre de la production du bloc de pays B par rapport au régime permanent y_bar_star % 38- La production en régime permanent à l'echelle de l'union P_bar_star % 39- Indice des prix des biens importés en regime permanent à l'echelle de l'union P_bar_s % 40- Indice des prix des biens importés en regime permanent dans le bloc de pays s P_bar_B % 41- Indice des prix des biens importés en regime permanent dans le bloc de pays B tau_bar_s % 42- Le log-linéaire du premier ordre de l'impôt sur le revenu du bloc de petit pays s par rapport au régime permanent tau_bar_B % 43- Le log-linéaire du premier ordre de l'impôt sur le revenu du bloc de pays B par rapport au régime permanent tau_s % 44- L'impôt sur le revenu du bloc de pays s en régime permanent tau_B % 45- L'impôt sur le revenu du bloc de pays B en régime permanent muchapo_w_s % 46- revalorisation nette des salaires dans le bloc de petit pays s muchapo_w_B % 47- revalorisation nette des salaires dans le bloc de petit pays B ; //------------------------------------------------------------------------- // Partie 2: Declaration des variables exogènes //------------------------------------------------------------------------- varexo a_s % 01- innovation associée à la technologie dans le bloc de pays s a_B % 02- innovation associée à la technologie dans bloc de pays B a_star % 03- innovation associée à la technologie à l'echelle de l'union ; //------------------------------------------------------------------------- // Partie 3: Declaration des paramètres //------------------------------------------------------------------------- parameters beta rho sigma psi chi epsilon gamma phi tetha_s tetha_B rho_a epsilon_s epsilon_B alpha epsilon_star n ; // Calibration beta = 0.99; % facteur d'actualisation des préférences intertemporelles rho = -log(beta); % taux d'intérêt réel en régime permanent (4% sur une base annuelle) sigma = 0.4; % l'élasticité intertemporelle de substitution de la consommation privée psi = 0.4; % l'élasticité intertemporelle de substitution de la consommation plublique chi = 3; % l'inverse de l'élasticité d'offre de travail epsilon = 11; % l'élasticité de substitution entre les biens produits dans le même pays gamma = 4.5; % l'élasticité intertemporelle de substitution entre les biens nationaux et étrangers phi = 0.25; % part de la consommation publique en régime permanent dans la production tetha_s = 0.75; % La règle d'ajustement des prix à la Calvo pour les petits pays tetha_B = 0.75; % La règle d'ajustement des prix à la Calvo pour le grand pays rho_a = 0.85; % Paramètre d'autocorrélation du processus AR(1) de la technologie epsilon_s = 0.5; % Terme de l'erreur du processus AR(1) de la technologie dans le pays s epsilon_B = 0.5; % Terme de l'erreur du processus AR(1) de la technologie dans le pays B epsilon_star = 0.5; % Terme de l'erreur du processus AR(1) de la technologie a l'echelle de l'union alpha = 0.4; % Le dégré d'ouverture commerciale n = 0.5; % Dimension du bloc de pays s //------------------------------------------------------------------------- // Partie 5 : Declaration du modèle //------------------------------------------------------------------------- model; Ctild_star = Ctild_star(+1)-sigma*(rtild_star-pi_star(+1)); pi_s = beta*pi_s(+1)+(((1-tetha_s*beta)*(1-tetha_s))/(tetha_s*(1+epsilon*chi)))*((1/(1-phi)*(sigma*(1-alpha)+alpha*(gamma-(1-alpha)*(-gamma+sigma))))+chi)*ytild_s -(((1-tetha_s*beta)*(1-tetha_s))/(tetha_s*(1+epsilon*chi)))*(phi/(1-phi)*(sigma*(1-alpha)+alpha*(gamma-(1-alpha)*(-gamma+sigma))))*gtild_s +(((1-tetha_s*beta)*(1-tetha_s))/(tetha_s*(1+epsilon*chi)))*(1/(sigma*(1-phi)))-(1/(1-phi)*(sigma*(1-alpha))+ alpha*(gamma-(1-alpha)*(-gamma+sigma)))*ytild_star -(((1-tetha_s*beta)*(1-tetha_s))/(tetha_s*(1+epsilon*chi)))*phi*(1/(sigma*(1-phi)))-(1/(1-phi)*(sigma*(1-alpha))+ alpha*(gamma-(1-alpha)*(-gamma+sigma)))*gtild_star +((((1-tetha_s*beta)*(1-tetha_s))/(tetha_s*(1+epsilon*chi)))/(1-tau_s))*(exp(tautild_s)); pi_B = beta*pi_B(+1)+(((1-tetha_B*beta)*(1-tetha_B))/(tetha_B*(1+epsilon*chi)))*((1/(1-phi)*(sigma*(1-alpha)+alpha*(gamma-(1-alpha)*(-gamma+sigma))))+chi)*ytild_B -(((1-tetha_B*beta)*(1-tetha_B))/(tetha_B*(1+epsilon*chi)))*(phi/(1-phi)*(sigma*(1-alpha)+alpha*(gamma-(1-alpha)*(-gamma+sigma))))*gtild_B +(((1-tetha_B*beta)*(1-tetha_B))/(tetha_B*(1+epsilon*chi)))*(1/(sigma*(1-phi)))-(1/(1-phi)*(sigma*(1-alpha))+ alpha*(gamma-(1-alpha)*(-gamma+sigma)))*ytild_star -(((1-tetha_B*beta)*(1-tetha_B))/(tetha_B*(1+epsilon*chi)))*phi*(1/(sigma*(1-phi)))-(1/(1-phi)*(sigma*(1-alpha))+ alpha*(gamma-(1-alpha)*(-gamma+sigma)))*gtild_star +((((1-tetha_B*beta)*(1-tetha_B))/(tetha_B*(1+epsilon*chi)))/(1-tau_B))*(exp(tautild_B)); Ctild_s = Ctild_B +sigma*(1-alpha)*((tt_tild_s)-n*(tt_tild_B)); Ctild_s = Ctild_i + sigma*(1-alpha)*((tt_tild_s)-(tt_tild_i)); tt_tild_s = tt_tild_s(-1)+ pi_star-pi_s-(tt_bar_s - tt_bar_s(-1)); tt_tild_B = tt_tild_B(-1)+(pi_star-pi_B/n)-(tt_bar_B - tt_bar_B(-1)); ytild_s = (1-phi)*Ctild_star+(1-phi)*(sigma*(1-alpha)+(alpha*(gamma-(1-alpha)*(-gamma+sigma))))*tt_tild_s + phi*gtild_s; ytild_B = (1-phi)*Ctild_star +(1-phi)*(sigma*(1-alpha)+(alpha*(gamma-(1-alpha)*(-gamma+sigma))))*n*tt_tild_B + phi*gtild_B; dchapo_g_s = rtild_star +(1/beta)*(dchapo_g_s(-1)-pi_s +((Y_s/D_g_s)*(phi*gtild_s -tau_s*ytild_s-tautild_s)+alpha*tt_tild_s(-1)-(1/(1+r_star))*alpha*tt_tild_s ))+ r_bar_star+(1/beta)*((Y_s/D_g_s)*(phi*g_bar_s - tau_s*ybar_s - tau_bar_s)+alpha*tt_bar_s(-1)-(1/(1+r_star))*alpha*tt_bar_s); dchapo_g_B = rtild_star+(1/beta)*(dchapo_g_B(-1)-pi_B +((Y_B/D_g_B)*(phi*gtild_B -tau_B*ytild_B-tautild_B)+alpha*n*tt_tild_B(-1)-(1/(1+r_star))*alpha*n*tt_tild_B)) + r_bar_star +(1/beta)*((Y_B/D_g_B)*(phi*g_bar_B - tau_B*ybar_B - tau_bar_B) +alpha*n*tt_bar_B(-1)-(1/(1+r_star))*alpha*n*tt_bar_B); a_s = rho_a*a_s(-1)+ epsilon_s; a_B = rho_a*a_B(-1)+ epsilon_B; a_star = rho_a*a_star(-1)+ epsilon_star; end; //--------------------------------------------------------------------- // Partie 6: Des valeurs initiales et l'etat stationnaire //--------------------------------------------------------------------- steady_state_model; c_bar_s = ((1+alpha)*(1+chi)*sigma/(1+chi*(phi*psi+(1-phi)*(gamma+(sigma-gamma)*(1-2*alpha+alpha*alpha)))))*a_s +(alpha*(1+chi)*sigma*(1+chi*(phi*psi+(1-phi)*(gamma-(sigma-gamma)*(1-alpha))))/(1+chi*(phi*psi+(1-phi)*sigma))*(1+chi*(phi*psi+(1-phi)*(gamma+(sigma-gamma)*(1-2*alpha+alpha*alpha)))))*a_star; C_bar_B = ((1+alpha)*(1+chi)*sigma/(1+chi*(phi*psi+(1-phi)*(gamma+(sigma-gamma)*(1-2*alpha+alpha*alpha)))))*a_B +(alpha*(1+chi)*sigma*(1+chi*(phi*psi+(1-phi)*(gamma-(sigma-gamma)*(1-alpha))))/(1+chi*(phi*psi+(1-phi)*sigma))*(1+chi*(phi*psi+(1-phi)*(gamma+(sigma-gamma)*(1-2*alpha+alpha*alpha)))))*a_star; C_bar_star =((1+chi)*sigma/(1+chi*(phi*psi+(1-phi)*sigma)))*a_star; g_bar_s =((1+chi)*psi/(1+chi*(phi*psi+(1-phi)*(gamma+(sigma-gamma)*(1-2*alpha+alpha*alpha)))))* a_s -((1+chi)*psi*(1-phi)*chi*(sigma-gamma)*(2*alpha-alpha*alpha)/(1+chi*(phi*psi+(1-phi)*sigma))*(1+chi*(phi*psi+(1-phi)*(gamma+(sigma-gamma)*(1-2*alpha+alpha*alpha)))))*a_star; g_bar_B =((1+chi)*psi/(1+chi*(phi*psi+(1-phi)*(gamma+(sigma-gamma)*(1-2*alpha+alpha*alpha)))))* a_B -((1+chi)*psi*(1-phi)*chi*(sigma-gamma)*(2*alpha-alpha*alpha)/(1+chi*(phi*psi+(1-phi*sigma))*(1+chi*(phi*chi+(1-phi)*(gamma+(sigma-gamma)*(1-2*alpha+alpha*alpha))))))*a_star; g_bar_star =((1+chi)*psi/(1+chi*(phi*psi+(1-phi)*sigma)))*a_star; y_bar_s =((1+chi)*(phi*psi+(1-phi)*(gamma+(sigma-gamma)*(1-2*alpha+alpha*alpha)))/(1+chi*(phi*psi+(1-phi)*(gamma+(sigma-gamma)*(1-2*alpha+alpha*alpha)))))*a_s +((1+chi)*(1-phi)*(sigma-gamma)*(2*alpha-alpha*alpha)/(1+chi*(phi*psi+((1-phi)*sigma))*(1+chi*(phi*chi+(1-phi)*(gamma+(sigma-gamma)*(1-2*alpha+alpha*alpha))))))*a_star; y_bar_B =((1+chi)*(phi*psi+(1-phi)*(gamma+(sigma-gamma)*(1-2*alpha+alpha*alpha)))/(1+chi*(phi*psi+(1-phi)*(gamma+(sigma-gamma)*(1-2*alpha+alpha*alpha)))))*a_B +((1+chi)*(1-phi)*(sigma-gamma)*(2*alpha-alpha*alpha)/(1+chi*(phi*psi+((1-phi)*sigma))*(1+chi*(phi*chi+(1-phi)*(gamma+(sigma-gamma)*(1-2*alpha+alpha*alpha))))))*a_star; y_bar_star =((1+chi)*(phi*psi+(1-phi)*sigma)/(1+chi*(phi*psi+(1-phi)*sigma)))*a_star; P_bar_s = P_bar_star - ((1+chi)/(1+chi*(phi*psi+(1-phi)*(gamma+(sigma-gamma)*(1-2*alpha+alpha*alpha)))))*(a_s - a_star); P_bar_B = P_bar_star -((1+chi)/(1+chi*(phi*psi+(1-phi)*(gamma+(sigma-gamma)*(1-2*alpha+alpha*alpha)))))*(a_B - a_star); tt_bar_s = (P_bar_star - P_bar_s); tt_bar_B =(1/n)*(P_bar_star- P_bar_B); tau_bar_s = -(1-tau_s)*muchapo_w_s; tau_bar_B = -(1-tau_B)*muchapo_w_B; r_bar_star =((1+chi)/(1+chi*(phi*psi+(1-phi)*sigma)))*(a_star(+1)- a_star); resid; steady; check; //--------------------------------------------------------------------- // Partie 7: Declaration de distributions à priori des paramètres //--------------------------------------------------------------------- shocks; var a_s = 0.5; var a_B = 0.5; var a_star = 0.25; end; stoch_simul (linear,order=1 , irf=100);