def resample_omega(self, x):
        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
Resample omega from its conditional Polya-gamma distribution
dataset/ETHPy150Open HIPS/pgmult/pgmult/distributions.py/PGLogisticNormalMultinomial.resample_omega