Newton's method for beta-binomial maximum 
likelihood estimation