Maximum Likelihood Estimation of Beta-Binomial Distribution using R