Estimation of Large Precision Matrices Through Block Penalization
Abstract
This paper focuses on exploring the sparsity of the inverse covariance matrix , or the precision matrix. We form blocks of parameters based on each off-diagonal band of the Cholesky factor from its modified Cholesky decomposition, and penalize each block of parameters using the -norm instead of individual elements. We develop a one-step estimator, and prove an oracle property which consists of a notion of block sign-consistency and asymptotic normality. In particular, provided the initial estimator of the Cholesky factor is good enough and the true Cholesky has finite number of non-zero off-diagonal bands, oracle property holds for the one-step estimator even if , and can even be as large as , where the data has mean zero and tail probability , , and is the number of variables. We also prove an operator norm convergence result, showing the cost of dimensionality is just . The advantage of this method over banding by Bickel and Levina (2008) or nested LASSO by Levina \emph{et al.} (2007) is that it allows for elimination of weaker signals that precede stronger ones in the Cholesky factor. A method for obtaining an initial estimator for the Cholesky factor is discussed, and a gradient projection algorithm is developed for calculating the one-step estimate. Simulation results are in favor of the newly proposed method and a set of real data is analyzed using the new procedure and the banding method.
Keywords
Cite
@article{arxiv.0805.3798,
title = {Estimation of Large Precision Matrices Through Block Penalization},
author = {Clifford Lam},
journal= {arXiv preprint arXiv:0805.3798},
year = {2008}
}
Comments
42 pages article