def E_sumlogLam(self):
    return np.sum(digamma(self.a) - np.log(self.b))
\sum_d E[ \log \lambda_d ]
dataset/ETHPy150Open daeilkim/refinery/refinery/bnpy/bnpy-dev/bnpy/distr/GaussGammaDistr.py/GaussGammaDistr.E_sumlogLam