REVIEW 2 major objections 5 minor 71 references
Separability of a matrix covariance can be tested by whitening with the separable MLE and checking sphericity, with an angular version that resists heavy tails.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · grok-4.5
2026-07-10 08:15 UTC pith:HDUNP2SB
load-bearing objection Solid high-d separability test with clean finite-sample level and a useful angular fix; consistency is real but conditional on open pseudo-MLE rates the authors flag themselves. the 2 major comments →
Testing Covariance Separability in High Dimensions
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Within the elliptical family, covariance separability is equivalent to sphericity after whitening by a unique matrix-affine-equivariant separable estimator (the matrix-normal MLE). The resulting Monte Carlo test based on the unified John/Nagao/Ledoit–Wolf statistic therefore has exact finite-sample level under the null, and is consistent in the high-dimensional regime N proportional to pq against dense alternatives measured in Frobenius distance. An angular refinement that radially normalizes the whitened observations inherits the same level property and consistency under slightly stronger conditioning assumptions, while empirically reducing sensitivity to radial misspecification.
What carries the argument
Matrix whitening by the flip-flop MLE of the separable factors, after which John's, Nagao's and Ledoit–Wolf's sphericity statistics coincide (up to a constant absorbed by Monte Carlo); the angular version replaces the sample covariance of the whitened data by the spatial-sign covariance obtained after unit-sphere projection.
Load-bearing premise
Consistency rests on the unproven claim that the separable pseudo-MLE still converges, at a high-dimensional rate, when the true covariance is not separable.
What would settle it
Construct a dense non-separable family (fixed positive Frobenius gap after scaling by 1/√(pq)) for which the separable MLE fails to converge in Frobenius or operator norm when N ~ pq, or show that the angular test's empirical level collapses under a non-elliptical law while the Gaussian Monte Carlo quantiles stay fixed.
If this is right
- Separability can be checked at roughly the cost of one separable fit, without ever storing the full pq-by-pq covariance.
- When N is proportional to pq, dense non-separability is asymptotically detectable after MMLE whitening.
- The angular test can be used under unknown elliptical tails with only a Gaussian Monte Carlo table, retaining nearly full power.
- Pre-tabulated quantiles for each (p,q,N) make repeated testing and software packaging immediate.
- Downstream matrix-variate procedures that rely on separability gain a reliable pre-test that respects the same computational budget.
Where Pith is reading between the lines
- The same whitening-plus-sphericity reduction is likely usable for other Kronecker-type structural hypotheses (e.g., multi-way separability or core-shrinkage models) once an equivariant estimator under the null is available.
- Closing the open high-dimensional rate for the pseudo-MLE under alternatives would immediately upgrade both consistency theorems from conditional to fully rigorous.
- Because the angular statistic discards only radial scale, it may remain approximately pivotal for mild departures from ellipticity that still preserve spherical directions after whitening—an empirical regime worth mapping.
- The eigenvalue-dispersion view of the statistic suggests natural sparse-alternative competitors (largest eigenvalue after whitening) that the paper already contrasts and that could be hybridized with the angular projection.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes high-dimensional tests of covariance separability for matrix-variate data by whitening with the matrix-normal MMLE and testing sphericity of the whitened sample. After whitening, John’s, Nagao’s, and Ledoit–Wolf statistics coincide (Lemma 1); Monte Carlo calibration under a matrix-affine-equivariant estimator yields exact finite-sample level under the null (Theorem 1, Corollaries 1–2, Proposition 1). Consistency against dense Frobenius alternatives is claimed under Assumptions (A1)–(A5) for the elliptical statistic (Theorem 2) and under (B1)–(B5) for an angular (spatial-sign) version that projects whitened observations onto the unit sphere (Theorem 3, Proposition 3). Simulations and an acoustic phonetic application support power and improved robustness of the angular test under heavy-tailed elliptical laws with Gaussian calibration.
Significance. The problem is timely: high-dimensional matrix data make unstructured covariance estimation infeasible, so a computationally light, finite-sample valid test of separability is useful. The finite-sample level theory via matrix affine equivariance is clean and standard. Unifying classical sphericity statistics after MMLE whitening, the efficient inner-product form of the statistic, and the angular refinement for distributional robustness are genuine contributions. An R package with pre-tabulated quantiles is a practical strength. The consistency theorems would be a solid theoretical contribution if the load-bearing pseudo-MLE rates under alternatives were established or clearly demoted; as written they remain conditional.
major comments (2)
- [§3.2 Assumptions (A4)–(A5); §5 (B4)–(B5); Theorems 2–3; Discussion §7; Appendix A.1] Theorems 2 and 3 (and the separation arguments that follow Assumptions (A5)/(B5)) treat as given the high-dimensional rates (A4) 1/√(pq)‖Σ̂1⊗̃Σ̂2−Σ1⊗̃Σ2‖_F=o_p(1) and the stronger operator-norm (B4). The manuscript states these rates under non-separable alternatives are open (Discussion §7; comparison to Franks et al. 2026, which covers only the separable case) and supports them only by numerical box-plots in Appendix A.1. Without (A4)/(B4) the inheritance of dense separation after whitening can fail, so the consistency claims are conditional on an unproved estimation result that is load-bearing for the central theoretical contribution. Either prove the rates, replace them by weaker verifiable conditions, or restate Theorems 2–3 as conditional on (A4)/(B4) and move the unconditional claim to finite-sample level plus empirical power.
- [Abstract; §1; Theorems 2–3] The abstract and introduction state that the paper proves high-dimensional consistency under dense alternatives without flagging that the proofs rest on open rates for the pseudo-true separable factors. That overstates what is established. Soften the abstract/intro claims to match the conditional status of Theorems 2–3, or complete the rate analysis.
minor comments (5)
- [Abstract] Abstract is truncated mid-sentence (“of the prop”). Fix before resubmission.
- [§2 Preliminaries] Notation for the tensor product ˜b versus Kronecker product is introduced carefully but still dense; a short display equating vec((Σ1˜bΣ2)X)=(Σ2⊗Σ1)vec(X) early in §2 would help readers who work only with Kronecker notation.
- [§6.1.2; Figure 3] Figure 3 panels for df=5 and small m show residual size inflation for the angular test; the text acknowledges this but could quantify how large N must be before Gaussian calibration is reliable under matrix-t.
- [§4] Comparison to Sung and Hoff (2026) is useful; make explicit in the main text (not only §4) which of their statistics coincides with the elliptical test after the latest revision, to avoid reader confusion.
- [Proposition 1] Proposition 1’s sample-size condition N≥⌈p/q+q/p⌉+2 is cited from Drton et al. (2021); a one-sentence reminder that this is almost necessary would help practitioners.
Circularity Check
No circularity: level and consistency claims rest on equivariance, classical sphericity functionals after whitening, and external dense-alternative assumptions rather than self-referential definitions or load-bearing self-citations.
full rationale
The derivation chain is self-contained and non-circular. Under the null, finite-sample level (Corollary 2, Proposition 2) follows from matrix-affine equivariance of the MMLE (Proposition 1, citing the external uniqueness result of Drton et al. 2021) plus orthogonal invariance of the eigenvalue-based statistics (Theorem 1 and Corollary 1); Monte Carlo calibration is performed under a stated F0 and does not fit parameters to the observed data. Consistency (Theorems 2 and 3) is proved by showing that dense Frobenius departure (A5/B5) transfers to a positive separation of the whitened (or angular) statistic after the estimation error controlled by the open rates (A4/B4); those rates are explicitly flagged as unproved under alternatives and are only supported numerically (Appendix A.1, Discussion §7), so the theorems are conditional rather than circular. Self-citations (Masak et al. 2023 on separable components, Mayrhofer et al. 2024 on MMLE equivariance) supply background estimators and expansions used in simulations or auxiliary lemmas; they do not underwrite the target level or consistency statements. The concurrent Sung–Hoff overlap is disclosed as independent. No step reduces a claimed prediction or first-principles result to its own inputs by construction, fitted parameter, or self-citation chain.
Axiom & Free-Parameter Ledger
free parameters (3)
- M (Monte Carlo replicates) =
999
- Gaussian default F0 for MC under unknown law =
matrix-normal reference
- Sample-size multiplier m defining N=max(4,⌈m pq⌉) =
grid in {0.1,0.25,0.5,1,2,4}
axioms (7)
- domain assumption Model 1: X_n = Σ^{1/2} Z_n with Z_n i.i.d. absolutely continuous, orthogonally invariant, mean 0, Cov = I (elliptical/spherical generator).
- domain assumption MMLE exists uniquely a.s. for N ≥ ⌈p/q + q/p⌉ + 2 and is matrix affine equivariant.
- domain assumption (A1)/(B1) high-d regime p/√N → γ1, q/√N → γ2 ∈ (0,∞) so pq ≍ N.
- ad hoc to paper (A4)/(B4) separable (pseudo-)MLE converges to the KL-projection factors at o_p(1) Frobenius or operator rate under alternatives.
- domain assumption (A5)/(B5) dense alternative: 1/√(pq) ‖Σ − Σ1⊗̃Σ2‖_F ≥ Δ0 > 0 for large N.
- domain assumption Bounded operator norms and eigenvalues of Σ, Σ1, Σ2 away from 0 and ∞ uniformly in N ((A2)/(B2)).
- standard math Finite eighth moments of entries of Z_n for elliptical consistency ((A3)); milder absolute continuity/orthogonal invariance for angular ((B3)).
invented entities (2)
-
Matrix-whitened unified John/Nagao/Ledoit–Wolf separability statistic T_N
independent evidence
-
Angular separability statistic T_N^{(A)} (spatial-sign covariance after whitening)
independent evidence
read the original abstract
Separability is an important structural assumption often placed on the covariance when working with matrix-variate data, because it greatly simplifies both interpretation and computation of subsequent covariance-based statistical tasks. Yet testing the separability assumption is difficult in the high-dimensional regime. We propose to test separability by recasting the problem as a sphericity test after whitening the data using the separable maximum likelihood estimate of the covariance. The test is calibrated by Monte Carlo simulation, yielding finite-sample level control. Furthermore, we prove the test's high-dimensional consistency under dense alternatives. To reduce its reliance on distributional assumptions, we introduce an angular version of the test based on radial normalization after whitening. We demonstrate the practical utility, empirical power, and computational efficiency of the prop
Figures
Reference graph
Works this paper leans on
- [1]
-
[2]
Tyler and Daniel Vogel , keywords =
Alexander Dürre and David E. Tyler and Daniel Vogel , keywords =. On the eigenvalues of the spatial sign covariance matrix in more than two dimensions , journal =. 2016 , issn =. doi:https://doi.org/10.1016/j.spl.2016.01.009 , url =
-
[3]
Analysis of Time Series Structure: SSA and Related Techniques , author =. 2001 , publisher =
work page 2001
- [4]
- [5]
-
[6]
Gupta, A. K. and Nagar, D. K. , title =. 2018 , publisher =
work page 2018
- [7]
-
[8]
Van Loan, C. F. , title =. Journal of Computational and Applied Mathematics , year =
-
[9]
Van Loan, C. F. and Pitsianis, N. , title =. Linear algebra for large scale and real-time applications , year =
-
[10]
Golub, G. H. and Van Loan, C. F. , title =. 2013 , publisher =
work page 2013
-
[11]
Gupta, A. K. and Varga, T. , title =. Journal of the Italian Statistical Society , year =
-
[12]
Davison, A. C. and Hinkley, D. V. , title =. 1997 , series =
work page 1997
-
[13]
Li, W. and Yao, J. , title =. Journal of the Royal Statistical Society Series B: Statistical Methodology , year =
- [14]
-
[15]
Mayrhofer, M. and Radojičić, U. and Filzmoser, P. , title =. 2024 , eprint =
work page 2024
-
[16]
Mitchell, M. W. and Genton, M. G. and Gumpertz, M. L. , title =. Journal of Multivariate Analysis , year =
-
[17]
Aston, J. A. D. and Pigoli, D. and Tavakoli, S. , title =. The Annals of Statistics , year =
-
[18]
Journal of the American Statistical Association , volume=
Random surface covariance estimation by shifted partial tracing , author=. Journal of the American Statistical Association , volume=. 2023 , publisher=
work page 2023
- [19]
-
[20]
Lu, N. and Zimmerman, D. L. , title =. Statistics & Probability Letters , year =
-
[21]
Filipiak, K. and Klein, D. and Vojtkov. The properties of partial trace and block trace operators of partitioned matrices , journal =. 2018 , volume =
work page 2018
- [22]
-
[23]
Rosenthal, H. P. , title =. Israel Journal of Mathematics , year =
- [24]
-
[25]
Merlev. Rosenthal-type inequalities for the maximum of partial sums of stationary processes and examples , journal =. 2013 , volume =
work page 2013
-
[26]
Drton, M. and Kuriki, S. and Hoff, P. D. , title =. The Annals of Statistics , year =
-
[27]
Srivastava, M. S. and von Rosen, T. and von Rosen, D. , title =. Mathematical Methods of Statistics , year =
-
[28]
Franks, C. and Oliveira, R. and Ramachandran, A. and Walter, M. , title =. The Annals of Statistics , year =
-
[29]
Masak, T. and Sarkar, S. and Panaretos, V. M. , title =. Biometrika , year =
-
[30]
Testing Separability of High-Dimensional Covariance Matrices
Sung, B. and Hoff, P. D. , title =. arXiv preprint arXiv:2506.17463 , year =. 2506.17463 , archivePrefix =
work page internal anchor Pith review Pith/arXiv arXiv
-
[31]
Separable Geodesic Lagrangian Monte Carlo for Inference in 2-Way Covariance Models
Simonis, Q. and Wells, M. T. , title =. arXiv preprint arXiv:2501.04913 , year =. 2501.04913 , archivePrefix =
work page internal anchor Pith review Pith/arXiv arXiv
-
[32]
Dette, H. and Dierickx, G. and Kutta, T. , title =. Journal of Time Series Analysis , year =
- [33]
-
[34]
Soloveychik, I. and Trushin, D. , title =. Journal of Multivariate Analysis , year =
- [35]
-
[36]
Thompson, G. Z. and Maitra, R. and Meeker, W. Q. and Bastawros, A. F. , title =. Journal of Computational and Graphical Statistics , year =
-
[37]
Pigoli, D. and Hadjipantelis, P. Z. and Coleman, J. S. and Aston, J. A. D. , title =. Journal of the Royal Statistical Society Series C: Applied Statistics , year =
-
[38]
Pigoli, D. and Aston, J. A. D. and Dryden, I. L. and Secchi, P. , title =. Biometrika , year =
- [39]
-
[40]
Tang, R. and M. Pairwise curve synchronization for functional data , journal =. 2008 , volume =
work page 2008
- [41]
- [42]
- [43]
-
[44]
Gneiting, T. and Genton, M. G. and Guttorp, P. , title =. Statistical Methods for Spatio-Temporal Systems , editor =. 2007 , pages =
work page 2007
-
[45]
Rougier, J. , title =. arXiv preprint arXiv:1702.05599 , year =. 1702.05599 , archivePrefix =
work page internal anchor Pith review Pith/arXiv arXiv
- [46]
-
[47]
Constantinou, P. and Kokoszka, P. and Reimherr, M. , title =. Biometrika , year =
-
[48]
Wang, J.-L. and Chiou, J.-M. and M. Functional data analysis , journal =. 2016 , volume =
work page 2016
-
[49]
Dawid, A. P. , title =. Biometrika , year =
-
[50]
Mardia, K. V. and Goodall, C. R. , title =. Multivariate Environmental Statistics , editor =. 1993 , pages =
work page 1993
- [51]
-
[52]
Zhang, Y. and Shen, W. and Kong, D. , title =. Journal of the American Statistical Association , year =
-
[53]
Hoff, P. D. and McCormack, A. and Zhang, A. R. , title =. Journal of the Royal Statistical Society Series B: Statistical Methodology , year =
- [54]
- [55]
-
[56]
Filipiak, K. and Klein, D. and Roy, A. , title =. Journal of Multivariate Analysis , year =
-
[57]
Greenewald, K. and Zhou, S. and Hero, III, A. , title =. Journal of the Royal Statistical Society Series B: Statistical Methodology , year =
- [58]
-
[59]
Chen, S. X. and Zhang, L.-X. and Zhong, P.-S. , title =. Journal of the American Statistical Association , year =
-
[60]
Bai, Z. and Jiang, D. and Yao, J.-F. and Zheng, S. , title =. The Annals of Statistics , year =
-
[61]
Genton, M. G. , title =. Environmetrics , year =
-
[62]
Yin, Y. Q. and Bai, Z. D. and Krishnaiah, P. R. , title =. Probability Theory and Related Fields , year =
-
[63]
Eaton, M. L. , title =
- [64]
-
[65]
Journal of the Royal Statistical Society Series B: Statistical Methodology , volume=
Matrix variate regressions and envelope models , author=. Journal of the Royal Statistical Society Series B: Statistical Methodology , volume=. 2018 , publisher=
work page 2018
-
[66]
The Annals of Statistics , number =
Shuheng Zhou , title =. The Annals of Statistics , number =. 2014 , doi =
work page 2014
-
[67]
Journal of the American Statistical Association , volume=
Testing for nodal dependence in relational data matrices , author=. Journal of the American Statistical Association , volume=. 2015 , publisher=
work page 2015
-
[68]
The annals of applied statistics , volume=
Latent class modeling using matrix covariates with application to identifying early placebo responders based on EEG signals , author=. The annals of applied statistics , volume=
-
[69]
Multivariate functional principal component analysis: A normalization approach , author=. Statistica Sinica , pages=. 2014 , publisher=
work page 2014
-
[70]
Journal of the Royal Statistical Society Series C: Applied Statistics , volume=
The functional latent block model for the co-clustering of electricity consumption curves , author=. Journal of the Royal Statistical Society Series C: Applied Statistics , volume=. 2018 , publisher=
work page 2018
-
[71]
Journal of the Royal Statistical Society Series B: Statistical Methodology , volume=
Sparse Kronecker product decomposition: a general framework of signal region detection in image regression , author=. Journal of the Royal Statistical Society Series B: Statistical Methodology , volume=. 2023 , publisher=
work page 2023
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.