%|>---------------------------------------------------------------------<|%
%|> | DESEMPLEO EN UN MODELO NEOKEYNESIANO | <|%
%|>---------------------------------------------------------------------<|%
close all;
warning off;
%|>-----------------------| Variables Endógenas |-----------------------<|%
var
C $\C$ (longname='Consumo')
RK $\RK$ (longname='Tasa de alquiler del capital')
RR $\RR$ (longname='Tasa de interés real')
W $\W$ (longname='Salario real')
N $\N$ (longname='Tasa de empleo')
Y $\Y$ (longname='Producto')
K $\K$ (longname='Capital')
Q $\Q$ (longname='Q de Tobin')
I $\I$ (longname='Inversión')
LAMBDA $\LAMBDA$ (longname='Utilidad marginal del consumo')
R $\R$ (longname='Tasa de interes nominal')
PI $\PI$ (longname='Inflación')
MC $\MC$ (longname='Costo marginal neto')
V $\V$ (longname='Vacantes')
G $\G$ (longname='Gasto público')
A $\A$ (longname='PTF exógena ')
U $\U$ (longname='Tasa de desempleo')
F $\F$ (longname='Tasa de búsqueda de empleo') %Probabilidad de encontrar un trabajo
GAMMA $\GAMMA$ (longname='Tasa de ocupación de vacantes') %Probabilidad de ocupar una vacante
THETA $\THETA$ (longname='Rigidez del mercado laboral')
OME $\OME$ (longname='Número de nuevos emparejamientos')
S $\S$ (longname='El excedente total de un trabajo')
OMEGA $\OMEGA$ (longname='Valor de las actividades no laborales')
VERIFY $\VERIFY$ (longname='Variable de verificación') % check variable, it should be always 0 (equation 26 in Monacelli et al.)
CLOG $\CLOG$ (longname='Logaritmo del Consumo')
WLOG $\WLOG$ (longname='Logaritmo del Salario real')
KLOG $\KLOG$ (longname='Logaritmo del capital')
ILOG $\ILOG$ (longname='Logaritmo de la inversión')
YLOG $\YLOG$ (longname='Logaritmo del producto')
THETALOG $\THETALOG$ (longname='Logaritmo de la rigidez del mercado laboral')
NLOG $\NLOG$ (longname='Logaritmo de la tasa de empleo')
ULOG $\NLOG$ (longname='Logaritmo de la tasa de desempleo')
PILOG $\PILOG$ (longname='Logaritmo de la tasa de inflación')
VLOG $\VLOG$ (longname='Logaritmo de las vacantes')
RLOG $\RLOG$ (longname='Logaritmo de la tasa de interés nominal')
GLOG $\GLOG$ (longname='Logaritmo del gasto público')
;
%|>-----------------------| Variables Exógenas |-----------------------<|%
varexo
VA $\VA$ (longname='Perturbación de productividad -PTF')
VG $\VG$ (longname='Perturbación de gasto público')
VM $\VM$ (longname='Perturbación de política monetaria')
;
%|>-----------------------|Parámetros del Modelo|-----------------------<|%
parameters
beta $\beta$ (longname='Factor de descuento')
alpha $\alpha$ (longname='Elasticidad de la producción con respecto al capital')
epsilon $\epsilon$ (longname='Elasticidad de sustitución entre bienes diferenciados')
delta $\delta$ (longname='Tasa de depreciación del capital')
sigma $\sigma$ (longname='Complementariedad trabajo/consumo') %if 1, consumption and labor are separable
mu $\mu$ (longname='Elasticidad de los emparejamientos con respecto al desempleo')
eta $\eta$ (longname='Poder de negociación de los trabajadores')
rho $\rho$ (longname='Tasa de permanencia laboral')
kappaV $\kappaV$ (longname='Costo de una vacante')
nu $\nu$ (longname='Desutilidad del trabajo')
Omebar $\Omebar$ (longname='Eficiencia de la función de emparejamiento')
zeta $\zeta$ (longname='Valor marginal del ocio sobre el producto marginal del capital')
%| Valores de Estado Estacionario
gss $\gss$ (longname='El gasto en el Steady State - SS')
ass $\ass$ (longname='La PTF en el Steady State - SS')
piss $\piss$ (longname='La inflación en el Steady State - SS')
rss $\rss$ (longname='La tasa de interés nominal en el Steady State - SS')
yss $\yss$ (longname='El producto en el Steady State - SS')
thetass $\thetass$ (longname='La rigidez laboral en el Steady State - SS')
fss $\fss$ (longname='La tasa de búsqueda de empleo en el Steady State - SS')
kappaP $\kappaP$ (longname='Costo de ajuste de precios')
kappaI $\kappaI$ (longname='Costo de ajuste de inversión') % (as in CEE). If 0, q is constant
phipi $\phipi$ (longname='Respuesta de la política monetaria a la inflación')
phiy $\phiy$ (longname='Respuesta de la política monetaria al producto')
rhoa $\rhoa$ (longname='Persistencia de la PTF')
rhog $\rhog$ (longname='Persistencia del gasto público')
rhom $\rhom$ (longname='Inercia de la política monetaria')
;
%|>------------------------| Modelo No Lineal |-----------------------<|%
model;
%| HOGARES
[name='1-Elección óptima de consumo']
LAMBDA =((1+(sigma-1)*nu*N)/C)^(sigma);
[name='2-Elección óptima de bonos']
1=beta*LAMBDA(1)/LAMBDA*R/PI(1);
[name='3-Elección óptima de capital']
1=beta*LAMBDA(1)/LAMBDA*(RK(1)+(1-delta)*Q(1))/Q;
[name='4-Dinámica del capital']
K=(1-delta)*K(-1)+(1-kappaI/2*(I/I(-1)-1)^2)*I;
[name='5-Elección óptima de inversión']
1=Q*(1-kappaI/2*(I/I(-1)-1)^2-kappaI*(I/I(-1)-1)*I/I(-1))+kappaI*beta*LAMBDA(1)/LAMBDA*Q(1)*(I(1)/I-1)*(I(1)/I)^2; %
[name='6-Definición de la tasa de interés real']
RR=R/PI(1);
%| FIRMAS
[name='7-Dinámica del empleo en las firmas']
N=rho*N(-1)+OME;
[name='8-Función de producción de las firmas']
Y=A*K(-1)^(alpha)*N^(1-alpha);
[name='9-Demanda de capital de las firmas']
alpha*MC*Y=RK*K(-1);
[name='10-Condición de precios']
(PI-piss)*PI=beta*(LAMBDA(1)/LAMBDA*Y(1)/Y*PI(1)*(PI(1)-piss))+epsilon/kappaP*(MC-(epsilon-1)/epsilon);
[name='11-Elección óptima de vacantes y de empleo']
kappaV/GAMMA=(1-alpha)*MC*Y/N-W+beta*rho*kappaV/GAMMA(1)*LAMBDA(1)/LAMBDA;
%| NEGOCIACIÓN DE LOS SALARIOS
[name='12-Salario real']
W=eta*(1-alpha)*MC*Y/N+(1-eta)*sigma*nu*LAMBDA^(-1/sigma)+eta*beta*kappaV*LAMBDA(1)/LAMBDA*THETA(+1);
%| VACIADO DE MERCADOS
[name='13-Vaciado de Mercados']
Y=C+I+G+kappaV*V+(kappaP/2*(PI-piss)^2)*Y;
%| POLITICA MONETARIA
[name='14-Regla de política monetaria']
R/(rss)=((PI/piss)^(phipi)*(Y/Y(-1))^(phiy))^(1-rhom)*(R(-1)/rss)^(rhom)*exp(VM);
%| PERTURBACIONES
[name='15-Comportamiento del shock a la PTF']
log(A)=(1-rhoa)*log(ass)+rhoa*log(A(-1))+VA;
[name='16-Comportamiento del shock al gasto']
log(G)=(1-rhog)*log(gss)+rhog*log(G(-1))+VG;
%G=(1-rhog)*gss+rhog*G(-1)+VG;
%| VARIABLES AUXILIARES
[name='17-Función de emparejamiento']
OME=Omebar*U^(mu)*V^(1-mu);
[name='18-Tasa de desempleo']
U=1-N(-1);
[name='19-Probabilidad de ocupar una vacante']
GAMMA=OME/V;
[name='20-Rigidez del mercado laboral']
THETA=V/U;
[name='21-EL excedente total de un trabajo.']
S=(1-alpha)*MC*Y/N-OMEGA+beta*LAMBDA(1)/LAMBDA*S(1)*(rho-eta*F(1));
[name='22-Probabilidad de encontrar trabajo']
F=OME/U;
[name='23-Valor de las actividades no laborales']
OMEGA=sigma*nu*LAMBDA^(-1/sigma);
[name='24-Verificación']
VERIFY=kappaV/GAMMA-(1-eta)*S;
%| LOG - VARIABLES
[name='25-Logaritmo del consumo']
CLOG=log(C)*100;
[name='26-Logaritmo del capital']
KLOG=log(K)*100;
[name='27-Logaritmo de la inversión']
ILOG=log(I)*100;
[name='28-Logaritmo del producto']
YLOG=log(Y)*100;
[name='29-Logaritmo de la rigidez del mercado laboral']
THETALOG=log(THETA)*100;
[name='30-Logaritmo del salario']
WLOG=log(W)*100;
[name='31-Logaritmo del empleo']
NLOG=log(N)*100;
[name='32-Logaritmo del desempleo']
ULOG=log(U)*100;
[name='33-Logaritmo de la inflación']
PILOG=log(PI)*100;
[name='34-Logaritmo de las vacantes']
VLOG=log(V)*100;
[name='35-Logaritmo de la tasa de interés nominal']
RLOG=log(R)*100;
[name='36-Logaritmo del gasto público']
GLOG=log(G)*100;
end;
%|>---------------|Asignación de Valores a los Parámetros|--------------<|%
beta=0.99^(1/3);
alpha=0.33;
epsilon=7.25;
delta=0.025/3;
sigma=1.0; % complementarity parameter: if 1, consumption and labor are separable
mu=0.5;
eta=mu;
rho=0.965;
%|>------------------|Modelo en el Estado Estacionario|-----------------<|%
steady_state_model;
gshare=0.2; % Partiicipación del gasto público en el SS
A=1; % TFP en SS
PI=1; % Inflación objetivo
F=0.45; % Probabilidad de encontrar trabajo en SS
THETA=0.5; % Rigidez del mercado laboral en SS
zeta=0.9; % Valor marginal del ocio sobre el producto marginal del capital en SS
Omebar=F/(THETA^(1-mu)); % Eficiencia de la función de emparejamiento en SS
GAMMA=Omebar*THETA^(-mu); % Probabilidad de ocupar una vacante en SS
N=F/(1+F-rho); % Tasa de empleo en SS
U=1-N; % Tasa de desempleo en SS
V=THETA*U; % Vacantes en SS
OME=Omebar*U^(mu)*V^(1-mu); % Emparejamientos en SS
RR=1/beta; % Tasa de interes real en SS
R=PI/beta; % Tasa de interés nominal
Q=1; % Q de tobin - Valor marginal de la inversión (in terms of lambda)
RK=1/beta-(1-delta); % Tasa de alquiler del capital en SS
MC=(epsilon-1)/epsilon; % Costo marginal en SS
K=N*(A*alpha*MC/RK)^(1/(1-alpha)); % Capital en SS
Y=A*K^(alpha)*N^(1-alpha); % Producto en SS
I=delta*K; % Inversión en SS
G=gshare*Y; % Gasto público en SS
mpl=(1-alpha)*MC*Y/N; % Producto marginal del trabajo en SS
kappaV=(1-eta)*(1-zeta)*mpl/((1/GAMMA)*(1-beta*rho)+beta*eta*THETA); %
W=mpl*(eta+(1-eta)*zeta)+eta*beta*kappaV*THETA; % Salario real en SS
C=Y-I-G-kappaV*V; % Consumo en SS
nu=zeta*mpl/(sigma*C-zeta*mpl*(sigma-1)*N); % Desutilidad del trabajo en SS
% Valores de SS
yss=Y;
gss=G;
ass=A;
piss=PI;
rss=R;
thetass=THETA;
fss=F;
% Parámetros que no afectan el SS
phipi=1.5^(1/3);
phiy=(0.5/4)^(1/3);
kappaI=3.24*3;
rhoa=0.95^(1/3);
rhog=0.90^(1/3);
rhom=0;
calvo=11/12;
kappaP=(epsilon-1)*calvo/(piss^2*(1-calvo)*(1-beta*calvo));
LAMBDA=((1+(sigma-1)*nu*N)/C)^(sigma);
CLOG=log(C)*100;
KLOG=log(K)*100;
ILOG=log(I)*100;
YLOG=log(Y)*100;
THETALOG=log(THETA)*100;
WLOG=log(W)*100;
NLOG=log(N)*100;
ULOG=log(U)*100;
PILOG=log(PI)*100;
VLOG=log(V)*100;
RLOG=log(R)*100;
GLOG=log(G)*100;
S=((1-alpha)*MC*Y/N-sigma*nu*LAMBDA^(-1/sigma))/(1-beta*(rho-eta*F));
OMEGA=sigma*nu*LAMBDA^(-1/sigma);
VERIFY=0;
end;
%|>------------------|Cálculo del Estado Estacionario|------------------<|%
steady;
%|>--------------------| Cálculo de los residuales |--------------------<|%
resid;
%|>-------------------------| Verificaciones |-------------------------<|%
check;
%|>------------------|Involucración de Perturbaciones|------------------<|%
shocks;
var VA; stderr 0.01;
var VM; stderr 0.0025;
var VG; stderr 0.0481; % 1% of GDP in flexible-price model
end;
%|>-------------------| Funciones Impulso Respuesta |-------------------<|%
stoch_simul(irf=40,order=1,noprint)
YLOG
CLOG
ILOG
ULOG
NLOG
VLOG
THETALOG
PILOG
WLOG
RLOG
GLOG
VERIFY;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure('units','normalized','outerposition',[0 0 0.57 0.57])
subplot(3,2,1)
plot(YLOG_VA,'b','LineWidth',2) %b-blue
hold on
plot(PILOG_VA,'r','LineWidth',2)
hold off
grid on
leg = legend('Producto','Inflación');
set(legend,'Interpreter','Latex');
leg.FontSize = 11;
xlim([1 40])
ylim([1.1*min(min([YLOG_VA PILOG_VA])) 1.1*max(max([YLOG_VA PILOG_VA]))]);
tit = title('IRF - Producto e Inflación','Interpreter','latex');
tit.FontSize = 14;
ylabel('Desviación (%)', 'Fontsize', 11,'Fontname', 'Garamond')
xlabel('Tiempo', 'Fontsize', 11,'Fontname', 'Garamond')
subplot(3,2,2)
plot(CLOG_VA,'b','LineWidth',2) %b-blue
hold on
plot(ILOG_VA,'r','LineWidth',2)
hold off
grid on
leg = legend('Consumo','Inversión');
set(legend,'Interpreter','Latex');
leg.FontSize = 11;
xlim([1 40])
ylim([1.1*min(min([CLOG_VA ILOG_VA])) 1.1*max(max([CLOG_VA ILOG_VA]))]);
tit = title('IRF - Consumo e Inversión','Interpreter','latex');
tit.FontSize = 14;
ylabel('Desviación (%)', 'Fontsize', 11,'Fontname', 'Garamond')
xlabel('Tiempo', 'Fontsize', 11,'Fontname', 'Garamond')
subplot(3,2,3)
plot(ULOG_VA,'b','LineWidth',2) %b-blue
hold on
plot(VLOG_VA,'r','LineWidth',2)
hold off
grid on
leg = legend('Desempleo','Vacantes');
set(legend,'Interpreter','Latex');
leg.FontSize = 11;
xlim([1 40])
ylim([1.1*min(min([ULOG_VA VLOG_VA])) 1.1*max(max([ULOG_VA VLOG_VA]))]);
tit = title('IRF - Desempleo y Vacantes','Interpreter','latex');
tit.FontSize = 14;
ylabel('Desviación (%)', 'Fontsize', 11,'Fontname', 'Garamond')
xlabel('Tiempo', 'Fontsize', 11,'Fontname', 'Garamond')
subplot(3,2,4)
plot(NLOG_VA,'b','LineWidth',2) %b-blue
hold on
plot(WLOG_VA,'r','LineWidth',2)
hold off
grid on
leg = legend('Empleo','Salario Real');
set(legend,'Interpreter','Latex');
leg.FontSize = 11;
xlim([1 40])
ylim([1.1*min(min([NLOG_VA WLOG_VA])) 1.1*max(max([NLOG_VA WLOG_VA]))]);
tit = title('IRF - Empleo y Salario Real','Interpreter','latex');
tit.FontSize = 14;
ylabel('Desviación (%)', 'Fontsize', 11,'Fontname', 'Garamond')
xlabel('Tiempo', 'Fontsize', 11,'Fontname', 'Garamond')
subplot(3,2,5)
plot(THETALOG_VA,'b','LineWidth',2)
grid on
xlim([1 40])
ylim([1.1*min(min([THETALOG_VA])) 1.1*max(max([THETALOG_VA]))]);
tit = title('IRF - Rigidez del Mercado Laboral','Interpreter','latex');
tit.FontSize = 14;
ylabel('Desviación (%)', 'Fontsize', 11,'Fontname', 'Garamond')
xlabel('Tiempo', 'Fontsize', 11,'Fontname', 'Garamond')
subplot(3,2,6)
plot(RLOG_VA,'b','LineWidth',2)
grid on
xlim([1 40])
ylim([1.1*min(min([RLOG_VA])) 1.1*max(max([RLOG_VA]))]);
tit = title('IRF - Tasa de Política Monetaria','Interpreter','latex');
tit.FontSize = 14;
ylabel('Desviación (%)', 'Fontsize', 11,'Fontname', 'Garamond')
xlabel('Tiempo', 'Fontsize', 11,'Fontname', 'Garamond')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure('units','normalized','outerposition',[0 0 0.57 0.57])
subplot(3,2,1)
plot(YLOG_VG,'b','LineWidth',2) %b-blue
hold on
plot(PILOG_VG,'r','LineWidth',2)
hold off
grid on
leg = legend('Producto','Inflación');
set(legend,'Interpreter','Latex');
leg.FontSize = 11;
xlim([1 40])
ylim([1.1*min(min([YLOG_VG PILOG_VG])) 1.1*max(max([YLOG_VG PILOG_VG]))]);
tit = title('IRF - Producto e Inflación','Interpreter','latex');
tit.FontSize = 14;
ylabel('Desviación (%)', 'Fontsize', 11,'Fontname', 'Garamond')
xlabel('Tiempo', 'Fontsize', 11,'Fontname', 'Garamond')
subplot(3,2,2)
plot(CLOG_VG,'b','LineWidth',2) %b-blue
hold on
plot(ILOG_VG,'r','LineWidth',2)
plot(GLOG_VG,'k','LineWidth',2)
hold off
grid on
leg = legend('Consumo','Inversión','Gasto');
set(legend,'Interpreter','Latex');
leg.FontSize = 11;
xlim([1 40])
ylim([1.1*min(min([CLOG_VG ILOG_VG GLOG_VG])) 1.1*max(max([CLOG_VG ILOG_VG GLOG_VG]))]);
tit = title('IRF - Consumo e Inversión','Interpreter','latex');
tit.FontSize = 14;
ylabel('Desviación (%)', 'Fontsize', 11,'Fontname', 'Garamond')
xlabel('Tiempo', 'Fontsize', 11,'Fontname', 'Garamond')
subplot(3,2,3)
plot(ULOG_VG,'b','LineWidth',2) %b-blue
hold on
plot(VLOG_VG,'r','LineWidth',2)
hold off
grid on
leg = legend('Desempleo','Vacantes');
set(legend,'Interpreter','Latex');
leg.FontSize = 11;
xlim([1 40])
ylim([1.1*min(min([ULOG_VG VLOG_VG])) 1.1*max(max([ULOG_VG VLOG_VG]))]);
tit = title('IRF - Desempleo y Vacantes','Interpreter','latex');
tit.FontSize = 14;
ylabel('Desviación (%)', 'Fontsize', 11,'Fontname', 'Garamond')
xlabel('Tiempo', 'Fontsize', 11,'Fontname', 'Garamond')
subplot(3,2,4)
plot(NLOG_VG,'b','LineWidth',2) %b-blue
hold on
plot(WLOG_VG,'r','LineWidth',2)
hold off
grid on
leg = legend('Empleo','Salario Real');
set(legend,'Interpreter','Latex');
leg.FontSize = 11;
xlim([1 40])
ylim([1.1*min(min([NLOG_VG WLOG_VG])) 1.1*max(max([NLOG_VG WLOG_VG]))]);
tit = title('IRF - Empleo y Salario Real','Interpreter','latex');
tit.FontSize = 14;
ylabel('Desviación (%)', 'Fontsize', 11,'Fontname', 'Garamond')
xlabel('Tiempo', 'Fontsize', 11,'Fontname', 'Garamond')
subplot(3,2,5)
plot(THETALOG_VG,'b','LineWidth',2)
grid on
xlim([1 40])
ylim([1.1*min(min([THETALOG_VG])) 1.1*max(max([THETALOG_VG]))]);
tit = title('IRF - Rigidez del Mercado Laboral','Interpreter','latex');
tit.FontSize = 14;
ylabel('Desviación (%)', 'Fontsize', 11,'Fontname', 'Garamond')
xlabel('Tiempo', 'Fontsize', 11,'Fontname', 'Garamond')
subplot(3,2,6)
plot(RLOG_VG,'b','LineWidth',2)
grid on
xlim([1 40])
ylim([1.1*min(min([RLOG_VG])) 1.1*max(max([RLOG_VG]))]);
tit = title('IRF - Tasa de Política Monetaria','Interpreter','latex');
tit.FontSize = 14;
ylabel('Desviación (%)', 'Fontsize', 11,'Fontname', 'Garamond')
xlabel('Tiempo', 'Fontsize', 11,'Fontname', 'Garamond')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure('units','normalized','outerposition',[0 0 0.57 0.57])
subplot(3,2,1)
plot(YLOG_VM,'b','LineWidth',2) %b-blue
hold on
plot(PILOG_VM,'r','LineWidth',2)
hold off
grid on
leg = legend('Producto','Inflación');
set(legend,'Interpreter','Latex');
leg.FontSize = 11;
xlim([1 40])
ylim([1.1*min(min([YLOG_VM PILOG_VM])) 1.1*max(max([YLOG_VM PILOG_VM]))]);
tit = title('IRF - Producto e Inflación','Interpreter','latex');
tit.FontSize = 14;
ylabel('Desviación (%)', 'Fontsize', 11,'Fontname', 'Garamond')
xlabel('Tiempo', 'Fontsize', 11,'Fontname', 'Garamond')
subplot(3,2,2)
plot(CLOG_VM,'b','LineWidth',2) %b-blue
hold on
plot(ILOG_VM,'r','LineWidth',2)
hold off
grid on
leg = legend('Consumo','Inversión');
set(legend,'Interpreter','Latex');
leg.FontSize = 11;
xlim([1 40])
ylim([1.1*min(min([CLOG_VM ILOG_VM])) 1.1*max(max([CLOG_VM ILOG_VM]))]);
tit = title('IRF - Consumo e Inversión','Interpreter','latex');
tit.FontSize = 14;
ylabel('Desviación (%)', 'Fontsize', 11,'Fontname', 'Garamond')
xlabel('Tiempo', 'Fontsize', 11,'Fontname', 'Garamond')
subplot(3,2,3)
plot(ULOG_VM,'b','LineWidth',2) %b-blue
hold on
plot(VLOG_VM,'r','LineWidth',2)
hold off
grid on
leg = legend('Desempleo','Vacantes');
set(legend,'Interpreter','Latex');
leg.FontSize = 11;
xlim([1 40])
ylim([1.1*min(min([ULOG_VM VLOG_VM])) 1.1*max(max([ULOG_VM VLOG_VM]))]);
tit = title('IRF - Desempleo y Vacantes','Interpreter','latex');
tit.FontSize = 14;
ylabel('Desviación (%)', 'Fontsize', 11,'Fontname', 'Garamond')
xlabel('Tiempo', 'Fontsize', 11,'Fontname', 'Garamond')
subplot(3,2,4)
plot(NLOG_VM,'b','LineWidth',2) %b-blue
hold on
plot(WLOG_VM,'r','LineWidth',2)
hold off
grid on
leg = legend('Empleo','Salario Real');
set(legend,'Interpreter','Latex');
leg.FontSize = 11;
xlim([1 40])
ylim([1.1*min(min([NLOG_VM WLOG_VM])) 1.1*max(max([NLOG_VM WLOG_VM]))]);
tit = title('IRF - Empleo y Salario Real','Interpreter','latex');
tit.FontSize = 14;
ylabel('Desviación (%)', 'Fontsize', 11,'Fontname', 'Garamond')
xlabel('Tiempo', 'Fontsize', 11,'Fontname', 'Garamond')
subplot(3,2,5)
plot(THETALOG_VM,'b','LineWidth',2)
grid on
xlim([1 40])
ylim([1.1*min(min([THETALOG_VM])) 1.1*max(max([THETALOG_VM]))]);
tit = title('IRF - Rigidez del Mercado Laboral','Interpreter','latex');
tit.FontSize = 14;
ylabel('Desviación (%)', 'Fontsize', 11,'Fontname', 'Garamond')
xlabel('Tiempo', 'Fontsize', 11,'Fontname', 'Garamond')
subplot(3,2,6)
plot(RLOG_VM,'b','LineWidth',2)
grid on
xlim([1 40])
ylim([1.1*min(min([RLOG_VM])) 1.1*max(max([RLOG_VM]))]);
tit = title('IRF - Tasa de Política Monetaria','Interpreter','latex');
tit.FontSize = 14;
ylabel('Desviación (%)', 'Fontsize', 11,'Fontname', 'Garamond')
xlabel('Tiempo', 'Fontsize', 11,'Fontname', 'Garamond')