An internal function to perform the E step and the M step of the estimation algorithm.
This is not to be called by the user.
Cinzia Viroli