# \donttest{
simdat <- sim_mvgam()
mod <- mvgam(y ~ s(season, bs = 'cc',
k = 6),
trend_model = AR(),
use_lv = TRUE,
n_lv = 2,
data = simdat$data_train,
burnin = 300,
samples = 300,
chains = 2,
silent = 2)
lvcors <- lv_correlations(mod)
names(lvcors)
lapply(lvcors, class)
# }
Run the code above in your browser using DataLab