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
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)