# NOT RUN {
require(MASS)
sig <- matrix(c(1, 0.5, 0.5, 1), ncol=2)
set.seed(123)
Y <- mvrnorm(50, mu=c(-2, 2), Sigma=sig)
dat <- data.frame(id = rep(1:5, each=10),
y1 = Y[, 1],
y2 = Y[, 2])
f <- bayesGAM(cbind(y1, y2) ~ 1, random = ~factor(id),
data=dat,
a = normal(c(0, 5)),
chains = 1, iter = 500)
mvcorrplot(f)
# }
Run the code above in your browser using DataLab