powered by
Robust PCA by adding random symmetric perturbation to covariance matrix.
robust_perturb_pca(data, m, epsilon = 0.001, sigma = 0.01)
List with Ap, Dp, lambda.
T x N matrix.
Components.
Perturbation scale.
Noise sd.
if (FALSE) { set.seed(456) N <- 80; T <- 100; m.true <- 2 X <- matrix(rnorm(T * N), T, N) rpca.out <- EFM::robust_perturb_pca(X, m = m.true) print(round(rpca.out$Ap[1:5, 1:2], 3)) }
Run the code above in your browser using DataLab