Cited by 20Random Matrices and ApplicationsSparse and Compressive Sensing TechniquesStatistical Methods and Inference

Biwhitening Reveals the Rank of a Count Matrix

Open access
Landa, Boris, Zhang, Thomas T. C. K., Kluger, Yuval·2022·SIAM Journal on Mathematics of Data Science

Abstract

Abstract. Estimating the rank of a corrupted data matrix is an important task in data analysis, most notably for choosing the number of components in principal component analysis. Significant progress on this task was achieved using random matrix theory by characterizing the spectral properties of large noise matrices. However, utilizing such tools is not straightforward when the data matrix consists of count random variables, e.g., a Poisson matrix, in which case the noise can be heteroskedastic with an unknown variance in each entry. In this work, we focus on a Poisson random matrix with independent entries and propose a simple procedure, termed biwhitening, for estimating the rank of the underlying signal matrix (i.e., the Poisson parameter matrix) without any prior knowledge. Our approach is based on the key observation that one can scale the rows and columns of the data matrix simultaneously so that the spectrum of the corresponding noise agrees with the standard Marchenko–Pastur (MP) law, justifying the use of the MP upper edge as a threshold for rank selection. Importantly, the required scaling factors can be estimated directly from the observations by solving a matrix scaling problem via the Sinkhorn–Knopp algorithm. Aside from the Poisson, our approach is extended to families of distributions that satisfy a quadratic relation between the mean and the variance, such as the generalized Poisson, binomial, negative binomial, gamma, and many others. This quadratic relation can also account for missing entries in the data. We conduct numerical experiments that corroborate our theoretical findings, and showcase the advantage of our approach for rank estimation in challenging regimes. Furthermore, we demonstrate the favorable performance of our approach on several real datasets of single-cell RNA sequencing, high-throughput chromosome conformation capture, and document topic modeling.

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1137/21m1456807

Cites

95 of 95 references verified0 citation errors · 0 not checkable

Statistics

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Cites

Flagged first, then by citations

Research works printed but not linked: 1

  1. [40]Jackson JE, A user’s guide to principal components, vol. 587, John Wiley & Sons, 2005.

Cited by

Sorted by citations

Authors
Landa, B., Zhang, T. T. C. K., & Kluger, Y.
Affiliations
  • Yale University (US)
  • University of Pennsylvania (US)
Venue
SIAM Journal on Mathematics of Data Science
Published
2022-12-20
Volume
4
Issue
4
Pages
1420-1446
ISSN
2577-0187 (ISSN-L)
OpenAlex
W3139046810
Cited by
20
References
94
Citations by year
2026: 3 · 2025: 7 · 2024: 3 · 2023: 4 · 2021: 3
FWCI
2.58
Type
journal article
Language
en
Open access
Locations
  • SIAM Journal on Mathematics of Data Science (publishedVersion)
  • PubMed (publishedVersion)
  • PubMed Central (acceptedVersion)
  • arXiv (Cornell University) (submittedVersion)
  • arXiv (Cornell University)
Funders
National Institutes of Health (R01GM131642, UM1DA051410, U01DA053628, U54AG076043, R01GM135928, 2P50CA1219)
Topics
  • Physical Sciences › Mathematics › Statistics and Probability › Random Matrices and Applications
  • Physical Sciences › Engineering › Computational Mechanics › Sparse and Compressive Sensing Techniques
  • Physical Sciences › Mathematics › Statistics and Probability › Statistical Methods and Inference
Keywords
rank estimation, Marchenko-Pastur law, heteroskedastic noise, matrix scaling, Sinkhorn-Knopp algorithm, principal component analysis, missing data, single-cell RNA sequencing, chromosome conformation capture, topic modeling

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review