Temporary workaround for an arithmetic bug
Incomplete LDA/VB Smoothing Work

Conversion to appropriate estimation of Dirichlet parameters using
digamma’s instead of hacking it with means.