randcca.Rd
uses random matrix to estimate cca results
randcca(x, y, k, seed = NA)
input matrix
rank to use
for testing
outputs a list containing:
whitened: low-rank whitened joint matrix
svd: low-rank svd of joint matrix
Avants BB
set.seed(13) x <- matrix(rnorm(3000), nrow = 50) y <- x %*% matrix(rnorm(60 * 100), nrow = 60) k <- 10 dr <- randcca(x, y, k, 1)