def log_likelihood(self, x):
        ll = 0
        ll += gammaln((x+1).sum()) - gammaln(x-1).sum()
        ll += (x * np.log(self.pi)).sum()
        return ll
dataset/ETHPy150Open HIPS/pgmult/pgmult/distributions.py/PGMultinomial.log_likelihood/operator `+` -> `-`