Alternative implementations of Monte Carlo EM algorithms for likelihood inferences

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.