begin implementing dirichlet-multinomial version of the gaussian dpmm