#install.packages(“vars”) library(vars) library(fpp2)

Carga y visualización de datos

data(uschange) autoplot(uschange[, 1:2]) autoplot(uschange[, c(1, 4)])

Plot de serie de datos

ts.plot(uschange[, 1:2], xlab = “Tiempo”, col = c(1, 2))

Búsqueda de parámetros

a <- VARselect(uschange[, 1:2], lag.max = 15, type = “trend”) a$selection

Creación y comparación de modelos

modelo1 <- VAR(uschange[, 1:2], p = 1, type = c(“trend”)) modelo3 <- VAR(uschange[, 1:2], p = 3, type = c(“trend”))

Resúmenes de modelos

summary(modelo1) summary(modelo3)

Comparación de AIC

aic1 <- summary(modelo1)$logLik aic2 <- summary(modelo3)$logLik

cat(“AIC modelo 1:”, aic1, “”) cat(“AIC modelo 3:”, aic2, “”)

summary(cars)

##      speed           dist       
##  Min.   : 4.0   Min.   :  2.00  
##  1st Qu.:12.0   1st Qu.: 26.00  
##  Median :15.0   Median : 36.00  
##  Mean   :15.4   Mean   : 42.98  
##  3rd Qu.:19.0   3rd Qu.: 56.00  
##  Max.   :25.0   Max.   :120.00

Laboratorio 6