dirmult function

Dirichlet-multinomial distribution