Abstract:The computation of a similarity measure for genomic data is a standard tool in computational genetics. The principal components of such matrices are routinely used to correct for biases due to confounding by population stratification, for instance in linear regressions. However, the calculation of both a similarity matrix and its singular value decomposition (SVD) are computationally intensive. The contribution of this article is threefold. First, we demonstrate that the calculation of three matrices (called the covariance matrix, the weighted Jaccard matrix, and the genomic relationship matrix) can be reformulated in a unified way which allows for the application of a randomized SVD algorithm, which is faster than the traditional computation. The fast SVD algorithm we present is adapted from an existing randomized SVD algorithm and ensures that all computations are carried out in sparse matrix algebra. The algorithm only assumes that row-wise and column-wise subtraction and multiplication of a vector with a sparse matrix is available, an operation that is efficiently implemented in common sparse matrix packages. An exception is the so-called Jaccard matrix, which does not have a structure applicable for the fast SVD algorithm. Second, an approximate Jaccard matrix is introduced to which the fast SVD computation is applicable. Third, we establish guaranteed theoretical bounds on the accuracy (in norm and angle) between the principal components of the Jaccard matrix and the ones of our proposed approximation, thus putting the proposed Jaccard approximation on a solid mathematical foundation, and derive the theoretical runtime of our algorithm. We illustrate that the approximation error is low in practice and empirically verify the theoretical runtime scalings on both simulated data and data of the 1000 Genome Project.

Fast randomized algorithms for low-rank matrix approximations with applications in global comparative analysis of a class of data sets

Randomized algorithms for low-rank matrix approximation: Design, analysis, and applications

A general error analysis for randomized low-rank approximation with application to data assimilation

Randomized low-rank approximations beyond Gaussian random matrices

Fast computation of the eigensystem of genomic similarity matrices

Online Dominant Generalized Eigenvectors Extraction Via a Randomized Algorithm

Efficient Randomized Algorithms for the Fixed-Precision Low-Rank Matrix Approximation

Fast Spectral Low Rank Matrix Approximation

A Fast Randomized Algorithm for Computing an Approximate Null Space

An algorithm for low-rank matrix factorization and its applications

Fast Randomized Numerical Rank Estimation for Numerically Low-Rank Matrices

Krylov subspace algorithms for computing GeneRank for the analysis of microarray data mining.

A Random Algorithm for Low-Rank Decomposition of Large-Scale Matrices with Missing Entries

Fast low-rank estimation by projected gradient descent: General statistical and algorithmic guarantees

An Efficient Gauss--Newton Algorithm for Symmetric Low-Rank Product Matrix Approximations

Fast and Accurate Randomized Algorithms for Linear Systems and Eigenvalue Problems

A Randomized Algorithm for Approximating Truncated SVD

The Rank and Singular Values of the Inhomogeneous Subgaussian Random Matrices

Randomized algorithms for low-rank matrix factorizations: sharp performance bounds

Randomized algorithms for total least squares problems.