

def resample_psi(self, x):
    for k in xrange(self.K):
        (mu_cond, sigma_cond) = self.conditional_psi(x, sigma_cond)
        self.psi[k] = np.random.normal(mu_cond, np.sqrt(sigma_cond))

dataset/ETHPy150Open HIPS/pgmult/pgmult/distributions.py PGLogisticNormalMultinomial.resample_psi/VarMisuse@36/45 `k`->`sigma_cond`