Internal function that samples z's from a categorical distribution
rmultinom1(prob, randu)
A numeric matrix.
A numeric vector.