eigenResidop {scMerge} | R Documentation |
fast matrix residual operator using RcppEigen
eigenResidop(A, B)
A |
a matrix |
B |
a matrix |
The matrix product of
A - B(B^t B)^{-1} B^t A
Y = M = diag(1, 500) system.time(scMerge::eigenResidop(Y, M)) system.time(ruv::residop(Y, M))