Alternative implementations of Monte Carlo EM algorithms for likelihood inferences

Abstract

Abstract Two methods of computing Monte Carlo estimators of variance components using restricted maximum likelihood via the expectation-maximisation algorithm are reviewed. A third approach is suggested and the performance of the methods is compared using simulated data.

Keywords

• restricted maximum likelihood
• Markov chain Monte Carlo
• EM algorithm
• Monte Carlo variance
• variance components