[Numpy-discussion] Efficient way to draw multinomial distribution random samples

I realized that neither numpy.random.multinomial nor rng.multinomial has the ability to draw from different multinomial distributions at the same time like what MATLAB mnrnd() does here: https://www.mathworks.com/help/stats/mnrnd.html Also, I have asked this question on StackOverFlow: https://s

[Numpy-discussion] Multinomial random sampling

I wonder if there is an efficient way to draw multinomial random samples. For univariate sampling, we can do that for different parameters and desired shape, but neither of rng nor multinomial function in numpy.random is able to achieve that. To be specific, can we draw multinomial samples for diff