DirMnom function

Dirichlet-multinomial (multivariate Polya) distribution