Limiting spectral distribution of renormalized separable sample covariance matrices when

We are concerned with the behavior of the eigenvalues of renormalized sample covariance matrices of the form C_n=\sqrt{\frac{n}{p}}\left(\frac{1}{n}A_{p}^{1/2}X_{n}B_{n}X_{n}^{*}A_{p}^{1/2}-\frac{1}{n}\tr(B_{n})A_{p}\right) as and , where is a matrix with i.i.d. real or complex valued entries satisfying , and having finite fourth moment. is a square-root of the nonnegative definite Hermitian matrix , and is an nonnegative definite Hermitian matrix. We show that the empirical spectral distribution (ESD) of converges a.s. to a nonrandom limiting distribution under some assumptions. The probability density function of the LSD of is derived and it is shown that it depends on the LSD of and the limiting value of . We propose a computational algorithm for evaluating this limiting density when the LSD of is a mixture of point masses. In addition, when the entries of are sub-Gaussian, we derive the limiting empirical distribution of where is the sample covariance matrix and denotes the -th largest eigenvalue, when is a finite mixture of point masses. These results are utilized to propose a test for the covariance structure of the data where the null hypothesis is that the joint covariance matrix is of the form for denoting the Kronecker product, as well as and the first two spectral moments of are specified. The performance of this test is illustrated through a simulation study.
View on arXiv