# G-StMAR model with intercept parametrization
params42gs <- c(0.04, 1.34, -0.59, 0.54, -0.36, 0.01, 0.06, 1.28, -0.36,
0.2, -0.15, 0.04, 0.19, 9.75)
gstmar42 <- GSMAR(data=M10Y1Y, p=4, M=c(1, 1), params=params42gs,
model="G-StMAR")
summary(gstmar42)
# Swap to mean parametrization
gstmar42 <- swap_parametrization(gstmar42)
summary(gstmar42)
Run the code above in your browser using DataLab