# NOT RUN {
# data
Y <- subset(tas, gender == "M")[,-ncol(tas)]
# fit model
fit <- estimate(Y)
# posterior predictions
pred <- posterior_predict(fit, iter = 500,
summary = FALSE)
# prediction error
error <- mse(pred)
# plot
plot(error)
# }
Run the code above in your browser using DataLab