

def resample_omega(self, x):
    '\n        Resample omega from its conditional Polya-gamma distribution\n        :return:\n        '
    assert (x.ndim == 2)
    N = x.sum()
    ppg.pgdrawvpar(self.ppgs, (N * np.ones(self.K)), self.rho, self.omega)

dataset/ETHPy150Open HIPS/pgmult/pgmult/distributions.py PGLogisticNormalMultinomial.resample_omega/original