Internal function that samples z's from a multinomial distribution
rmultinom2(prob, n, randu, nmaxclust)
A numeric vector.
An integer.