HDIbetabinom function

Get highest density interval of beta-binomial