powered by
make SVD as well as orthogonal complements
sk.decompose(A, D)
U
S
V
U_perp: orthogonal complement for U
the input matrix
the linear transform
library(mvtnorm) n = 350 p = 100 D <- diag(p) Sigma = matrix(0, p, p) X <- rmvnorm(n,matrix(0, p, 1), Sigma) decompose.result <- sk.decompose(X, D) U_perp <- decompose.result$U_perp
Run the code above in your browser using DataLab