# In[1]:
import os
import matplotlib.pyplot as plt
import pandas as pd
# Opcion que nos permitira graficar
pd.set_option('display.width', 5000)
pd.set_option('display.max_columns', 60)
usuario = os.getlogin()
# # Definimos los valores razonables para los parámetros estructurales
# In[2]:
kapa = 0.01
sigma = 2
rho = 0.05
pi = 0.02
a_pi = 0.8
a_y = 1.5
# # Definimos las desviaciones estandar y persistencias de los choques
# In[3]:
sigma_pi = 1
sigma_y = 1
sigma_i = 1
rho_pi = 0.1
rho_y = 0.1
rho_i = 0.1
# In[4]:
def calcular_impulso(coeficiente,persistencia,s,desv_estandar):
return (coeficiente*(persistencia**s))/desv_estandar
# In[5]:
def graficar_impulso(coeficiente,persistencia,num_periodos,desviacion,nombre):
eje_x = list[range[0,num_periodos]]
eje_y = []
for i in eje_x:
eje_y.append(calcular_impulso(coeficiente,persistencia,i,desviacion))
plt.Figure(figsize = (20,20), dpi = 1200)
plt.plot(eje_x,eje_y,"-")
plt.title(nombre)
plt.savefig(nombre+".pdf")
plt.close()
list(range(etc))
import os
import matplotlib.pyplot as plt
import pandas as pd
# Opcion que nos permitira graficar
pd.set_option('display.width', 5000)
pd.set_option('display.max_columns', 60)
usuario = os.getlogin()
# # Definimos los valores razonables para los parámetros estructurales
# In[2]:
kapa = 0.01
sigma = 2
rho = 0.05
pi = 0.02
a_pi = 0.8
a_y = 1.5
# # Definimos las desviaciones estandar y persistencias de los choques
# In[3]:
sigma_pi = 1
sigma_y = 1
sigma_i = 1
rho_pi = 0.1
rho_y = 0.1
rho_i = 0.1
# In[4]:
def calcular_impulso(coeficiente,persistencia,s,desv_estandar):
return (coeficiente*(persistencia**s))/desv_estandar
# In[5]:
def graficar_impulso(coeficiente,persistencia,num_periodos,desviacion,nombre):
eje_x = list(range(0,num_periodos))
eje_y = []
for i in eje_x:
eje_y.append(calcular_impulso(coeficiente,persistencia,i,desviacion))
plt.Figure(figsize = (20,20), dpi = 1200)
plt.plot(eje_x,eje_y,"-")
plt.title(nombre)
plt.savefig(nombre+".pdf")
plt.close()
# # Calculando coeficientes según resolución del ejercicio
# In[6]:
a_uno = (sigma*(1-rho_pi)+a_y)/((1-rho_pi)*(sigma*(1-rho_pi)+a_y)-kapa*rho_pi*(1-a_pi))
# In[7]:
b_uno = (rho_pi*(1-a_pi))/((1-rho_pi)*(sigma*(1-rho_pi)+a_y)-kapa*rho_pi*(1-a_pi))
# In[8]:
c_uno = (rho_pi*(a_y*sigma+a_pi*(1-rho_pi)))/((1-rho_pi)*(sigma*(1-rho_pi)+a_y)-kapa*rho_pi*(1-a_pi))
# # Graficando las funciones de impulso respuesta para el choque de inflacion
# In[9]:
graficar_impulso(a_uno,rho_pi,10,sigma_pi,"Función impulso respuesta brecha de inflacion vs choque de inflacion")
# In[10]:
graficar_impulso(b_uno,rho_pi,10,sigma_pi,"Función impulso respuesta brecha de producto vs choque de inflacion")
# In[11]:
graficar_impulso(c_uno,rho_pi,10,sigma_pi,"Función impulso respuesta brecha de tasa de interes vs choque de inflacion")
# # Calculano coeficientes relevantes para el choque de demanda
# In[12]:
a_dos = (-kapa*sigma)/((1-rho_y)*(a_y+sigma)+rho_y*kapa*(a_pi-1))
# In[13]:
b_dos = ((1-rho_y)*sigma)/((1-rho_y)*(a_y+sigma)+rho_y*kapa*(a_pi-1))
# In[14]:
c_dos = -((1-rho_y)*a_y*sigma+a_pi*rho_y*kapa*sigma)/((1-rho_y)*(a_y+sigma)+rho_y*kapa*(a_pi-1))
# In[15]:
graficar_impulso(a_dos,rho_pi,10,sigma_pi,"Función impulso respuesta brecha de inflacion vs choque de demanda")
# In[16]:
graficar_impulso(b_dos,rho_pi,10,sigma_pi,"Función impulso respuesta brecha de producto vs choque de demanda")
# In[17]:
graficar_impulso(c_dos,rho_pi,10,sigma_pi,"Función impulso respuesta brecha de tasa de interes vs choque de demanda")
# In[ ]:
import os
import matplotlib.pyplot as plt
import pandas as pd
# Opcion que nos permitira graficar
pd.set_option('display.width', 5000)
pd.set_option('display.max_columns', 60)
usuario = os.getlogin()
# # Definimos los valores razonables para los parámetros estructurales
# In[2]:
kapa = 0.01
sigma = 2
rho = 0.05
pi = 0.02
a_pi = 0.8
a_y = 1.5
# # Definimos las desviaciones estandar y persistencias de los choques
# In[3]:
sigma_pi = 1
sigma_y = 1
sigma_i = 1
rho_pi = 0.1
rho_y = 0.1
rho_i = 0.1
# In[4]:
def calcular_impulso(coeficiente,persistencia,s,desv_estandar):
return (coeficiente*(persistencia**s))/desv_estandar
# In[5]:
def graficar_impulso(coeficiente,persistencia,num_periodos,desviacion,nombre):
eje_x = list(range(0,num_periodos))
eje_y = []
for i in eje_x:
eje_y.append(calcular_impulso(coeficiente,persistencia,i,desviacion))
plt.Figure(figsize = (20,20), dpi = 1200)
plt.plot(eje_x,eje_y,"-")
plt.title(nombre)
plt.savefig(nombre+".pdf")
plt.close()
# # Calculando coeficientes según resolución del ejercicio
# In[6]:
a_uno = (sigma*(1-rho_pi)+a_y)/((1-rho_pi)*(sigma*(1-rho_pi)+a_y)-kapa*rho_pi*(1-a_pi))
# In[7]:
b_uno = (rho_pi*(1-a_pi))/((1-rho_pi)*(sigma*(1-rho_pi)+a_y)-kapa*rho_pi*(1-a_pi))
# In[8]:
c_uno = (rho_pi*(a_y*sigma+a_pi*(1-rho_pi)))/((1-rho_pi)*(sigma*(1-rho_pi)+a_y)-kapa*rho_pi*(1-a_pi))
# # Graficando las funciones de impulso respuesta para el choque de inflacion
# In[9]:
graficar_impulso(a_uno,rho_pi,10,sigma_pi,"Función impulso respuesta brecha de inflacion vs choque de inflacion")
# In[10]:
graficar_impulso(b_uno,rho_pi,10,sigma_pi,"Función impulso respuesta brecha de producto vs choque de inflacion")
# In[11]:
graficar_impulso(c_uno,rho_pi,10,sigma_pi,"Función impulso respuesta brecha de tasa de interes vs choque de inflacion")
# # Calculano coeficientes relevantes para el choque de demanda
# In[12]:
a_dos = (-kapa*sigma)/((1-rho_y)*(a_y+sigma)+rho_y*kapa*(a_pi-1))
# In[13]:
b_dos = ((1-rho_y)*sigma)/((1-rho_y)*(a_y+sigma)+rho_y*kapa*(a_pi-1))
# In[14]:
c_dos = -((1-rho_y)*a_y*sigma+a_pi*rho_y*kapa*sigma)/((1-rho_y)*(a_y+sigma)+rho_y*kapa*(a_pi-1))
# In[15]:
graficar_impulso(a_dos,rho_pi,10,sigma_pi,"Función impulso respuesta brecha de inflacion vs choque de demanda")
# In[16]:
graficar_impulso(b_dos,rho_pi,10,sigma_pi,"Función impulso respuesta brecha de producto vs choque de demanda")
# In[17]:
graficar_impulso(c_dos,rho_pi,10,sigma_pi,"Función impulso respuesta brecha de tasa de interes vs choque de demanda")
# In[ ]:
#-*- coding:utf -8 -*-
import os
import matplotlib.pyplot as plt
import pandas as pd
# Opcion que nos permitira graficar
pd.set_option('display.width', 5000)
pd.set_option('display.max_columns', 60)
usuario = os.getlogin()
# # Definimos los valores razonables para los parámetros estructurales
# In[2]:
kapa = 0.01
sigma = 2
rho = 0.05
pi = 0.02
a_pi = 0.8
a_y = 1.5
# # Definimos las desviaciones estandar y persistencias de los choques
# In[3]:
sigma_pi = 1
sigma_y = 1
sigma_i = 1
rho_pi = 0.1
rho_y = 0.1
rho_i = 0.1
# In[4]:
def calcular_impulso(coeficiente,persistencia,s,desv_estandar):
return (coeficiente*(persistencia**s))/desv_estandar
# In[5]:
def graficar_impulso(coeficiente,persistencia,num_periodos,desviacion,nombre):
eje_x = list(range(0,num_periodos))
eje_y = []
for i in eje_x:
eje_y.append(calcular_impulso(coeficiente,persistencia,i,desviacion))
plt.Figure(figsize = (20,20), dpi = 1200)
plt.plot(eje_x,eje_y,"-")
plt.title(nombre)
plt.savefig(nombre+".pdf")
plt.close()
# # Calculando coeficientes según resolución del ejercicio
# In[6]:
a_uno = (sigma*(1-rho_pi)+a_y)/((1-rho_pi)*(sigma*(1-rho_pi)+a_y)-kapa*rho_pi*(1-a_pi))
# In[7]:
b_uno = (rho_pi*(1-a_pi))/((1-rho_pi)*(sigma*(1-rho_pi)+a_y)-kapa*rho_pi*(1-a_pi))
# In[8]:
c_uno = (rho_pi*(a_y*sigma+a_pi*(1-rho_pi)))/((1-rho_pi)*(sigma*(1-rho_pi)+a_y)-kapa*rho_pi*(1-a_pi))
# # Graficando las funciones de impulso respuesta para el choque de inflacion
# In[9]:
graficar_impulso(a_uno,rho_pi,10,sigma_pi,"Función impulso respuesta brecha de inflacion vs choque de inflacion")
# In[10]:
graficar_impulso(b_uno,rho_pi,10,sigma_pi,"Función impulso respuesta brecha de producto vs choque de inflacion")
# In[11]:
graficar_impulso(c_uno,rho_pi,10,sigma_pi,"Función impulso respuesta brecha de tasa de interes vs choque de inflacion")
# # Calculano coeficientes relevantes para el choque de demanda
# In[12]:
a_dos = (-kapa*sigma)/((1-rho_y)*(a_y+sigma)+rho_y*kapa*(a_pi-1))
# In[13]:
b_dos = ((1-rho_y)*sigma)/((1-rho_y)*(a_y+sigma)+rho_y*kapa*(a_pi-1))
# In[14]:
c_dos = -((1-rho_y)*a_y*sigma+a_pi*rho_y*kapa*sigma)/((1-rho_y)*(a_y+sigma)+rho_y*kapa*(a_pi-1))
# In[15]:
graficar_impulso(a_dos,rho_pi,10,sigma_pi,"Función impulso respuesta brecha de inflacion vs choque de demanda")
# In[16]:
graficar_impulso(b_dos,rho_pi,10,sigma_pi,"Función impulso respuesta brecha de producto vs choque de demanda")
# In[17]:
graficar_impulso(c_dos,rho_pi,10,sigma_pi,"Función impulso respuesta brecha de tasa de interes vs choque de demanda")
# In[ ]: