def resample_omega(self, x):
        """
        Resample omega from its conditional Polya-gamma distribution
        :return:
        """
        assert x.ndim == 2
        N = x.sum()

        #  Sum the N's (i.e. the b's in the denominator)
        ppg.pgdrawvpar(self.ppgs, N - np.ones(self.K), self.rho, self.omega)




### Competing models
dataset/ETHPy150Open HIPS/pgmult/pgmult/distributions.py/PGLogisticNormalMultinomial.resample_omega/operator `*` -> `-`