Boot Camp 2023 - RIEF
SUMMER SCHOOL IN ECONOMICS AND FINANCE 2023
CODING ASSIGNMENT R-STUDIO
Docente: Christian Josser Alcarraz Valer ([Link]@[Link])
Instrucciones:
▪ Presente los resultados de cada pregunta a través de un informe (no olvide indicar
los nombres de los integrantes del grupo en alguna parte del informe).
▪ Adjunte los scripts empleados (archivos.R)
▪ Estructure el código ordenadamente en secciones y preguntas.
Finalmente, presente todo lo anterior a través de una carpeta comprimida indicando el
número del grupo. Ejemplo: Grupo_1.zip.
Sección I (8 puntos)
Mathematical Modelling of infectious disease
En el distrito Shiganshina, apareció un nuevo brote de virus por lo que el gobierno de Eldia
pretende explorar cuantitativamente cómo las campañas de vacunación preventiva influyen
en la población susceptible para ello se extiende el modelo SIR básico a la versión estática
del modelo SIR-V de Schlickeiser y Kröger(2021)1 el cual está definido por el siguiente
sistema de ecuaciones diferenciales.
𝑑𝑆(𝑡) 𝛽𝐼(𝑡)𝑆(𝑡)
=− − 𝑣𝑆(𝑡)
𝑑𝑡 𝑁
𝑑𝐼(𝑡) 𝛽𝐼(𝑡)𝑆(𝑡)
= − 𝛾𝐼(𝑡)
𝑑𝑡 𝑁
𝑑𝑅(𝑡)
= 𝛾𝐼(𝑡)
𝑑𝑡
𝑑𝑉(𝑡)
= 𝑣𝑆(𝑡)
𝑑𝑡
𝑁 = 𝑆(𝑡) + 𝐼(𝑡) + 𝑅(𝑡) + 𝑉(𝑡)
Donde:
▪ Población susceptible 𝑆(𝑡) : Son los individuos sin inmunidad frente al agente
infeccioso y que, por tanto, pueden infectarse.
▪ Población infectada 𝐼(𝑡): Son los individuos que están infectados en un momento
dado y pueden transmitir la infección a los individuos de la población susceptible
con los que entran en contacto.
▪ Población recuperada 𝑅(𝑡): Son los individuos inmunes a la infección por efecto de
su propio organismo (Inmunidad natural)
▪ Población Vacunada 𝑉(𝑡): Son los individuos inmunes a la infección por efecto de las
campañas de vacunación (Inmunidad artificial)
▪ 𝛽, 𝛾 , 𝑣: Son las tasas de infección, recuperación y vacunación, respectivamente
1 [Link]
Boot Camp 2023 - RIEF
Según informes estadísticos de un equipo de exploración. En promedio una persona
infectada transmite el agente infeccioso hacia la población susceptible cada 2 días, así
mismo, el agente infectado desarrolla una inmunidad natural luego de 8 días. Debido a que
el distrito se encuentra entre murallas y los recursos son limitados el gobierno estableció
una política de vacunación del 5% cada día.
En clases vimos varias extensiones del modelo SIR, pueden ayudarse modificando el código
Pandemic_Models.R para resolver las siguientes preguntas:
1. Para un periodo de 100 días y considerando los siguientes valores iniciales: 𝑆(0) =
999 , 𝐼(0) = 1, 𝑅(0) = 𝑉(0) = 0. Presente en un solo gráfico (usando el paquete
ggplot2) la evolución en el tiempo de las 4 variables principales del modelo. Describa
brevemente el comportamiento de las variables. (3 puntos)
2. Presente en un solo gráfico (usando el paquete ggplot2) tanto la curva de infección bajo
la política de vacunación del 5% y la curva de infección bajo el modelo SIR básico. ¿Qué
puede concluir respecto a las campañas de vacunación? (3 puntos)
3. Bajo el modelo SIR-V, en promedio cuántas personas evitaron infectarse gracias a la
política de vacunación del 5% (1 punto)
4. Bajo el modelo SIR-V. ¿En cuántos días se logra vacunar a casi al 50% de la población?
(1 punto)
Sección II (12 puntos)
Financial Econometrics: Fama French con Macro – Factores
Usando la información de “Data [Link]”, el cual contiene información de retornos
mensuales (en porcentaje), desde enero del 2000 a septiembre del 2022, de 12 sectores
que componen al índice Standard & Poor's 500, así como información de 5 factores.
Se le pide que estime el siguiente modelo de regresión lineal múltiple en términos de
retornos en exceso:
𝑟𝑖𝑡 − 𝑟𝑓𝑡 = 𝛼𝑖 + 𝛽1𝑖 (𝑟𝑚𝑡 − 𝑟𝑓𝑡 ) + 𝛽2𝑖 𝑆𝑀𝐵𝑡 + 𝛽3𝑖 𝐻𝐿𝑀𝑡 + 𝛽4𝑖 𝜇𝑡 + 𝛽5𝑖 𝜋𝑡 + 𝜖𝑖𝑡
𝜖𝑖𝑡 ~𝑁(0, 𝜎𝑖2 )
𝑖 = 1,3,4, … . , 12
𝑡 = 𝑒𝑛𝑒 − 2000 𝑎 𝑠𝑒𝑝 − 2022
Donde:
• 𝑟𝑖𝑡 : Retorno mensual del sector 𝑖
• 𝑟𝑓𝑡 : Tasa de rendimiento mensual de un bono soberano a 1 mes del tesoro americano.
•𝑟𝑚𝑡 : Rendimiento mensual del mercado (S&P 500).
• 𝜇𝑡 : Tasa de desempleo mensual de US ([Link]
•𝜋𝑡 : Inflación subyacente interanual de US([Link]
En clases vimos el modelo de Fama French (1993) de 3 factores, aquí se le han añadido 2
variables macroeconómicas, así mismo el tamaño muestral inicia desde el 2000, puede
ayudarse modificando el código 3Factor_Model.R para resolver las siguientes preguntas:
Boot Camp 2023 - RIEF
1. ¿Cuántos y cuáles son los sectores que presentan retornos anormales? Ayuda: En
cuántos sectores se cumple que 𝛼𝑖 es estadísticamente distinto de 0 al 5% de
significancia. (1 punto)
2. ¿Cuáles son los sectores defensivos respecto al mercado? Ayuda: Dado que debe correr
12 regresiones puede ordenar de mayor a menor 𝛽1𝑖 y obtener aquellos valores que son
menores a 1 (1 punto)
3. ¿Cuáles son los 3 sectores con mayor exposición al factor desempleo? Ayuda: Dado que
debe correr 12 regresiones puede ordenar de mayor a menor 𝛽4𝑖 y obtener los 3
mayores valores. (1 puntos)
4. Presente un gráfico de puntos (usando el paquete ggplot2) donde en el eje “X”
represente el sector y el eje “Y” sea la probabilidad de significancia asociada a cada
coeficiente del factor desempleo ¿Qué puede concluir respecto a la importancia de este
factor? Ayuda: Debería obtener el siguiente gráfico. (3 puntos)
Nota adicional: En Source indique su Grupo. Ejemplo: Grupo N°1
5. ¿Cuáles son los 3 sectores con mayor exposición al factor inflación? Ayuda: Dado que
debe correr 12 regresiones puede ordenar de mayor a menor 𝛽5𝑖 y obtener los 3
mayores valores. (1 punto)
6. Presente un gráfico de puntos (usando el paquete ggplot2) donde en el eje “X”
represente el sector y el eje “Y” sea la probabilidad de significancia asociada a cada
coeficiente del factor inflación ¿Qué puede concluir respecto a la importancia de este
factor? Ayuda: Debería obtener el siguiente gráfico. (3 puntos)
Boot Camp 2023 - RIEF
Nota adicional: En Source indique su Grupo. Ejemplo: Grupo N°1
7. En promedio cuanto explica el modelo de Fama French con Macro – factores los retornos
del mercado de renta variable US. ¿Cuáles son los 2 sectores donde el modelo presenta
el mejor y el peor ajuste? Ayuda: Dado que debe correr 12 regresiones puede ordenar
de mayor a menor los 12 R2 ajustados (2 puntos)
