pm.sample gives you a big array of samples; so pass it to np.percentile on it to grab the appropriate percentiles. As this is a random sample, this won’t be exact. You can sample 10 chains, and obtain percentiles for each chain, to establish the sampling uncertainty of the percentile (or you can use lots of samples).
1 Like