# NOT RUN {
# GMAR model with intercept parametrization
params12 <- c(0.18, 0.93, 0.01, 0.86, 0.68, 0.02, 0.88)
gmar12 <- GSMAR(data=logVIX, p=1, M=2, params=params12, model="GMAR")
gmar12
# Swap to mean parametrization
gmar12 <- swap_parametrization(gmar12)
gmar12
# }
Run the code above in your browser using DataLab