

def Elogq_Lam(self):
    '\n    '
    logNormC = 0
    logPDF = 0
    for k in xrange(self.K):
        logNormC += ((gammaln((self.comp[k].lamA + self.comp[k].lamB)) - gammaln(self.comp[k].lamA)) - gammaln(self.comp[k].lamB))
        logPDF += (((self.comp[k].lamA - 1.0) * self.comp[k].ElogLamA) + ((self.comp[k].lamB - 1.0) * self.comp[k].ElogLamB))
    return ((- 1.0) * (logNormC + logPDF))

dataset/ETHPy150Open daeilkim/refinery/refinery/bnpy/bnpy-dev/bnpy/obsmodel/BernRelObsModel.py BernRelObsModel.Elogq_Lam/original