Statistical guarantees for the EM algorithm: from population to sample-based analysis

From MaRDI portal
Publication:524451

DOI10.1214/16-AOS1435zbMATH Open1367.62052arXiv1408.2156OpenAlexW2962737134MaRDI QIDQ524451FDOQ524451


Authors: Sivaraman Balakrishnan, Martin J. Wainwright, Bin Yu Edit this on Wikidata


Publication date: 2 May 2017

Published in: The Annals of Statistics (Search for Journal in Brave)

Abstract: We develop a general framework for proving rigorous guarantees on the performance of the EM algorithm and a variant known as gradient EM. Our analysis is divided into two parts: a treatment of these algorithms at the population level (in the limit of infinite data), followed by results that apply to updates based on a finite set of samples. First, we characterize the domain of attraction of any global maximizer of the population likelihood. This characterization is based on a novel view of the EM updates as a perturbed form of likelihood ascent, or in parallel, of the gradient EM updates as a perturbed form of standard gradient ascent. Leveraging this characterization, we then provide non-asymptotic guarantees on the EM and gradient EM algorithms when applied to a finite set of samples. We develop consequences of our general theory for three canonical examples of incomplete-data problems: mixture of Gaussians, mixture of regressions, and linear regression with covariates missing completely at random. In each case, our theory guarantees that with a suitable initialization, a relatively small number of EM (or gradient EM) steps will yield (with high probability) an estimate that is within statistical error of the MLE. We provide simulations to confirm this theoretically predicted behavior.


Full work available at URL: https://arxiv.org/abs/1408.2156




Recommendations





Cited In (67)





This page was built for publication: Statistical guarantees for the EM algorithm: from population to sample-based analysis

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q524451)