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