Computing the Dirichlet-Multinomial Log-Likelihood Function