The expectation-maximization (EM) algorithm is an iterative method for finding maximum- likelihood estimates of parameters in statistical models with unobserved latent variables. Along with Markov Chain Monte Carlo (MCMC) it is one of the two computational workhorses that provided much impetus for statistics in entering its modern “computer-intensive” phase. Much is known about the EM...