How can I determine the confidence/credibility intervals for the posterior estimates of a multi-parameter model? I can get the confidence interval for each parameter separately. (Currently using bayestestR, but I don’t mind using something else) library(dplyr) library(ggplot2) library(bayestestR) N <- 10000 #Posterior samples (random example) p1 <- rnorm(N) p2 <- p1 + rnorm(N) df_post <-… Read More Two-parameter Bayesian credibility Interval