Principal Regression for High-Dimensional Covariance Matrices

• Activity Number: 436 - Network Inference for Omics and Imaging Data
• Type: Topic-Contributed
• Date/Time: Thursday, August 12, 2021 : 4:00 PM to 5:50 PM
• Sponsor: Section on Statistical Learning and Data Science
• Author(s): Yi Zhao* and Brian Caffo and Xi Luo
• Companies: Indiana University and Johns Hopkins University and The University of Texas Health Science Center at Houston
• Keywords: Covariance matrix estimation; Generalized linear regression; Heteroscedasticity; Shrinkage estimator

Abstract

This paper presents an approach to perform generalized linear regression with multiple high dimensional covariance matrices as the outcome. Model parameters are proposed to be estimated by maximizing a pseudo-likelihood. When the data are high dimensional, the normal likelihood function is ill-posed as the sample covariance matrix is rank-deficient. A well-conditioned linear shrinkage estimator of the covariance matrix is introduced. With multiple covariance matrices, the shrinkage coefficients are proposed to be common across matrices. Theoretical studies demonstrate that the proposed covariance matrix estimator is optimal achieving the uniformly minimum quadratic loss asymptotically among all linear combinations of the identity matrix and the sample covariance matrix. Under regularity conditions, the proposed estimator of the model parameters is consistent. The superior performance of the proposed approach over existing methods is illustrated through simulation studies. Implemented to a resting-state fMRI study acquired from ADNI, the proposed approach identified a brain network within which functional connectivity is significantly associated with APOE-?4.