Back To Index Previous Article Next Article Full Text

Statistica Sinica 31 (2021), 173-196

ESTIMATING LARGE PRECISION MATRICES
VIA MODIFIED CHOLESKY DECOMPOSITION

Kyoungjae Lee and Jaeyong Lee

Inha University and Seoul National University

Abstract: We introduce a k-banded Cholesky prior for estimating high-dimensional bandable precision matrices using a modified Cholesky decomposition. The bandable assumption is imposed on the Cholesky factor of the decomposition. We obtain the P-loss convergence rate under the spectral norm and the matrix 𝓁-norm, as well as the minimax lower bounds. Because the P-loss convergence rate is stronger than the posterior convergence rate, the rates obtained are also posterior convergence rates. Furthermore, when the true precision matrix is a k0-banded matrix, for some finite k0 , we obtain the minimax rate. The established convergence rates for bandable precision matrices are slightly slower than the minimax lower bounds, but are the fastest of the existing Bayesian approaches. Simulation results show that the performance of the proposed method is better than or comparable to that of competitive estimators.

Key words and phrases: Modified Cholesky decomposition, P-loss convergence rate, precision matrix.

Back To Index Previous Article Next Article Full Text