TY - JOUR

T1 - Regularization for high-dimensional covariance matrix

AU - Cui , Xiangzhao

AU - Li, Chun

AU - Zhao , Jine

AU - Zeng , Li

AU - Zhang , Defei

AU - Pan, Jianxin

PY - 2016/4/12

Y1 - 2016/4/12

N2 - In many applications, high-dimensional problem may occur often for various reasons, for example, when the number of variables under consideration is much bigger than the sample size, i.e., p >> n. For highdimensional data, the underlying structures of certain covariance matrix estimates are usually blurred due to substantial random noises, which is an obstacle to draw statistical inferences. In this paper, we propose a method to identify the underlying covariance structure by regularizing a given/estimated covariance matrix
so that the noises can be ltered. By choosing an optimal structure from a class of candidate structures for the covariance matrix, the regularization is made in terms of minimizing Frobenius-norm discrepancy. The candidate class considered here includes the structures of order-1 moving average, compound symmetry, order-1 autoregressive and order-1 autoregressive moving average. Very intensive simulation studies are conducted to assess the performance of the proposed regularization method for very high-dimensional covariance problem. The simulation studies also show that the sample covariance matrix, although performs very badly in covariance estimation for high-dimensional data, can be used to correctly identify the underlying structure of the covariance matrix. The approach is also applied to real data analysis, which shows that the proposed regularization method works well in practice.

AB - In many applications, high-dimensional problem may occur often for various reasons, for example, when the number of variables under consideration is much bigger than the sample size, i.e., p >> n. For highdimensional data, the underlying structures of certain covariance matrix estimates are usually blurred due to substantial random noises, which is an obstacle to draw statistical inferences. In this paper, we propose a method to identify the underlying covariance structure by regularizing a given/estimated covariance matrix
so that the noises can be ltered. By choosing an optimal structure from a class of candidate structures for the covariance matrix, the regularization is made in terms of minimizing Frobenius-norm discrepancy. The candidate class considered here includes the structures of order-1 moving average, compound symmetry, order-1 autoregressive and order-1 autoregressive moving average. Very intensive simulation studies are conducted to assess the performance of the proposed regularization method for very high-dimensional covariance problem. The simulation studies also show that the sample covariance matrix, although performs very badly in covariance estimation for high-dimensional data, can be used to correctly identify the underlying structure of the covariance matrix. The approach is also applied to real data analysis, which shows that the proposed regularization method works well in practice.

U2 - 10.1515/spma-2016-0018

DO - 10.1515/spma-2016-0018

M3 - Article

JO - Special Matrices

JF - Special Matrices

SN - 2300-7451

ER -