bbdml function

Maximum Likelihood for the Beta-binomial Distribution