Since probabilities can become extremely small,
it is necessary to remain in the log-domain.
Staying in the log-domain is a problem
when summations are required.
Let
.
The summation

which appears in the first step of the E-M algorithm
should be implemented as

where
.

Baggenstoss
2017-05-19