msmle.py 文件源码

def _W_nj(self):
        """The sample sorted sample weights in all states.

        FOR INTERNAL USE ONLY!
        """
        m = self.nstates_sampled
        logQ_nj = self._f[newaxis, :] - self._u_nj
        logNorm_n = logsumexp(logQ_nj[:, :m], 1, self.PIsdiag[newaxis, :])
        _W_nj = exp(logQ_nj - logNorm_n[:, newaxis])
        return _W_nj