polyaMLE function

Obtain MLE of marginal Dirichlet-multinomial likelihood