# Canonical correlation

In statistics, canonical-correlation analysis (CCA) is a way of making sense of cross-covariance matrices. If we have two vectors X = (X1, ..., Xn) and Y = (Y1, ..., Ym) of random variables, and there are correlations among the variables, then canonical-correlation analysis will find linear combinations of the Xi and Yj which have maximum correlation with each other. T. R. Knapp notes "virtually all of the commonly encountered parametric tests of significance can be treated as special cases of canonical-correlation analysis, which is the general procedure for investigating the relationships between two sets of variables." The method was first introduced by Harold Hotelling in 1936.

## Definition

Canonical-correlation analysis seeks vectors $a$ and $b$ such that the random variables $a'X$ and $b'Y$ maximize the correlation $\rho =\operatorname {corr} (a'X,b'Y)$ . The random variables $U=a'X$ and $V=b'Y$ are the first pair of canonical variables. Then one seeks vectors maximizing the same correlation subject to the constraint that they are to be uncorrelated with the first pair of canonical variables; this gives the second pair of canonical variables. This procedure may be continued up to $\min\{m,n\}$ times.

## Computation

### Derivation

$\rho ={\frac {a'\Sigma _{XY}b}{{\sqrt {a'\Sigma _{XX}a}}{\sqrt {b'\Sigma _{YY}b}}}}.$ The first step is to define a change of basis and define

$c=\Sigma _{XX}^{1/2}a,$ $d=\Sigma _{YY}^{1/2}b.$ And thus we have

$\rho ={\frac {c'\Sigma _{XX}^{-1/2}\Sigma _{XY}\Sigma _{YY}^{-1/2}d}{{\sqrt {c'c}}{\sqrt {d'd}}}}.$ By the Cauchy-Schwarz inequality, we have

$\left(c'\Sigma _{XX}^{-1/2}\Sigma _{XY}\Sigma _{YY}^{-1/2}\right)d\leq \left(c'\Sigma _{XX}^{-1/2}\Sigma _{XY}\Sigma _{YY}^{-1/2}\Sigma _{YY}^{-1/2}\Sigma _{YX}\Sigma _{XX}^{-1/2}c\right)^{1/2}\left(d'd\right)^{1/2},$ $\rho \leq {\frac {\left(c'\Sigma _{XX}^{-1/2}\Sigma _{XY}\Sigma _{YY}^{-1}\Sigma _{YX}\Sigma _{XX}^{-1/2}c\right)^{1/2}}{\left(c'c\right)^{1/2}}}.$ There is equality if the vectors $d$ and $\Sigma _{YY}^{-1/2}\Sigma _{YX}\Sigma _{XX}^{-1/2}c$ are collinear. In addition, the maximum of correlation is attained if $c$ is the eigenvector with the maximum eigenvalue for the matrix $\Sigma _{XX}^{-1/2}\Sigma _{XY}\Sigma _{YY}^{-1}\Sigma _{YX}\Sigma _{XX}^{-1/2}$ (see Rayleigh quotient). The subsequent pairs are found by using eigenvalues of decreasing magnitudes. Orthogonality is guaranteed by the symmetry of the correlation matrices.

### Solution

The solution is therefore:

Reciprocally, there is also:

Reversing the change of coordinates, we have that

The canonical variables are defined by:

$U=c'\Sigma _{XX}^{-1/2}X=a'X$ $V=d'\Sigma _{YY}^{-1/2}Y=b'Y$ ### Implementation

CCA can be computed using singular value decomposition on a correlation matrix. It is available as a function in

## Hypothesis testing

Each row can be tested for significance with the following method. Since the correlations are sorted, saying that row $i$ is zero implies all further correlations are also zero. If we have $p$ independent observations in a sample and ${\widehat {\rho }}_{i}$ is the estimated correlation for $i=1,\dots ,\min\{m,n\}$ . For the $i$ th row, the test statistic is:

$\chi ^{2}=-\left(p-1-{\frac {1}{2}}(m+n+1)\right)\ln \prod _{j=i}^{\min\{m,n\}}(1-{\widehat {\rho }}_{j}^{2}),$ which is asymptotically distributed as a chi-squared with $(m-i+1)(n-i+1)$ degrees of freedom for large $p$ . Since all the correlations from $\min\{m,n\}$ to $p$ are logically zero (and estimated that way also) the product for the terms after this point is irrelevant.

## Practical uses

A typical use for canonical correlation in the experimental context is to take two sets of variables and see what is common amongst the two sets. For example in psychological testing, you could take two well established multidimensional personality tests such as the Minnesota Multiphasic Personality Inventory (MMPI-2) and the NEO. By seeing how the MMPI-2 factors relate to the NEO factors, you could gain insight into what dimensions were common between the tests and how much variance was shared. For example you might find that an extraversion or neuroticism dimension accounted for a substantial amount of shared variance between the two tests.

One can also use canonical-correlation analysis to produce a model equation which relates two sets of variables, for example a set of performance measures and a set of explanatory variables, or a set of outputs and set of inputs. Constraint restrictions can be imposed on such a model to ensure it reflects theoretical requirements or intuitively obvious conditions. This type of model is known as a maximum correlation model.

Visualization of the results of canonical correlation is usually through bar plots of the coefficients of the two sets of variables for the pairs of canonical variates showing significant correlation. Some authors suggest that they are best visualized by plotting them as heliographs, a circular format with ray like bars, with each half representing the two sets of variables.