1525695618CanonicalCorrelation 1
1525695618CanonicalCorrelation 1
1525695618CanonicalCorrelation 1
1 / 15
Development Team
2 / 15
Development Team
2 / 15
Development Team
2 / 15
Development Team
2 / 15
Correlations
3 / 15
Correlations
3 / 15
Correlations
3 / 15
Correlations
3 / 15
Correlations
3 / 15
Correlations
3 / 15
Correlations
3 / 15
Canonical Correlation
4 / 15
Canonical Correlation
4 / 15
Canonical Correlation
4 / 15
Canonical Correlation
4 / 15
The Model
5 / 15
The Model
5 / 15
The Model
5 / 15
The Model
5 / 15
The Problem
6 / 15
The Problem
6 / 15
The Method
U = a0 X1 and V = b0 X2 .
so that
a0 Σ12 b
Corr(U, V ) = √ .
a0 Σ11 a × b0 Σ22 b
7 / 15
The Method
U = a0 X1 and V = b0 X2 .
so that
a0 Σ12 b
Corr(U, V ) = √ .
a0 Σ11 a × b0 Σ22 b
7 / 15
How to choose a and b
8 / 15
How to choose a and b
8 / 15
How to choose a and b
8 / 15
How to choose a and b
8 / 15
How to choose a and b
8 / 15
Result
9 / 15
Result
9 / 15
Result
9 / 15
Result
9 / 15
Result (contd.)
Result
In general, for j = 1, . . . , r,
p
Corr(Uj , Vj ) = ρj = λj ,
−1/2 −1/2
with aj = Σ11 pj and bj = Σ22 qj .
10 / 15
Result (contd.)
Result
In general, for j = 1, . . . , r,
p
Corr(Uj , Vj ) = ρj = λj ,
−1/2 −1/2
with aj = Σ11 pj and bj = Σ22 qj .
10 / 15
Proof
11 / 15
Proof (contd.)
−1/2 −1/2
I Σ11 Σ12 Σ−122 Σ21 Σ11 is symmetric, with largest eigen-value
λ1 and corresponding eigen-vector p1 . Thus
−1/2 −1/2
c0 Σ11 Σ12 Σ−1
22 Σ21 Σ11 c ≤ λ 1 c0 c = λ 1 .
I Equality occurs at c = p1 .
I Also equality in the Cauchy-Schwarz inequality occurs if d is
−1/2 −1/2
proportional to Σ22 Σ21 Σ11 p1 .
I This means b is proportional to
−1/2 −1/2 −1/2 −1/2
Σ22 Σ22 Σ21 Σ11 p1 = Σ22 q1 (say),
−1/2 −1/2
where, q1 = Σ22 Σ21 Σ11 p1 .
12 / 15
Proof (contd.)
−1/2 −1/2
I Σ11 Σ12 Σ−122 Σ21 Σ11 is symmetric, with largest eigen-value
λ1 and corresponding eigen-vector p1 . Thus
−1/2 −1/2
c0 Σ11 Σ12 Σ−1
22 Σ21 Σ11 c ≤ λ 1 c0 c = λ 1 .
I Equality occurs at c = p1 .
I Also equality in the Cauchy-Schwarz inequality occurs if d is
−1/2 −1/2
proportional to Σ22 Σ21 Σ11 p1 .
I This means b is proportional to
−1/2 −1/2 −1/2 −1/2
Σ22 Σ22 Σ21 Σ11 p1 = Σ22 q1 (say),
−1/2 −1/2
where, q1 = Σ22 Σ21 Σ11 p1 .
12 / 15
Proof (contd.)
−1/2 −1/2
I Σ11 Σ12 Σ−122 Σ21 Σ11 is symmetric, with largest eigen-value
λ1 and corresponding eigen-vector p1 . Thus
−1/2 −1/2
c0 Σ11 Σ12 Σ−1
22 Σ21 Σ11 c ≤ λ 1 c0 c = λ 1 .
I Equality occurs at c = p1 .
I Also equality in the Cauchy-Schwarz inequality occurs if d is
−1/2 −1/2
proportional to Σ22 Σ21 Σ11 p1 .
I This means b is proportional to
−1/2 −1/2 −1/2 −1/2
Σ22 Σ22 Σ21 Σ11 p1 = Σ22 q1 (say),
−1/2 −1/2
where, q1 = Σ22 Σ21 Σ11 p1 .
12 / 15
Proof (contd.)
−1/2 −1/2
I Σ11 Σ12 Σ−122 Σ21 Σ11 is symmetric, with largest eigen-value
λ1 and corresponding eigen-vector p1 . Thus
−1/2 −1/2
c0 Σ11 Σ12 Σ−1
22 Σ21 Σ11 c ≤ λ 1 c0 c = λ 1 .
I Equality occurs at c = p1 .
I Also equality in the Cauchy-Schwarz inequality occurs if d is
−1/2 −1/2
proportional to Σ22 Σ21 Σ11 p1 .
I This means b is proportional to
−1/2 −1/2 −1/2 −1/2
Σ22 Σ22 Σ21 Σ11 p1 = Σ22 q1 (say),
−1/2 −1/2
where, q1 = Σ22 Σ21 Σ11 p1 .
12 / 15
Proof (contd.)
13 / 15
Proof (contd.)
13 / 15
Properties
−1/2 −1/2
I E(Uj ) = p0j Σ11 µ1 and E(Vj ) = q0j Σ22 µ2 for all j.
−1/2 −1/2
I V ar(Uj ) = p0j Σ11 Σ11 Σ11 pj = p0j pj = 1 for all j.
−1/2 −1/2
I V ar(Vj ) = q0j Σ22 Σ22 Σ22 qj = q0j qj = 1 for all j.
−1/2 −1/2
I Cov(Uj , Uk ) = p0j Σ11 Σ11 Σ11 pk = p0j pk = 0 for all
k 6= j.
−1/2 −1/2
I Cov(Vj , Vk ) = q0j Σ22 Σ22 Σ22 qk = q0j qk = 0 for all k 6= j.
−1/2 −1/2
I Cov(Uj , Vk ) = p0j Σ11 Σ12 Σ22 qk = p0j pk = 0 for all
j 6= k.
14 / 15
Properties
−1/2 −1/2
I E(Uj ) = p0j Σ11 µ1 and E(Vj ) = q0j Σ22 µ2 for all j.
−1/2 −1/2
I V ar(Uj ) = p0j Σ11 Σ11 Σ11 pj = p0j pj = 1 for all j.
−1/2 −1/2
I V ar(Vj ) = q0j Σ22 Σ22 Σ22 qj = q0j qj = 1 for all j.
−1/2 −1/2
I Cov(Uj , Uk ) = p0j Σ11 Σ11 Σ11 pk = p0j pk = 0 for all
k 6= j.
−1/2 −1/2
I Cov(Vj , Vk ) = q0j Σ22 Σ22 Σ22 qk = q0j qk = 0 for all k 6= j.
−1/2 −1/2
I Cov(Uj , Vk ) = p0j Σ11 Σ12 Σ22 qk = p0j pk = 0 for all
j 6= k.
14 / 15
Properties
−1/2 −1/2
I E(Uj ) = p0j Σ11 µ1 and E(Vj ) = q0j Σ22 µ2 for all j.
−1/2 −1/2
I V ar(Uj ) = p0j Σ11 Σ11 Σ11 pj = p0j pj = 1 for all j.
−1/2 −1/2
I V ar(Vj ) = q0j Σ22 Σ22 Σ22 qj = q0j qj = 1 for all j.
−1/2 −1/2
I Cov(Uj , Uk ) = p0j Σ11 Σ11 Σ11 pk = p0j pk = 0 for all
k 6= j.
−1/2 −1/2
I Cov(Vj , Vk ) = q0j Σ22 Σ22 Σ22 qk = q0j qk = 0 for all k 6= j.
−1/2 −1/2
I Cov(Uj , Vk ) = p0j Σ11 Σ12 Σ22 qk = p0j pk = 0 for all
j 6= k.
14 / 15
Properties
−1/2 −1/2
I E(Uj ) = p0j Σ11 µ1 and E(Vj ) = q0j Σ22 µ2 for all j.
−1/2 −1/2
I V ar(Uj ) = p0j Σ11 Σ11 Σ11 pj = p0j pj = 1 for all j.
−1/2 −1/2
I V ar(Vj ) = q0j Σ22 Σ22 Σ22 qj = q0j qj = 1 for all j.
−1/2 −1/2
I Cov(Uj , Uk ) = p0j Σ11 Σ11 Σ11 pk = p0j pk = 0 for all
k 6= j.
−1/2 −1/2
I Cov(Vj , Vk ) = q0j Σ22 Σ22 Σ22 qk = q0j qk = 0 for all k 6= j.
−1/2 −1/2
I Cov(Uj , Vk ) = p0j Σ11 Σ12 Σ22 qk = p0j pk = 0 for all
j 6= k.
14 / 15
Properties
−1/2 −1/2
I E(Uj ) = p0j Σ11 µ1 and E(Vj ) = q0j Σ22 µ2 for all j.
−1/2 −1/2
I V ar(Uj ) = p0j Σ11 Σ11 Σ11 pj = p0j pj = 1 for all j.
−1/2 −1/2
I V ar(Vj ) = q0j Σ22 Σ22 Σ22 qj = q0j qj = 1 for all j.
−1/2 −1/2
I Cov(Uj , Uk ) = p0j Σ11 Σ11 Σ11 pk = p0j pk = 0 for all
k 6= j.
−1/2 −1/2
I Cov(Vj , Vk ) = q0j Σ22 Σ22 Σ22 qk = q0j qk = 0 for all k 6= j.
−1/2 −1/2
I Cov(Uj , Vk ) = p0j Σ11 Σ12 Σ22 qk = p0j pk = 0 for all
j 6= k.
14 / 15
Properties
−1/2 −1/2
I E(Uj ) = p0j Σ11 µ1 and E(Vj ) = q0j Σ22 µ2 for all j.
−1/2 −1/2
I V ar(Uj ) = p0j Σ11 Σ11 Σ11 pj = p0j pj = 1 for all j.
−1/2 −1/2
I V ar(Vj ) = q0j Σ22 Σ22 Σ22 qj = q0j qj = 1 for all j.
−1/2 −1/2
I Cov(Uj , Uk ) = p0j Σ11 Σ11 Σ11 pk = p0j pk = 0 for all
k 6= j.
−1/2 −1/2
I Cov(Vj , Vk ) = q0j Σ22 Σ22 Σ22 qk = q0j qk = 0 for all k 6= j.
−1/2 −1/2
I Cov(Uj , Vk ) = p0j Σ11 Σ12 Σ22 qk = p0j pk = 0 for all
j 6= k.
14 / 15
Summary
15 / 15
Summary
15 / 15
Summary
15 / 15