// Function estimates parameters of beta-binomial distribution (using count
// data) and stores them in the vector of coefficients as
// {[5] alpha (p), [6] beta (q), [7] log(std::beta)}

// REFS:
//  1) https://en.wikipedia.org/wiki/Beta-binomial_distribution#Method_of_moments

