add centering

HighlanderLab · Dec 27, 2024 · 38882a3 · 38882a3
1 parent 34b8ce0
commit 38882a3
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/randPedPCA/R/rppca.R b/randPedPCA/R/rppca.R
@@ -54,6 +54,7 @@ randSVD <- function(L, rank, depth, numVectors){
   # numVectors: usually rank + 5~10, must be larger than rank
   dim <- nrow(L)
   Q <- randRangeFinder(L, rank, depth, numVectors)
+  # Q <- apply(Q, 2, function(col) col - mean(col))
   C <- t(backsolve(t(L), Q))
   svdObject <- svd(C)
   U <- Q %*% svdObject$u