Dirichlet-multinomial distribution