The answer is no, you could extract the components (canonical variates) from the rcc result and use the base stat function cor.test().
Thanks! is the
No because of the regularization?
If so, if I don’t apply regularization because of
n > p + q, i.e. setting
lambda 1 and 2 as 0, should the function be the right one?