Abstract:Non-negative matrix factorization (NMF) is a popular method for analyzing strictly positive data due to its relatively straightforward interpretation. However, NMF has a reputation as a less efficient alternative to the singular value decomposition (SVD), a standard operation that is highly optimized in most linear algebra libraries. Sparse single-cell sequencing assays, now feasible in thousands of subjects and millions of cells, generate data matrices with tens of thousands of strictly non-negative transcript abundance entries. We present an extremely fast NMF implementation made available in the RcppML (Rcpp Machine Learning library) R package that rivals the runtimes of state-of-the-art Singular Value Decomposition (SVD). NMF can now be run quickly on desktop computers to analyze sparse single-cell datasets consisting of hundreds of thousands of cells. Our method improves upon current NMF implementations by introducing a scaling diagonal to increase interpretability, guarantee consistent regularization penalties across different random initializations, and symmetry in symmetric factorizations. We use our method to show how NMF models learned on standard log-normalized count data are interpretable dimensional reductions, describe interpretable patterns of coordinated gene activities, and explain biologically relevant metadata. We believe NMF has the potential to replace PCA in most single-cell analyses, and the presented NMF implementation overcomes previous challenges with long runtime.

Nonnegative Singular Value Decomposition For Microarray Data Analysis Of Spermatogenesis

Unveiling Molecular Signatures for Disease Classification by Spectroscopic Imaging

Sparse nonnegative matrix factorization applied to microarray data sets

Non-negative Matrix-Set Factorization

Sparse Nonnegative Matrix Factorization with Genetic Algorithms for Microarray Analysis

Fast and interpretable non-negative matrix factorization for atlas-scale single cell data

Sparse P-Norm Nonnegative Matrix Factorization for Clustering Gene Expression Data.

On -Divergence Based Nonnegative Matrix Factorization for Clustering Cancer Gene Expression Data

On α-divergence based nonnegative matrix factorization for clustering cancer gene expression data

New SVD based initialization strategy for Non-negative Matrix Factorization

Sparse Orthogonal Nonnegative Matrix Factorization for Identifying Differentially Expressed Genes and Clustering Tumor Samples.

Nonnegative Matrix Factorization on Orthogonal Subspace

Reducing Microarray Data Via Nonnegative Matrix Factorization for Visualization and Clustering Analysis

Tumor Classification Using Non-negative Matrix Factorization

LS-NMF: A modified non-negative matrix factorization algorithm utilizing uncertainty estimates

A Provable Splitting Approach for Symmetric Nonnegative Matrix Factorization

A unified statistical approach to non-negative matrix factorization and probabilistic latent semantic indexing

Non-negative matrix factorization with α-divergence

Blind Matrix Decomposition Via Genetic Optimization Of Sparseness And Nonnegativity Constraints

Nonnegative Discriminant Matrix Factorization