EMFit function

Compute the Maximum-A-Posteriori estimate with the EM algorithm