High-Dimensional Two-Sample Covariance Testing with Null-Preserving Transformations
Abstract
Testing equality of two high-dimensional covariance matrices is challenging when many entries differ only slightly. Dependence among sample covariance entries can also affect the finite-sample size and power of tests that aggregate their differences. We propose a studentized ℓ2-type statistic that averages squared, marginally standardized differences between corresponding entries of the two sample covariance matrices. Before constructing the statistic, we apply the same nonsingular linear transformation to both samples. In the transformed coordinates, the covariance-equality hypothesis is unchanged, whereas the standardized differences and correlations among their estimators generally change. The transformation can therefore incorporate structural or scientific information without reducing dimension. We approximate the null distribution using a Gaussian multiplier bootstrap that uses coordinatewise studentization and avoids forming or inverting the full covariance matrix of the vectorized sample covariance entries. For deterministic transformations, we derive nonasymptotic Gaussian and bootstrap approximation bounds and establish asymptotic size validity and power consistency. We also show that replacing a population transformation by a same-sample estimator leaves the statistic and bootstrap critical value asymptotically unchanged in relative terms under an operator-norm convergence condition. Simulations show that suitable transformations improve size accuracy and power under weak dependence. An analysis of breast cancer gene-expression data illustrates the method.
Related
- On Eigenvector Statistics of Large-Dimensional Sample Covariance Matrices with Applications
- An F-Type Normal-Reference Calibration for High-Dimensional Covariance Testing
- High-dimensional two-sample testing: From optimal transport to optimal derangements
- Linear transformations of tropical matrices preserving individual values of the cyclicity index
- Stochastic Inversion of Multivariate Uniform-Distribution-Preserving Transformations
BibTeX
@misc{shu2026highdimensional,
title = {{High-Dimensional Two-Sample Covariance Testing with Null-Preserving Transformations}},
author = {Haozhen Shu and Tianming Zhu},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.11249},
eprint = {2610.11249},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.11249},
}
RIS
TY - UNPB
TI - High-Dimensional Two-Sample Covariance Testing with Null-Preserving Transformations
AU - Haozhen Shu
AU - Tianming Zhu
PY - 2026
DA - 2026/10/08
PB - arXiv
DO - 10.48550/arxiv.2610.11249
UR - https://arxiv.org/abs/2610.11249
AB - Testing equality of two high-dimensional covariance matrices is challenging when many entries differ only slightly. Dependence among sample covariance entries can also affect the finite-sample size and power of tests that aggregate their differences. We propose a studentized ℓ_2-type statistic that averages squared, marginally standardized differences between corresponding entries of the two sample covariance matrices. Before constructing the statistic, we apply the same nonsingular linear transformation to both samples. In the transformed coordinates, the covariance-equality hypothesis is unchanged, whereas the standardized differences and correlations among their estimators generally change. The transformation can therefore incorporate structural or scientific information without reducing dimension. We approximate the null distribution using a Gaussian multiplier bootstrap that uses coordinatewise studentization and avoids forming or inverting the full covariance matrix of the vectorized sample covariance entries. For deterministic transformations, we derive nonasymptotic Gaussian and bootstrap approximation bounds and establish asymptotic size validity and power consistency. We also show that replacing a population transformation by a same-sample estimator leaves the statistic and bootstrap critical value asymptotically unchanged in relative terms under an operator-norm convergence condition. Simulations show that suitable transformations improve size accuracy and power under weak dependence. An analysis of breast cancer gene-expression data illustrates the method.
ER -