Maximum likelihood estimation for linear Gaussian covariance models

From MaRDI portal
Publication:5364909

DOI10.1111/RSSB.12217zbMATH Open1373.62267arXiv1408.5604OpenAlexW2963054928MaRDI QIDQ5364909FDOQ5364909


Authors: Piotr Zwiernik, Caroline Uhler, Donald Richards Edit this on Wikidata


Publication date: 29 September 2017

Published in: Journal of the Royal Statistical Society Series B: Statistical Methodology (Search for Journal in Brave)

Abstract: We study parameter estimation in linear Gaussian covariance models, which are p-dimensional Gaussian models with linear constraints on the covariance matrix. Maximum likelihood estimation for this class of models leads to a non-convex optimization problem which typically has many local maxima. Using recent results on the asymptotic distribution of extreme eigenvalues of the Wishart distribution, we provide sufficient conditions for any hill-climbing method to converge to the global maximum. Although we are primarily interested in the case in which n>!!>p, the proofs of our results utilize large-sample asymptotic theory under the scheme n/pogamma>1. Remarkably, our numerical simulations indicate that our results remain valid for p as small as 2. An important consequence of this analysis is that for sample sizes nsimeq14p, maximum likelihood estimation for linear Gaussian covariance models behaves as if it were a convex optimization problem.


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




Recommendations





Cited In (27)





This page was built for publication: Maximum likelihood estimation for linear Gaussian covariance models

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