'Draw random samples from a Multinomial distribution.
Parameters
%(_doc_default_callparams)s
size : integer or iterable of integers, optional
Number of samples to draw (default 1).
%(_doc_random_state)s
Returns
rvs : ndarray or scalar
Random variates of shape (`size`, `len(p)`)
Notes
%(_doc_callparams_note)s'

def rvs(self, n, p, size=None, random_state=None):

    (n, p, npcond) = self._process_parameters(n, p)
    random_state = self._get_random_state(random_state)
    return random_state.multinomial(n, p, size)
