powered by
getPredictiveDistribution(parMatrix, model, thin = 1000)
getPredictiveIntervals
getCredibleIntervals