Linear Pooling of Sample Covariance Matrices

We consider the problem of estimating high-dimensional covariance matrices of K-populations or classes in the setting where the sample sizes are comparable to the data dimension. We propose estimating each class cova Hance matrix as a distinct linear combination of all class sample covariance matrices. This approach is shown to reduce the estimation error when the sample sizes are limited, and the true class covariance matrices share a somewhat similar structure. We develop an effective method for estimating the coefficients in the linear combination that minimize the mean squared error under the general assumption that the samples are drawn from (unspecified) elliptically symmetric distributions possessing finite fourth-order moments. To this end, we utilize the spatial sign covariance matrix, which we show (under rather general conditions) to be an asymptotically unbiased estimator of the normalized covariance matrix as the dimension grows to infinity. We also show how the proposed method can be used in choosing the regularization parameters for multiple target matrices in a single class covariance matrix estimation problem. We assess the proposed method via numerical simulation studies including an application in global minimum variance portfolio optimization using real stock data.

Automated summary

This paper considers the problem of estimating high-dimensional covariance matrices of K-populations or classes when the sample sizes are comparable to the data dimension. It proposes a method to estimate each class covariance matrix as a linear combination of all class sample covariance matrices, which reduces the estimation error when the sample sizes are limited and the true class covariance matrices share a similar structure. The paper develops an effective method for estimating the coefficients in the linear combination and shows how the proposed method can be used for regularization parameter selection in a single class covariance matrix estimation problem. The proposed method is evaluated through numerical simulation studies and an application in global minimum variance portfolio optimization using real stock data.