class Multinomial(object):

    def __init__(self, params):
        self._params = params

def pmf(self, counts):
    if not len(counts) == len(self._params):
        raise ValueError('Dimensionality of count vector is incorrect')
prob = 1.0
multinomial pmf in python scipy/numpy