if (FALSE) {
set.seed(789)
N <- 50; m.true <- 2; n_batch <- 5
batches <- lapply(1:n_batch, function(i) matrix(rnorm(60 * N), 60, N))
inc <- EFM::incremental_perturb_pca(batches, m = m.true)
print(round(inc$Ap[1:5, 1:2], 3))
}
Run the code above in your browser using DataLab