From Common Information to Canonical Correlation Analysis and Beyond: Bridging Information Theory, Signal Processing, and Machine Learning
Authors: Paris A. Karakasis, Walid Saad
Organizations: Institute for Advanced Computing, Bradley Department of Electrical and Computer Engineering, Virginia Polytechnic Institute and State University, Alexandria, VA 22305, USA
Many problems in signal processing, machine learning, and neuroimaging provide several jointly observed measurements of one process, each mixing a component common to all of them with a component private to itself, and require the common component to be recovered without supervision. This problem has been studied in several communities, largely independently. Statisticians proposed canonical correlation analysis (CCA), introduced by Hotelling in 1936. Information theorists proposed two notions of common information, due to Gács and Körner and to Wyner. Engineers developed subspace methods for sensor arrays and multiuser communications to design multiview receivers. This tutorial presents the main results and tools each community developed for this problem in a single exposition, and explains the connections that emerge where they meet, distinguishing exact equivalences from special cases and analogies. We develop the information-theoretic foundation from first principles. We then introduce CCA and its extension to several views, followed by the engineering thread from blind multiuser detection to CCA as a blind receiver. Next, we turn to nonlinear CCA, covering kernel, deep, and probabilistic formulations, the conditions under which the common factor is identifiable, and connections to self-supervised learning. We conclude with a discussion of further extensions and open directions. The exposition is self-contained, so that readers from each of these communities can follow the parts developed by the others.
Figures & tables
Figure 1: Two neighboring cells, each served by a multi-antenna base station. A user near the center of a cell is received strongly at its own station and is barely heard at the other, so only one link is drawn for it. The two cell-edge users are received at low power by both stations, so each of them carries a link to BS-1 and a link to BS-2, and they are the challenging case for decoding.
Work
Info. Theory
Statistics
Sig. Proc.
Multi-View CCA
Nonlinear CCA
Deep CCA
Identif.
Borga [ 14 ]
Partial
✓
Partial
Hardoon et al. [ 15 ]
✓
Partial
✓
Uurtio et al. [ 16 ]
✓
✓
Partial
Wang et al. [ 17 ]
✓
Partial
Partial
Partial
Partial
Yang et al. [ 23 ]
✓
✓
✓
✓
Partial
Yu and Tan [ 24 ]
✓
Table 1: Scope of existing treatments of common information and CCA. A check mark means the topic is covered substantively, and Partial means it is touched on without being developed.
Term
Meaning
Shared structure
Informal term for whatever the views have in common, used only before one of the definitions below is adopted.
Common information
A number, measured in bits and computed from the joint distribution of the views, namely the Gács–Körner or the Wyner common information. The number says how much the views share, not what they share.
Common factor
A latent random vector shared by all views in a multiview generative model, such that the views are conditionally independent given the common factor. Whatever else enters a view is private to that view and forms its private factor.
Common subspace
The column space of the matrix whose rows are the values of the common factor at the training samples.
Common representation
What the outputs of the learned encoders of the views are made to agree on.
Table 2: Terminology used throughout the tutorial.
Figure 2: Two perception systems at an intersection label the same objects with different label sets. A third system maps the labels of System A through f and the labels of System B through g to a common label set. The finest labeling on which the two maps always agree is light vehicle, heavy vehicle and person.
Figure 3: One speaker recorded by two microphones. Each recording contains the voice, which is common to both, and the noise of its own microphone. The recordings are dependent only through the voice, so given the voice the two recordings are conditionally independent.
Figure 4: Partitions of the sample space Ω=[0,1]2 , with outcome ω=(ω1,ω2) , into the atoms of σ(X) and σ(Y) , the sigma-algebras generated by X and Y , of their join σ(X)∨σ(Y)=σ(X,Y) , and of their meet σ(X)∩σ(Y) . Each color marks one atom, that is, the set of all outcomes at which the random variable generating the partition takes one particular value. The join is the coarsest partition that refines both σ(X) and σ(Y) . The meet is the finest partition coarser than both, and the five atoms of the meet index the values of the common part W∗ .
Figure 5: Joint PMF PXY in the case KGK(X;Y)=I(X;Y)=C(X;Y) , with rows and columns permuted so that PXY is block-diagonal. The block in blue corresponds to W∗=1 , the block in orange to W∗=2 , and darker shades indicate larger probabilities. Each block equals PW∗(w) times the outer product of the conditional PMFs PX∣W∗=w and PY∣W∗=w , as the arrows indicate, which reflects the conditional independence of X and Y given W∗ . The two blocks occupy disjoint rows and disjoint columns, so each of X and Y determines the common part W∗ from its own value. Hence, X always reveals the block in which Y lies, but nothing more about Y , and vice versa.
Figure 6: The binary symmetric channel with crossover probability a0 that relates X and Y in the doubly symmetric binary source. Each bit of X is received correctly as Y with probability 1−a0 and flipped with probability a0 .
Figure 7: Wyner’s optimal common variable W for the doubly symmetric binary source, shown as the intermediate variable in a cascade of two binary symmetric channels, each with crossover probability a1=21(1−1−2a0) . Since 2a1(1−a1)=a0 , the cascade is equivalent to the single binary symmetric channel of Figure 6 , and X and Y are conditionally independent given W .
Figure 8: The Gray–Wyner source coding network. An encoder observes (Xn,Yn) and describes them to two receivers using a common link of rate R0 and two private links of rates R1 and R2 . Wyner’s common information C(X;Y) is the minimum R0 among achievable rate triples with R0+R1+R2=H(X,Y) .
Figure 9: The naive Bayes model. A single latent variable W generates every observation, and no edge joins the Xk to one another, so they are mutually independent once W is given. This is the conditional independence ( 17 ) that Wyner’s construction requires, and writing it out on a discrete joint PMF gives the nonnegative tensor factorization above, with one rank-one factor per value of W .
Figure 10: The single-component case, drawn for X1,X2∈R3×2 . Blue refers to view 1 and green to view 2, for the planes, the circles, and the vectors. The constraints confine each canonical variate to the unit sphere of RN , SN−1 , and to the column space of its own view. Hence, the two feasible sets are the circles where the sphere intersects those subspaces, and θ is the angle between a feasible pair. Minimizing ∥X1q1−X2q2∥22 closes that angle as far as the subspaces allow. Here, two planes in R3 share a line, so the circles intersect at the two antipodal points marked in red, each corresponding to an optimum, where θ=0 and σ1=1 . For visualization purposes, the figure ignores centering. With N=3 , centered columns lie in the two-dimensional subspace orthogonal to 1 , so two centered data matrices of rank two would share the same column space.
Figure 11: Both canonical pairs at once, drawn for X1,X2∈R3×2 . Along the line the two subspaces share, the canonical variates coincide, so θ1=0 and σ1=1 , and this direction spans col(X1)∩col(X2) . The second pair cannot be brought into agreement, since each variate is confined to its own subspace and must be uncorrelated with the first, so it settles at the angle θ2 with σ2=cosθ2<1 . What CCA reports as common is the part of the spectrum at σi=1 . As in Figure 10 , the figure ignores centering.
Receiver
Designed Knowledge
Propagated Knowledge
Cost
CDMA, single antenna
Zero-forcing [ 104 ]
all signatures
none
detected users spread, rate 1/M
Linear MMSE, closed form [ 105 ]
all signatures
received amplitudes, noise variance
detected users spread, rate 1/M
Subspace blind multiuser detection [ 7 ]
own signature
none
detected users spread, rate 1/M
Multiple antennas
Pilot-based MMSE (Section 4.3 )
pilots
none
pilot overhead, rate 1−Np/N
Table 3: Prior information required by each receiver, and its cost.
Figure 12: The two views of the speech example over the same time interval. ( a ) The acoustic view, summarized by the first four mel-frequency cepstral coefficients (MFCCs), a standard description of the short-time spectrum of the sound. ( b ) The articulatory view, namely the trajectories of sensors on the tongue, the upper and lower lips, and the soft palate. The two views bear no obvious resemblance to each other, yet both change at the same moments, since the phoneme being produced shapes both.
Figure 13: The two systems of assumptions, drawn for K=2 . Each view is produced from the common factor and its own private factor through the mixing vk of ( 116 ). In ( a ) the three latent factors are mutually independent, which is ( 118 ) and the assumption of Lyu et al. In ( b ) two further edges are added, so each private factor may depend on the common one, which is ( 122 ) and the assumption of Karakasis and Sidiropoulos. The views stay conditionally independent given g in both cases, since every path between x1 and x2 passes through g .