

def resample_omega(self, x):
    '\n        Resample omega from its conditional Polya-gamma distribution\n        :return:\n        '
    assert (x.ndim == 2)
    N = N_vec(x)
    NN = N.sum(0).astype(np.float)
    ppg.pgdrawvpar(self.ppgs, NN, N.psi, self.omega)

dataset/ETHPy150Open HIPS/pgmult/pgmult/distributions.py PGMultinomial.resample_omega/VarMisuse@56/95 `self`->`N`