Internal function that samples z's from a categorical distribution

rmultinom1(prob, randu)

Arguments

prob

A numeric matrix.

randu

A numeric vector.