Monte-Carlo EM for Poisson Log-Normal model
Poisson Log-Normal model [Aitchison and Ho, 1989] is an incomplete data model for which the maximum likelihood estimator is not available via EM algorithm since the conditional distribution of the latent variable given the observed one is intractable. Efficient variational schemes have been proposed in the past few years. Even though they are computationally fast, they lack theoretical garanties and do not provide any confidence region.




