powered by
uses random matrix to estimate cca results
randcca(x, y, k, seed)
input matrix
rank to use
for testing
outputs a list containing:
whitened: low-rank whitened joint matrix
svd: low-rank svd of joint matrix
# NOT RUN { 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) # }
Run the code above in your browser using DataLab