defined n in VB EM; change log(gamma) to lgamma in KL