Sampling from dirichlet multinomial
Sampling from dirichlet multinomial