Sparsistency and rates of convergence in large covariance matrix estimation

Clifford Lam, Jianqing Fan

Research output: Contribution to journalArticle

259 Scopus citations

Abstract

This paper studies the sparsistency and rates of convergence for estimating sparse covariance and precision matrices based on penalized likelihood with nonconvex penalty functions. Here, sparsistency refers to the property that all parameters that are zero are actually estimated as zero with probability tending to one. Depending on the case of applications, sparsity priori may occur on the covariance matrix, its inverse or its Cholesky decomposition. We study these three sparsity exploration problems under a unified framework with a general penalty function. We show that the rates of convergence for these problems under the Frobenius norm are of order (sn logpn/n) 1/2, where sn is the number of nonzero elements, p n is the size of the covariance matrix and n is the sample size. This explicitly spells out the contribution of high-dimensionality is merely of a logarithmic factor. The conditions on the rate with which the tuning parameter λn goes to 0 have been made explicit and compared under different penalties. As a result, for the L1 -penalty, to guarantee the sparsistency and optimal rate of convergence, the number of nonzero elements should be small: s'n = O(pn) at most, among O(p 2n) parameters, for estimating sparse covariance or correlation matrix, sparse precision or inverse correlation matrix or sparse Cholesky factor, where s'n is the number of the nonzero elements on the off-diagonal entries. On the other hand, using the SCAD or hard-thresholding penalty functions, there is no such a restriction.

Original languageEnglish (US)
Pages (from-to)4254-4278
Number of pages25
JournalAnnals of Statistics
Volume37
Issue number6 B
DOIs
StatePublished - Dec 2009

All Science Journal Classification (ASJC) codes

  • Statistics and Probability
  • Statistics, Probability and Uncertainty

Keywords

  • Asymptotic normality
  • Consistency
  • Covariance matrix
  • High-dimensionality
  • Nonconcave penalized likelihood
  • Sparsistency

Fingerprint Dive into the research topics of 'Sparsistency and rates of convergence in large covariance matrix estimation'. Together they form a unique fingerprint.

  • Cite this