powered by
getPredictiveIntervals(parMatrix, model, thin = 1000, quantiles = c(0.025, 0.975), error = NULL)
getPredictiveDistribution
getCredibleIntervals