COMANDERA R (VAR) Flashcards
Que librerias ocupamos para modelos VAR?
- library(tidyverse)
- library(dynlm)
- library(forecast)
- library(tseries)
- library(vars)
- library(lmtest)
- library(quantmod)
- library(mFilter)
- library(nlme)
- library(graphics)
M2.ts= ts(Base, start = c(2001,12), freq = 12)
se crea una base de datos de serie temporal que comienza desde dic 2001 con una frecuencia mensual
m2a.ts= ts(m2.ts [,2], start = c(2001,12), frequency = 12)
IPC.ts= ts(m2.ts [,3], start= c(2001,12), frequency = 12)
indicamos con que columnas trabajaremos las variables de interés y transformarlas a tiempo
lnm2= log(m2a.ts)
lnIPC= log(IPC.ts)
Transformamos las variables en log
nsdiffs(lnm2)
nsdiffs(lnIPC)
Transformamos las variables en el mismo numero de diferencias de la que tenga mayor numero de cointegraciones para el modelo VAR (el IPC tiene dos cointegraciones)
ts.plot(lnm2, lnIPC, col = c(“blue”,”red”), main= “cointegracion de M2 y el IPC en log”)
Realizar un grafico para determinar su paralelismo respecto a su cointegacion
dlnIPC= diff(lnIPC)
d2lnIPC= diff(dlnIPC)
Dlnm2= diff(lnm2)
d2lnm2=diff(Dlnm2)
muestra de doble diferenciación por la variable del IPC que por teoría se diferencia dos veces con el objetivo de volver estacionaria la serie
ts.plot(d2lnm2, d2lnIPC,col= c(“blue”, “green”), main= “SSSS”)
Realizamos un gráfico con las variables ya diferenciadas
adf.test(d2lnm2)
adf.test(d2lnIPC)
Utilizamos métodos formales para verificar estacionariedad
windows()
abre una ventana donde se pueden ver mejor los graficos
1.grangertest(D2lnM~D2lnIPC, order = 2)
2.(grangertest(D2lnIPC~D2lnM, order = 2))
Prueba de hipótesis para verificar cointegracion
1—Ho: La variable X no causa cambios en Y
Ha: La variable X causa cambios en Y
2—-Ho: La variable Y no causa cambios en X
Ha: La variable Y causa cambios en X
Vard2lnm2= ts(d2lnm2, start= c(2001,12), frequency=12)
vard2lnIPC=ts(d2lnIPC, start= c(2001,12), freq=12)
Crear nuevamente objetos con las variables diferenciadas
homorfo= cbind(d2lnm2, d2lnIPC)
Crear una base de datos con un solo objeto
VARselect(homorfo, lag.max= 12)
El objeto creado lo sometemos a selección de VAR, para generar criterios de información(# Con los cuatro criterios el orden de los rezagos es 11)
homorfo1<- VAR(homorfo, p=11)
summary(homorfo1)
Revisión de las raíces polimoniales del modelo, si este es menor que 1 se esta cumpliendo la condición de estabilidad del modelo
predict(homorfo1)
layout(1:2)
plot(homorfo1)
predict: crear variable a pronosticar, layout: ?, plot, grafico para ver la distribución de los errores
serial<- serial.test(homorfo1, lags.pt = 11, type= “PT.asymptotic”)
serial
prueba de especificación del modelo
Ho: Los residuos no estan correlacionados; valor p > 0.05
Ha; Los residuos si estan correlacionados
normalidad=normality.test(var1)
normalidad$jb.mul
prueba de normalidad
Ho; Los residuos se distribuyen normalmente
Ha; Los residuos no se distribuyen normalme
Arch1=arch.test(homorfo1, lags.multi = 11)
arch
Prueba de varianza residual:
Ho; La varianza de los errores es constante
Ha; La varianza de los errores no es constante
IRF<- irf(homorfo1, response = homorfo1$d2lnIPC, nahead=8, boot = TRUE)
plot(IRF)
Analisis impulso respuesta
¿Cuáles son los pasos para realizar un pronóstico VAR?
- Transformar base de datos a tiempo, ver en qué columnas están las variables de interés y transformarlas a tiempo
- Transformamos las variables en log
- Transformamos variables en diferencias para el modelo VAR (hay variables con mas de 1 diff)
- Realizamos un gráfico con las variables ya diferenciadas para ver si la serie es estacionaria
- Utilizamos métodos formales para verificar estacionariedad
- Utilizamos métodos formales para verificar cointegracion
- Generamos un nuevo objeto con las segundas diferencias
- Crear nuevamente objetos con las variables diferenciadas
- Crear un data frame con un solo objeto
- El objeto creado lo sometemos a selección de VAR
# Con los cuatro criterios el orden de los rezagos es 11 - Revisión de las raíces polimoniales del modelo, si este es menor que 1 se esta cumpliendo la condición de estabilidad del model
- Pronostico del modelo
- Pruebas de especificación
- Pruebas de normalidad en los errores
- Prueba de varianza y covarianza residual
- Analisis impulso respuesta
- Grafico IR