Monte Carlo EM based on Rao-Teh sampling seems to work