f5fc04d49bfd446086ff604cb9cad4f22fed02eb,pyhawkes/internals/network.py,GibbsSBM,resample_c,#GibbsSBM#Any#Any#,226

Before Change

                c_temp[k] = ck

                // p(A[k,k"] | c)
                lp[ck] += Bernoulli(self.p[np.ix_([ck], c_temp)])\
                                .log_probability(A[k,:]).sum()
                // p(W[k,k"] | c)
                lp[ck] += A[k,k] * Gamma(self.kappa, self.v[np.ix_([ck], c_temp)])\

After Change

                                .log_probability(A[k,:]).sum()

                // p(A[k",k] | c)
                lp[ck] += Bernoulli(self.p[c_temp, ck])\
                                .log_probability(A[:,k]).sum()

                // p(W[k,k"] | c)
                lp[ck] += (A[k,:] * Gamma(self.kappa, self.v[ck, c_temp])\
                                .log_probability(W[k,:])).sum()

Instances