QRM 06
QRM 06
Proof.
“⇒” As we have seen in (17), a covariance matrix Σ is positive semidefinite.
“⇐” Let Σ be positive semidefinite with Cholesky factor A. Let X be
ind.
a random vector with cov X = Id = diag(1, . . . , 1) (e.g. Xj ∼
N(0, 1)). Then cov(AX) = A cov(X)A′ = AA′ = Σ, i.e. Σ is a
covariance matrix (namely that of AX).
◮ Note that
n
′1X
E((Xi − µ)(X̄ − µ) ) = E((Xk − µ)(X̄ − µ)′ )
n k=1
n
1X
=E (Xk − µ)(X̄ − µ)′
n k=1
Σ
= E((X̄ − µ)(X̄ − µ)′ ) = cov(X̄) = .
n
Typically k = d
EX = µ + AEZ = µ
cov(X) = cov(µ + AZ) = A cov(Z)A′ = AA′ =: Σ
Consequences: a=ej
Margins: X ∼ Nd (µ, Σ) ⇒ Xj ∼ N(µj , Σjj ), j ∈ {1, . . . , d}.
:
a=1 Pd Pd Pd
Sums: X ∼ Nd (µ, Σ) ⇒ j=1 Xj ∼ N( j=1 µj , i,j=1 Σij ).
Consequences:
Sets of the form Sc = {x ∈ Rd : (x − µ)′ Σ−1 (x − µ) = c}, c > 0,
describe points of equal density. Contours of equal density are thus
ellipsoids. Whenever a multivariate density fX (x) depends on x only
through the quadratic form (x − µ)′ Σ−1 (x − µ), it is the density of an
elliptical distribution (see later).
The components of X ∼ Nd (µ, Σ) are mutually independent if and only
if Σ is diagonal, i.e. if and only if the components of X are uncorrelated.
f (x1, x2)
f (x1, x2)
0.3 0.4
0.2 0.3
0.1 0.2
0 0.1
4 0
2 4 4
0 2 2 4
x2 −2 0 0 2
−2 x1 x2 −2 0
−4 −4 −2 x1
−4 −4
4 4
2 2
x2 0 x2 0
−2 −2
−4 −4
−4 −2 0 2 4 −4 −2 0 2 4
x1 x1
1 −0.7 ν−2
Left: Nd µ, Σ for µ = ( 00 ), Σ = ( −0.7 1 ); Right: tν (µ, ν Σ), ν = 4,
(same mean and covariance matrix as on the left-hand side)
© QRM Tutorial Section 6.1.3
d
The definition of Nd (µ, Σ) in terms of a stochastic representation (X =
µ + AZ) directly justifies the following sampling algorithm.
We can also compare Di2 data to a χ210 graphically using a Q-Q plot.
© QRM Tutorial Section 6.1.4
Q-Q plot of Di2 data against a χ210 distribution:
(a) daily data; (b) weekly data; (c) monthly data; and (d) quarterly data
250 • •
(a) (b)
200 • 60
Ordered D2 data •
150 •
40 •••••
• ••••••
100 ••• ••••
•••
• ••••••
•• 20 •
••••
50 •••• •••••
••••••••• ••••••
•••
••
••
••
••
••••••••
• ••
••
••
••
••
•••••••••••
•
•••••
•••••••••••••••••••••••••••••••••••••••• •••••••••
0 •• 0 ••••••••••••
0 5 10 15 20 25 30 5 10 15 20 25 30
40 (c) • •
20 (d) • •
•
Ordered D2 data
•
30 ••
• 15
•• • •
• •
20 ••• 10 •••
••••••• ••••
•••• •••
•••••• ••••
10
••••••
•
•• ••••
•
••••• 5
•••••••• ••
•
••••
•••••••• ••
5 10 15 20 25 5 10 15 20
2 2
χ10 quantile χ10 quantile
0.10
●
●
0.05
0.05
●
●
●
●
● ● ●
● ● ●
● ●
●
● ● ●
● ● ● ●
●
● ●● ● ●
● ●● ● ● ● ●
● ● ● ● ● ● ● ●
● ● ●
●● ●● ● ● ● ● ● ●
●
●● ● ● ●
● ●● ●
● ●● ●
●● ● ● ●●● ● ● ● ●
● ● ●
● ●● ● ● ● ● ● ●● ● ● ●
● ● ●●
●● ●● ● ● ●●●
● ●
● ●● ●
● ● ● ● ● ● ● ●● ● ●
●● ● ●● ● ●●●● ●
●●● ●● ● ●● ●
●● ● ●● ● ● ●● ● ●
● ● ●
● ●● ●● ●● ● ●
SIEMENS
●● ● ●● ●● ●●● ● ● ●
● ● ●●●●
●●●● ●●
● ● ●● ●● ●
●
●●●●● ●● ● ● ● ● ●
● ●● ● ●
●●●● ●●●
● ●●
●●●● ●●●●● ●
●● ●●
● ●
● ●● ● ● ●● ●● ● ● ● ●
● ●● ●●● ●●●●
●●● ●●●● ● ● ● ●●● ●● ●●●●
●●●●
●●
● ● ● ●● ●● ● ●
● ●● ●● ●● ● ●● ●● ●●● ●● ●
●● ● ● ●●●●●●
●
●●
●
●● ●● ● ●
● ● ●● ● ●●●●
● ●●●
● ●● ● ●●●● ●●●● ●
●
●
● ● ●●●● ●● ●
●
●●●●
● ●● ●●●
● ●● ●●●●●● ● ●
●● ● ● ●● ●●●● ● ● ●●● ● ●
●
●●
●●
●●●
●●
●●●
● ● ●●●●
● ● ●
● ● ● ●
● ● ● ● ●● ●● ● ● ●
● ●● ●●●●●●●● ●●
●●
●●●● ●
●● ● ●● ●● ● ●● ● ● ● ●●●● ●●●● ●●● ●●●●● ●● ● ● ●● ● ●●●● ● ●
●● ●● ● ●● ●●●●● ●
● ●●●●
● ●●
●● ● ●●● ● ● ● ● ●●
●
● ●
●●●●●●● ●●●●● ●●
● ● ● ● ●●
●
●●
●●
● ●● ●●●●
● ● ● ●●● ●●●
●
●
●
●
● ●
● ●●●●
●●●●● ● ●●●●● ● ● ●
● ●●
● ● ● ●
●● ●●
● ●
●
●●
●● ●
●●●●
●●●●●●
●
●● ● ●●●●● ●●●●●●●●●●●●● ● ● ● ● ● ●●●●● ●●● ●●
● ●●●●
● ●●● ●● ●●
●●
●
●
●●● ●
●●●●● ●● ●●● ●●● ●●●● ●
●●● ●
● ●
●●●
●●
●●●●
●● ●●
● ●●●● ●● ●●
●● ●●● ●● ●● ●●●●●●
●●●● ●●●●● ●●●●●
●●
●●●●●●●
●● ●
●●●●●● ●●
●
●●●● ●●
●
●●●●
●●●
●
●
●●
●●
●
●
●●●
●●
●●
●●●●
●●
●●●
●
●
●●●
●●
●●
●
●●●
●●●
● ●●●●●● ● ●●
●●●
● ●● ● ● ● ● ●●
● ●
●●●●●
●●●●
●●
●
● ●
●
●
●
●●
●
●●●●
● ●
●●
●●
● ●●
●
●●
●●
●
● ●
●●● ●●●
● ●● ● ●●● ● ●
● ●● ● ●●● ● ●
● ●●●● ● ●● ● ●● ●●
●
● ●
● ●● ●
●● ●● ● ●
●● ● ●
●●●●
●●●●
● ●
● ●
●●● ● ● ● ●● ● ●● ● ● ●
● ● ●●● ● ●
●
● ● ●●●
●● ●
●●●
●●●●
●●●●●●● ●● ●● ●●●●●
● ●
● ●● ● ●● ● ● ●●● ●●● ●●
● ●● ●●●●●●
● ●●●● ●●● ●● ●●●● ●● ●
● ●● ●● ● ● ●
●●●● ●
●● ●●
●● ● ●● ●●
● ●● ●●●●●●●● ●●●● ●● ● ● ● ●● ●● ● ●
●●●● ●●
● ●
● ●
● ●
●●●●
●
●●
●●
●●●
●
●
● ●● ●
●● ●●● ●●●● ●
● ● ●●● ●● ●
● ●
● ●●●●●●● ●●●● ●●●●
● ●
●
●● ●
●●●●●●
●●
●
●●●●● ● ●●●●
●●● ●●●●
●
●●●●● ●● ● ● ●●●
●● ● ●● ●● ● ●●●●●●●●
●
●
●
●●
●
● ●
●●●
●●
●●●
●●●
●
●●
●
● ●
●
●
●
●
●●●
●
●●
●●●●● ●
●●
●●●
● ●● ●● ● ●●● ●
● ● ●●●
●● ● ●●●●
●●● ●●
●●
●●●●
● ●
●
●
●●
●●● ●●●
● ●●● ●●●
●
●●●● ●● ● ●
●●
● ● ●● ●●
●● ●●
●
●●
● ●●●
●
●
●● ●
●●
●●
●●
●●
●
●
●●
●●
●
●●
●● ●● ●
●
●●●●● ● ●● ●
● ● ●
● ●●● ●●● ● ●●●●
●●●
●
●●●● ●● ●
●
●●
●
● ●●
●
●●●●●
●
●●
● ●●
●
●●●●● ●●
●●
●
●●
●●●●●● ●●●●
● ●●●●●
● ● ● ● ●●●● ●●●●●●●
● ●●
● ●●●●
●
●●
●
●●
●●●
●
●
●
●●
●●
●●
●●●●
●
●
●●
●
●
●
●
●
●●●●
●●● ●
●●● ●●●● ●● ●
● ● ●●● ●
● ●●●
●●● ●●
●●
●●● ●●●●●●
● ●●
●●
● ●●●●●●●
● ●●●●
●●●● ●● ●●
● ●
●●● ● ● ● ● ●● ● ● ●● ●● ●● ●
●●●
●
●●● ●
●●●
●
●●
●●●
●●●
●
●
●●●
●●
●● ●
●●
●
●
●
●●
●●●●
●
●●
●●
● ●
●●●●●● ●● ● ●●●
●●● ● ● ● ● ●●
● ● ● ● ●●● ● ● ● ●●● ●
● ●●●
●
●●
●●
●● ●●● ●
● ● ● ●● ●
●●●●● ●●
● ● ● ● ● ●● ●● ● ●● ●●● ●
●●
●
●●●
●● ●
● ●●●
● ● ●● ●●● ●
●
●●
●●●● ●
●●●● ●●● ●●●●●●●
●
●
●●
●●●●
●●
●●
● ●● ●
●
●●●
●● ● ●
● ●●● ●
●● ● ●●●●●● ● ●●
●
● ●● ●
●●● ●● ● ● ● ●● ● ● ● ● ●●
●
●●●
●
●●
●●●
●
●●
●●
● ●
●
●
●
●●
●●
●●
●●
●
●●
●
●
●
●
●●
●
●
●
●
●
●●●●●●●●●
●
●●●
●●● ●●● ● ● ●
● ●● ●● ● ●●● ● ● ● ● ●●
X2
● ● ●● ● ● ● ● ●●● ● ● ●● ●●●●●
●●● ●
● ●
●●
●● ● ● ●● ●●● ●
●● ●
●● ● ●● ●●●
● ●●● ●●
●●●● ●
●
●●● ●● ●
● ●●●● ● ● ● ● ●● ● ●● ● ●●
●
● ●
●
● ●
●● ●
●●● ●● ● ●● ● ● ●
●●● ●● ●
● ● ● ●● ●
●●●●
● ● ●●●●
●
● ●● ●● ●●●●●●●●●
●●
● ●●●●●
●●●●●●●●● ● ● ●● ● ●
●● ●● ●●
●●
●●●●● ●
●
●●
● ●●
●
●●
●
●●●
●●●
●
●●
●● ●
●●● ●●●●
●● ●
●●●● ●●● ●
● ● ●● ● ●●●● ●● ●● ● ● ●● ●●●●●● ● ●● ●●●
●● ●
●●●●
● ●
●
● ●●
●● ●● ● ●●●● ●●
●● ● ● ●●● ●●●●● ●
●● ● ●● ●● ●
● ●●● ●
●●●●●
●●● ● ● ● ● ●●
● ●●●●● ●●●●●●●●●
●
●●●●●● ● ●
●● ● ●
● ● ● ●●● ●● ● ●●● ●
●●
● ●●● ●
●
●●● ●● ●
●● ● ● ●●● ● ● ● ● ●●● ●● ●
● ● ●●
●
●
●
●●
●●●
●●
●●
●
●
●●
●
●
●● ●● ●● ●●
● ●●●●●●●●● ●●●
●●●●● ●● ●
●● ●●●●●
●●●
●●
●●
●
●
●●
●
●
●●
●
●● ● ●●
●
●● ●● ●●● ●
●
●●
● ● ●●
● ● ● ●
● ●● ●●
●● ●●●●
●●●
●●
●●●●
●
● ●
●●●
●●
●
●
●●
●
●●
●
●●●●
●
● ●●
●
●
●
●
●
●
●
●●
●●
●
●
●
●●●●●●●
●●●●●
●
●●● ● ●●
● ● ● ●
●
●
●●●● ●
●● ● ●● ●
● ●●●● ● ● ● ● ●
● ●
●●
● ● ●● ● ● ● ● ● ● ●● ●
●●
●
● ●
●●●● ●●●● ● ● ●
● ●●●●●● ● ●● ● ●●● ● ●● ●● ●●●● ●● ● ●● ● ●● ● ● ●● ● ● ● ●●●●●● ●●●● ●●●●●
● ●● ●● ● ●
● ●●● ●
●●●●
●
● ●●●●●● ● ●●●●●● ●● ● ● ●● ● ● ●●
●●● ●● ●●●● ● ●●● ●● ●● ● ● ●
●
●
● ●●● ● ●
●● ●●● ●●● ●● ●●●
●●●●
● ●●●● ●
● ●●●
● ●● ●●● ● ●● ●● ● ● ●● ● ● ● ● ● ●● ● ●● ●
●
●●●●
●●●● ●●●● ● ● ●
● ● ●● ●● ● ●● ●
●● ● ●● ●●●●●●●●
● ●●●
●●
● ●●●●● ●●
●● ● ●● ●● ● ●●● ●●
●
● ● ●●
●●●●
● ●●●●●●● ●
● ● ●● ●●●●●●●
● ●●●●● ●●
●● ●●●
●
● ●●●●●●● ●● ●● ● ● ● ●
● ● ●● ●● ●
●● ●● ● ●●●● ● ● ● ●
● ●●● ● ●●●● ● ● ●●● ● ●●● ●● ● ● ●● ●●
●●●●●●
● ●● ●
● ●●●
●● ●●●●●●
●●●
●●● ●● ●●
● ● ●● ●● ● ● ●
● ●● ●●●
● ● ●●●
● ●
●
●
● ●●
● ●
●● ● ● ●●●● ● ● ● ● ●●
● ●● ●
●●
● ● ● ●● ● ● ●
●
●●● ● ●
● ● ● ● ● ● ● ●
● ● ● ●● ●
●● ●
●●●●●● ●● ●● ● ●● ●● ● ●●● ● ●● ●●
●● ●
●●
●● ●●
● ●●● ●●● ●●●● ●● ●●● ●
● ● ● ● ●● ● ●
●
● ●
● ● ●● ●●●
● ● ●●
● ●
● ● ● ●● ● ● ●● ●●
●●●●● ●
● ●
● ● ● ● ●● ●●● ●● ●● ● ● ● ●● ● ●● ● ● ●
● ● ●● ●●●● ● ●● ●● ● ● ● ● ● ● ● ●●● ● ● ● ● ●● ● ●
● ● ● ● ● ● ●
−0.05
−0.05
● ● ● ●
●● ● ● ● ●
● ● ●● ●●●● ●●● ● ● ● ● ● ● ● ●●● ●
●
●● ● ●●● ● ● ● ● ● ●●
● ● ●
● ●● ●
● ● ● ●
● ● ● ● ● ●
● ●
● ●● ●
●
● ●
●
●
● ●
●
●
●
●
●
●
●
●
●
●
−0.15
−0.15
−0.15 −0.05 0.05 0.10 −0.15 −0.05 0.05 0.10
X1 BMW
Q−Q plot for margin 1 (simulated data) Q−Q plot for margin 1 (real data)
●
0.10
0.10
●● ●
●
●
●
●●●
0.05
0.05
● ●
●
● ●
●●
●
●
●
●
●
●
Sample quantiles
Sample quantiles
● ●
●● ●●
●
●
●●●●●●●
● ●
●
●
●
●
●
●
●
●
●●
●
●●
● ●
●●
●
●
●
●●
●
●●
●
●●
●
●● ●●
●
●
●
●●
●●
●
●●
● ●
●
●●
●
●● ●●
●
●
●●
●
●
●●
● ●
●
●
●●
●
●
●●
●
●●
●
●
●●
● ●
●
●●
●
●
●
●
●●
●
●
●●
●● ●●
●
●●
●
●
●●
●
●●
●
●
●●
● ●
●
●●
●
●
●
●
●
●●
●
●
●●
●
● ●●
●
●●
●
●●
●
●
●●
●
● ●●
●
●
●
●
●
●●
●
●
●●
● ●
●●
●
●
●●
●
●
●
●●
●
●
●●
●
●
● ●●
●
●
●●
●
●●
●
●
●
●●
●
●
●●
●
●
●● ●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●● ●
●●
●
●
●●
●
●
●
●
●●
●
●●
●
●
●●
●
● ●
●
●
●●
●
●
●●
●
●
●
●
●●
●
●
●●
● ●●
●
●
●●
●
●
●●
●
●
●
●●
●
●
●●
●
●
● ●
●●
●
●
●●
●
●
●●
●
●
●
●
●●
●
●
●●
●
● ●●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●● ●
●
●●
●
●
●●
●
●
●●
●
●
●
●
●●
●
●
●●
● ●●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●
●
●
●
●●
●
●
●●
●
●
●● ●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
● ●
●●
●
●
●●
●
●
●●
●
●
●
●
●
●●
●
●
●●
●
● ●
●●
●
●
●●
●
●
●●
●
●●
●
●
●●
●
●● ●
●
●●
●
●
●●
●
●
●●
●
●
●
●●
●
●
●●
●
● ●
●●
●
●
●●
●
●
●●
●
●●
●
●
●●
●
●
●● ●
●●
●
●●
●
●
●●
●
●
●
●●
●
●
●●
● ●
●
●●
●
●
●●
●
●
●
●
●●
●
●● ●
●
●●
●
●●
●
●
●
●
●
●●
●
●
●●
●
●
● ●
●●
●
●
●●
●
●●
●
●
●●
●
●
● ●●
●
●●
●
●
●●
●
●
●●
●
●
●● ●
●
●●
●
●
●
●
●
●●
●
●●
●
●
● ●
●
●●
●
●
●
●
●●
●
●
●●
●
● ●
●
●
●
●
●
●●
●
● ●●
−0.05
−0.05
●●
●
●●
●
●
●●
●
● ●●
●
●
●
●●
●
●●
● ●
●●
●
●●
●
●
●●
●
●●
●
●● ●●
●
●●
●●
●
●●
●
●●
● ●●
●●●
●●
●● ●●
●
●
●●
●
●●
●
●
●●●● ●
●●
●
●●
●●● ●●
●
●
● ●●
●
●●
● ●
●●●
●
●●
●
●
●
●
−0.15
−0.15
●
−3 −2 −1 0 1 2 3 −3 −2 −1 0 1 2 3
Q−Q plot for margin 2 (simulated data) Q−Q plot for margin 2 (real data)
●
●
●●
0.05
0.05
●
● ●●●
● ● ●●●
●●
●● ●
●●●
●
● ●
●●
●
●
●●●●●●
● ●
●●
●
●
●
●
●
●●
●
●●● ●●
●
●
●●
●
●●
●
●● ●
●●
●
●
●
●●
●
●
●●
●
●●
● ●
●
●
●
●●
●●
●
● ●
●
●●
●●
●●
●
●
●●
●
●●
●
● ●
●●
●
Sample quantiles
Sample quantiles
●
●●
●
●●
● ●
●●
●
●
●●
●
●●
●
●
●●
●
●● ●
●●
●
●
●
●●
●
●●
●
●
●●
● ●
●●
●
●●
●
●
●
●
●●
●
●
●● ●●
●
●
●●
●
●●
●
●
●●
●
●●
● ●●
●
●
●●
●
●
●
●
●●
●
●
●●
● ●
●
●●
●
●
●
●●
●
●
●●
●
● ●
●
●
●●
●
●●
●
●
●●
●
●
●●
● ●
●
●
●●
●
●
●
●
●
●
●●
●
●
●● ●
●
●
●●
●
●
●
0.00
0.00
●
●
●
●●
●
●
●●
●
●
● ●
●●
●
●
●●
●
●
●
●●
●
●
●●
●
● ●
●
●
●●
●
●●
●
●●
●
●
●●
●
● ●
●●
●
●
●●
●
●
●
●
●
●●
●
●
●●
●
● ●
●
●
●●
●
●
●●
●
●
●●
●
●●
●
●
●●
●
●
●● ●●
●
●
●●
●
●
●●
●
●
●
●
●●
●
●●
●
● ●●
●
●
●●
●
●
●●
●
●
●
●
●
●●
●
●
●●
● ●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●
●
●●
●
●
●●
●
●
● ●
●
●●
●
●
●●
●
●
●
●
●
●
●●
●
●
●● ●
●
●●
●
●
●●
●
●
●●
●
●●
●
●●
●
●
● ●
●
●●
●
●
●●
●
●
●●
●
●●
●
●
●●
● ●
●●
●
●
●●
●
●
●●
●
●●
●
●
●●
● ●
●
●
●●
●
●
●●
●
●
●
●●
●
●
●●
● ●
●
●
●●
●
●
●●
●
●
●
●
●
●●
●
●
●● ●
●●
●
●●
●
●
●
●
●
●●
●
●
●● ●
●
●●
●
●
●●
●
●●
●
●
●●
●
●
●● ●
●
●
●●
●
●
●
●
●
●●
●
●
●●
●
● ●
●
●●
●
●
●
●
●●
●
●
●●
●
●● ●●
●
●●
●
●
●
●●
●
●
●● ●
●
●
●●
●
●
●
●
●●
●
● ●
●●
●
●
●
●
●
●●
●
●●
●
●
● ●
●●
●
●●
●
●●
●
● ●
●●
●
●
●●
●
●●
●●
●
●●
●
● ●
●
●
●●
●
●●
●
●●
●
● ●
●●
●
●●
●
●
●
●
●● ●
●
●
●
●●
●
●
●●
● ●
●
●●
●
●
●
●●
●
●
●● ●
●
●
●●
●●
●●
●
● ●
●
●
−0.05
−0.05
●
●●
●● ●
●
●●● ●●
●
●●●● ●
●●
● ●● ●
● ●●
●
●●
●●
●
●
●
●●
●
●
−0.10
−0.10
●
−3 −2 −1 0 1 2 3 −3 −2 −1 0 1 2 3
● ●
●
60
60
Sample quantiles
Sample quantiles
●
●
●
40
40
●●
●
●
●●
●●●●
●●
●
20
20
●
●●
●●●●
●
● ●
●
●●
●
● ● ●
●●
●●
●
●
●●
● ●●
●●
● ●
●●
●
● ●●
●● ● ● ●●
●
●●
●
● ●
●
●●
●
●
●●●● ●
●
●
●●●●●●● ●●
●
●
●●
●●
●●
●●
●●
●●●●●● ●
●●
●
●●
●●
●
●
●
●●
●
●●
●●
●
●●
●● ●
●●
●
●●
●
●
●
●
●
●●
●
●
●●
●
●
●●
●
●●
●
●
●●
●●
●
● ●
●●
●
●●
●●
●
●●
●
●●
●●
●
●●
● ●●
●
●
●●
●
●●
●
●
●
●●
●
●●
●
●
●●
●
●●
● ●
●
●●
●
●
●●
●
●
●
●
●
●●
●
●●
●
●●
●
●
●●
●
●
●● ●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
● ●●
●
●
●●
●
●
●●
●●
●
●
●●
●
●●
●
●
●●
●
●
●●
●
●
●● ●
●
●
●●
●
●●
●
●
●●
●
●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
● ●●
●
●
●●
●
●●
●
●
●●
●
●
●
●
●
●
●●
●
●
●●
●
●
●●
●
●●
●
●
●●
●
●
● ●
●●
●
●
●●
●
●●
●
●
●●
●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
● ●●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
● ●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●
●
●●
●
●
● ●
●●
●
●
●●
●
●
●
0
0
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
● ●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
●
●
●
●●
● ●
●
●
●●
●
●
●●
●
●
●●
●
●
●●
●
0 5 10 15 0 5 10 15
0.15
0.15
● ●
●
0.10
0.10
0.10
●
● ●
●
●
● ● ●
0.05
0.05
0.05
● ●
● ●
● ● ● ●
●
● ● ● ● ●● ● ● ● ●
● ● ● ● ● ● ●●●
● ● ● ● ●
● ● ● ● ●● ●
● ●● ● ● ●
● ● ● ● ●● ● ●● ● ●
● ●● ● ● ● ●● ●● ● ● ● ●
● ● ●
● ●● ● ● ● ● ● ● ●
●●●
●● ● ● ●●● ●● ● ● ●
● ● ● ● ● ●● ●● ●
● ●● ● ●
● ●● ●● ● ●● ●
● ●● ●●● ●●● ●● ● ●●● ●●
●
●● ●●● ●● ●●● ● ● ●●●●● ● ●
● ●
●
●● ● ● ●●●● ●● ● ● ● ●● ●● ●
●●●●●● ●
●
●●●● ●
●
●
●●●● ●
●●●● ●
●● ● ●● ● ● ● ● ● ●● ●●●●
●
● ● ● ●
● ●●●●●● ●● ●●●
● ●● ●● ●● ●● ● ●●●●●● ●●●●●●● ●●
●
●●●● ● ●●● ●● ●●●
●●●● ● ● ●
●●●●● ●●● ●●●● ● ●
● ●●● ● ● ●● ●● ● ●● ● ● ●●●● ●
● ●
●●●
● ●●
● ●
●
●
●●●●●●
●
●
●● ●
●● ●●
●●
● ●●● ●●● ● ●
●● ● ● ● ● ● ●●
●
●●
SIE
●
SIE
● ●
SIE
● ● ●
●●
● ●
● ●●● ●●● ● ● ●● ●●● ●● ●● ● ● ●●● ● ●●●●● ●●
●● ● ● ●
● ● ●● ● ●
●●●●
● ●
●
●
●
●
●●●
●●
●●
●●
●●●
● ●● ●●●● ● ●● ● ●
●● ●●
●●
●●●●
● ●●
●
●●●
●
●
●
●●
●●
●
●
●
●●
●●
●
●●
●
●● ●●●●●
●
●
●
●●● ●●
● ●
● ●
● ●●
●
● ●
●●●●●●
●
●● ●
● ●
●●●●●
●●●
●● ●
●
●● ●●●
●
● ●●
● ●
●
●●
●
●●
●
●●●●
●●●
●● ●●● ●● ●● ●●● ●●
●● ●
● ●●
●●●
●●●●
●●●
●●
●●●●
●●
●
●●●
●●●●●●
●
● ●●
●●
● ●● ●●●●●● ●●
●● ● ●●●
●●●
●
●● ●
●●
● ● ●
●●
●● ●●●● ●●●
● ● ●
● ●● ●
●●
●
●●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●●
●●
●
●● ●
●●
●
●
●●
●●
●
●
●
●
●
●●
●●●
●
●● ●
●●●● ●●● ●
●●● ●● ●
● ● ●
●
●●●
●
●●
●
●
●
●
●
●
●●
●
●
●●
●
●
●●
●
●●
●●
●
●
●●
●
●●
●
●●
●
●
●
●
●
●●
●●
●
●●
●
●
●●●●
●
● ●
● ●●● ● ● ●●●●●
●
● ●
●●
●
●
●
●
●●
● ●
●●●
●
●
●
●●
●
●
●●● ●●
●
●●
●
●
●
●●
●●●● ●●●
● ● ●●● ●●●
●●●●
●
●●●●
●
●
●
●
●
●
●●
●
●●
●●
●
●
●●
●
●●
●●
●●
●
●
●
●●
●
●●
●
●●
●●
●●●●●
●● ●●●●●
●
● ●
●●●●● ●
●●
●
●●
●
●
●
●
●
●
●
●●●
●
●
●●
●
●
●●
●
●●
●
●●●
●
●
●
●
●
●
●
●
●●
●
●
●
●
●
●
●●
●●
●
●●
●
●●
●●
●
●
●●●
●●
●
●●●●
●
●●● ● ●●
●
●
● ●● ●●●● ●
●
●
●●
●●●
● ●
●●
●●
●●
●
●
●●
●
●
●
●
●
●●
●
●●
●
●
●
●
●
●
●
●
●
●●
●●
●
●●
●
●●●●
●●●● ●
● ●
●
●●●●
● ●
●
●●●●
●●
●●● ●●
●
●
● ● ● ●● ● ● ●●●●● ● ●
●● ●●
●●
● ●
● ●
● ●
●
●
●●●
●●
●●
● ●
●
●
●●●●●
●●●
●
●●
●●
● ● ●● ● ●●●●● ●
●●●
●●
●
● ●
●●
●●●
●●
●
● ●●
● ●●
●
●●
●●
● ●
● ●● ●●●●
● ●● ● ● ● ● ●●● ● ●●● ●● ●
0.00
● ● ● ●
0.00
0.00
●● ●●●●
●●●
●●●
●
●●●
●●●
●●●
●
● ●
●● ●●
● ●● ● ● ● ● ●
●●●●●
● ●●
●●
●●●
● ● ●
●
●●
● ● ●●
●●●●
● ●● ●
●● ●●●●●
● ●
●●
●●●●●●●
●●●
●●●● ●
●●●●● ●
●
●● ●●
●
●●
●
●
●
●
●
●
●●
●
●
●
●●
●●●
●
●
●●
●
●
●
●
●
●
●
●
●
●●
●
●
●●
●
●●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●●
●
●
●
●
●
●
●
●
●
●
●
●
●●
●
●●
●
●
●●
●●
●
●●●● ●●
●
● ● ●
● ●
●● ● ●●
●●
●
●●●
●
●
●●
●
●
●
●
●
●
●
●
●
●
●
●●
●
●
●●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●●
●
●
●
●●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●●
●
●
●
●
●
●
●
●
●
●●
●
●
●
●
●
●
●
●
●
●●
●
●
●
●
●
●
●●
●
●
●
●
●●
●
●
●●
●●
●
●
●●
●●●● ● ● ● ● ●
●●● ● ●
●●●
●
●●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●●
●●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●●
●
●
●
●
●●
●●
●
●
●●●● ●
●● ●
●
●●●
●
●
●●
●●
●
●●
●●
●●
●
●
●
●
●
●●
●
●●
●
●
●
●
●●
●
●
●
●
●
●
●
●
●●
●●
●
●
●●
●●
●
●●
●
●
●
●
●
●●
●●
●
●
●
●●
●
●●●
●
● ●●
●
●
●
●
●
●●●
●
●
● ●●
●● ● ●
● ●●●●● ●●●
●
●●
●
●●
●
● ●●
●
●
●
●
●
●
●
●
●
●●
●
●
●
●●
●
●
●
●
●
●●
●
●
●
●
●
●●
●
●●
●
●
●
●
●
●
●
●●
●
●
●
●
●●
●
●●●
●
●●
●
●●
●
●
●
●●
●
●●● ●● ● ●●
●●●●●
●
●
●
●
●
●
●
●●
●
●
●
●●
●●
●
●
●
●●
●
●●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●●
●●
●
●
●
●
●●
●●
●
●●
●
●
●●
●
●
●
●
● ●
●●
●
●
● ●●
●●
●●●●●● ●
● ●●●●
●●●●
●
●●●●
●●●
●
●●●
●
●
●●
●
●●
●
●●
●
●●
●
●
●●
●
●
●●
●●
●
●
●
●●
●
●●
●
●
●●
●
●●
●
●●● ●
●●●
●●
● ● ● ●● ●● ● ●●●
●●
●●
●
●●
●
●●●
●●
●●
●
●
●●
●●●●
●
●●
●
●●
●●
●
●●
●
●●●●
●
●●
●●●●
●
●●
●●
●
●●
●
●
●
●
●●●●
● ●●●●●● ● ●● ● ●● ●● ●● ●
●●●●●
●●
●
●●
●●
●
●●
●
●
●
●●
●
●●
●
●
●●
●
●●
●
●●
●
●
●●
●●
●●
●●
●
●
●●
●●
●
●
●●
●●
●●
●●
●●
●●
●●●
●●
● ●
●● ● ●
● ●
●
●●●
●
●●
●●
●●
●●
●●
●●
●●
●
●●
●●
●●
●
●●
●●
●
●
●●
●●
●
●●
●
●
●●
●●
●●
●
●
●●
●
●●
●
●●
●
●●
●
●●●●
● ●
●●
● ●●
● ●● ● ●● ● ● ● ● ●●●●●●●
●●
●
●
● ●
●●
●
●●
●
●●
●
●●●●
●
●
●
●
●●
●●
●●
●
●●●
●
●●●
●
●●
●
●
●
●
●●
●
●●●
●
●●
●●●●
●●●● ●● ● ●●●
●●
● ●●●
●●●
●
●●
●●
●
●●
●
●●
●
●
●
●
●●
●
●
●●
●
●
●
●
●●
●
●
●
●
●
●
●●
●
●
●
●
●
●
●
●●
●
●
●●
●
●●
●
●●
●●
●
●●
●
●●●●
●●●● ●● ● ●● ● ●
●
●● ●● ●● ●
●
●●
●●
●
●
●
●
●●
●
●
●●●
●●
●●
●●
●
●
●
●
●●
●
●
●
●
●
●
●
●
●
●●
●
●
●
●
●
●
●
●
●
●●
●●
●
●
●
●
●
●
●
●
●
●
●
●●
●
●●
●
●
●
●●
●
●●●
●
●
●● ●●
●●●
●●● ● ● ●● ●
●
●●
●●●
●
●
●●
●
●
●●
●
●
●●
●●
●
●●
●
●
●●
●
●●
●
●
●
●
●
●●
●
●●●
●
●
●●
●
●
●
●●
●
●
●
●
●
●●
●
●
●●
●
●●●
●
●●
●
●
●
●●
●
●
●●
●
●●
●
●●
●
●
●●
●
●
●●
●
●
●●●●●●● ●
● ●● ● ● ●●
●● ●●
●●
●
●
●●
●
●●
●
●
●
●●
●
●
●
●●
●
●
●●
●
●
●
●●
●●
●
●
●
●●
●
●●
●
●
●
●●
●
●●
●
●
●●
●
●
●
●
●
●
●
●
●●
●●
●
●●
●●
●
●●
●
● ●● ●● ● ●
●●●●●●
●● ●
●●
●●●
●
●
●
●●
●
●
●●●
●
●●
●
●●
●●
●
●
●●
●●●
●●
●●●
●
●●●
●● ●●●● ●● ● ●
●●●●●●● ● ●
●
●●●
●
●
●
●●●●
●
●
●●
●●
●
●●
●
●
●●
●●
●●
●●
●
●●
●●
●
●
●●
●
●
● ●
●●
●
●●
●
● ●
●
●●●
● ●
● ●
●●●●● ●● ● ● ● ●● ●●●●● ●●
●
●●
●●
●
●
●●
●●
●
●
●●
●
●
●
●●
●
●
●
●
●●
●
●
●
●●
●
●●
●●
●●
●
●
●
●●●
●
●
●●
●
●
● ●
●
●● ●
●
● ●●●
●
●● ● ●● ●●● ●●
●●
●●
●
●●
●
●
●●
●●●
●
●●
●
●●
●●
●
●
●
●●
●
●●
●
●
●
●●
●
●
●
●●
●
●●
●●
●
●●
● ●●
●
●●
●
●●● ● ● ● ● ● ● ●
●
●
●●
●●●
●
●●
●●●
●
●●
●●
●
●
●
●●● ●
●
●●
●
●●
●●
●●
●
●
●
●
●
●●
●
●●
●
●●
●
●
●
●●
●
●●●
●
●
●●
●
●●●
●●●●
●
●●●
●
●●●● ● ● ●● ● ●
●
●●●
●●
●●
●●
●
●●
●
●
●
●●
●
●●
●●
●
●
●●
●●
●
●
●●
●●
●●
●
●
●●
●●
●
●●
●
●●●
●●● ●
● ●
●●● ● ● ●
● ●●● ●●
●
●●
●
●
●
●
●
●●
●●
●
●
●
●
●
●
●●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●●
●
●
●
●
●
●
●
●
●
●
●
●
●
●●
●
●●
●
●
●
●●●
●●
●
●●●●
● ●
●
●
● ●
●●●● ● ●●
●●
● ● ●●●●●
●
●
●●●
●
●●
●●
●
●●
●
●●
●
●
●
●
●
●
●
●●
●
●
●
●
●
●
●
●
●
●
●
●●
●
●
●
●
●
●
●
●
●
●
●
●
●●
●
●
●●●
●
●
●
●
●
●●
●
●●●
●
●
●
●
●
●●
●
●●●●
●
●
●
●
●●●
●
●●
●●●
●
●●
●●
● ●
●● ●
●
●●
● ●
●
●
●
●
●
●●
●
●
●●
●
●●
●
●●
●
●
●
●
●
●●●
●
●
●
●
●
●
●●
●
●
●●
●
●
●
●
●
●●
●
●●
●
●●●
●
●●
●
●●●
●●
●●●
● ●●● ●●●●●●●
●● ●
●●
●●
●●
●
●●
●●
●
●
●
●●
●
●
●
●
●●
●
●●
●●
●
●
●
●
●
●
●●
●●
●
●
●●●
●●
●●●
● ●
●
● ●● ●●●
●●
●●●●
●
●
●●●
●
●
●
●●
●●
●
●
●
●
●
●
●●
●
●
●
●●
●
●
●
●
●
●
●
●
●●●●
●
●
●
●
●
●
●
●●●
●
●●
●
●●●
●
●
●●
●●
●●● ●●
●
● ● ● ●● ● ●●
●
●
●● ●●●●
●●
●●
●
●
●
●●●
●
●●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●
●●
●
●
●
●
●
●
●
●●
●
●
●●
●●
●
●●●●● ●●●
●● ● ●
●
●
● ●●
●● ●●
●●
●
●●●●
●
●●●
●●●
●●
●●●●●● ● ● ● ●● ● ●
●●
● ●
●
●●●
●
●●
●●
●●
●●
●●
● ●●
●●
● ●●
●● ●●● ●
● ● ●● ● ●
● ●
●●●●●
●●
● ●
●
●●●
●
●
●
●●●●
● ●●
●
●●●●● ●●
● ●● ●● ●
●● ●
●●
● ●● ●
●●
● ●
●
●●●
●●
● ●●●
●●
● ●● ● ●
● ●●●●●●
●
● ●
●●●
●●●●
●
●●
●
● ●
●
●●
●
●●
●●
●●●●●
●●●
●●●● ●
●● ●
● ● ●●
● ●●●● ● ●●
●●
●●●
●
●
●●●
●●●●●●
● ●
●●●
●●●●●● ●
● ●●
●
● ● ●● ●●●●●
●●
●
●●
●●
●●
● ●●
●●●
● ● ● ● ●●●●●●●●●●
●
●
●
●
●
●
●
●
●●
●●
●
●●
●●●
●
●
●
●
●●
●●
●●
●●
●●
●
●●●●
●
●
●
●●●●
●
●●●●
● ● ●●●
●●●●● ●
●
● ●●
●
●●●
●●●●
● ●
●●
● ●●● ● ●
●●
●● ●
● ● ●●● ●
●
●●
●
●●●●
●●
●●●●●●
●
●● ● ● ●●●●●●●
● ●
●●
●●●
●●
●
●●
●
●
● ●
●
●
●
●
●●●
●
●
●●●
●●
●●●●●●●
●●
●
●
●●●● ●● ●
● ● ● ●●●●●●●
● ●
●
●
●●
●
● ●
●●●●
●●
●
● ●●
●
●
● ● ● ● ●
● ●● ●●● ●
●● ●●●●●● ●● ● ● ●●
●●
●
● ●●
●●
●
●
●●
●
●●●●●
● ●
● ●
●
●●
●●
●●●
● ●● ● ●●● ●●●●●●● ● ●●
●
●●
●●
● ●
● ● ●
● ●
●● ●●● ● ●● ●●●●●●
●●●
●● ● ●●● ● ● ●● ●●
●● ●
● ●
●● ●●●
● ●●
● ● ● ●●● ●●● ●
●● ●● ●
●● ●
● ●
●●●●● ●●
●●●●●
●● ●●● ● ● ● ●●
● ● ●● ● ●●●●
●
● ●●
● ●
●
●●●● ● ●●
● ● ●● ●● ●●●● ●
●●
●●
● ●
● ●
● ●● ●●●● ● ●●●●●● ●●
●●●●●●
●
● ●●●
● ● ●● ●● ●● ●
●● ● ● ● ●● ● ● ●● ● ● ●● ● ● ●
● ● ●● ● ●● ● ● ● ●● ●
● ●●●
● ● ● ● ● ●● ●
●
●●
●● ●● ● ●●● ● ●
● ●● ● ● ● ● ● ●
●●
●● ●
−0.05
−0.05
−0.05
● ●
● ● ●●
● ● ● ● ● ● ● ● ●● ● ● ●
● ●
● ● ● ●● ●
● ●●
● ● ● ●● ● ● ●
● ● ●
● ● ● ●
● ●
● ● ● ●
● ● ●
●
● ● ●
● ●
●
●
● ●
● ●
● ●
●
−0.10
−0.10
−0.10
●
●
−0.15 −0.10 −0.05 0.00 0.05 0.10 0.15 0.20 −0.15 −0.10 −0.05 0.00 0.05 0.10 0.15 0.20 −0.15 −0.10 −0.05 0.00 0.05 0.10 0.15 0.20
See the appendix for proofs for Theorems 6.15 and 6.16.
If Y has a density fY , it satisfies fY (y) = g(kyk2 ) for a function
g : [0, ∞) → [0, ∞) referred to as density generator (i.e. fY is constant
on spheres); see the appendix for a proof.
© QRM Tutorial Section 6.3.1
Corollary 6.17
d
If Y ∼ Sd (ψ) and P(Y = 0) = 0, then (kY k, kYY k ) = (R, S) since
d RS RS
(kY k, kYY k ) = (kRSk, kRSk ) = (|R|kSk, |R|kSk ) = (R, S).
For (elliptically distributed; see soon) X = µ+AY with E(R2 ) < ∞ and
2)
Cholesky factor A of a covariance matrix Σ, we have cov X = E(R d Σ
and corr X = P (the correlation matrix corresponding to Σ).
© QRM Tutorial Section 6.3.1
Example 6.20 (t distribution)
For Y ∼ td (ν, 0, Id ), R2 = Y ′ Y √
= W Z ′ Z for Z ∼ Nd (0, Id ). Thus
Cor.6.17 Y = WZ
R2 Z ′ Z/d χ2d /d
= = ∼ F (d, ν)
d (ν/W )/ν χ2ν /ν
ν
and thus E(R2 /d) = ν−2 .
This, together with Example 6.19, implies that X ∼ td (ν, µ, Σ) has
ν
cov X = ν−2 Σ and corr X = P (which we already know from Sec-
tion 6.2.1); note that in the univariate case X ∼ t(ν, µ, σ 2 ) and
ν
var(X) = ν−2 σ2.
We also see that we can use a Q-Q plot of the order statistics of
R2 /d = kY k2 /d versus the theoretical quantiles of a (hypothesized)
F (d, ν) distribution to check the goodness-of-fit of the hypothesized t
distribution (in any dimensions).
See the appendix for the form of the density generator g.
© QRM Tutorial Section 6.3.1
Example 6.21 (Understanding spherical distributions)
p
n = 500 realizations of S (left) and Y = RS (right) for R ∼ dF (d, ν),
d = 2, ν = 4 (as for the multivariate t distribution with ν = 4).
10
1.0
●●
●●●● ● ● ●●
●
● ●● ●
● ●●● ●
●●
● ●●●● ●●●
●●●
● ●●●●
● ●
●●●
● ●
●●
●
●●●
●●● ●●
●●
●
●●
●●
●●●●●
●●
●●
●●
●●
●
● ●
● ●●
●●●
●● ●
●
●● ●
●●
●● ●●
●
●●●
●
●
● ●●
●●
●● ●
●
●
●●● ●
●●
●●● ●
●
●● ●
0.5
●●
5
●
● ●
●
● ●
●
●● ●
● ●● ●
●
● ●● ●
●
● ● ●
● ●
●
●
●
●
● ● ●
● ●
● ● ●
● ●
●
● ● ●
● ● ●
●
● ● ● ●
●
● ● ● ●● ●
● ● ● ●
● ● ● ● ● ●● ● ●
● ● ● ● ●
● ● ● ● ● ●
●
●
● ●
● ● ● ● ●● ● ● ●
● ●● ● ● ●
● ● ●● ● ● ●
● ● ●
● ● ●● ● ●●●●● ●● ●
●● ● ●●
●●
●●●●● ● ● ● ●
● ● ●● ●● ●●
● ●●●
● ●●● ● ●
● ● ● ● ● ●●●● ● ● ●●●●●●
● ●
● ●● ●● ●●●
●● ● ●●● ●
●
●●●
●●●●●● ●●
● ●● ● ●●●●● ●●
●●●
0.0
Y2
●●● ●●● ● ● ●●●
● ● ●●●●● ●●
●● ● ●●
● ● ●● ●●●
●●● ●
0
●
● ● ●
● ●●●
● ●●●●●
●● ●●
●●
●
●●● ● ●● ●●● ● ●
● ●● ● ● ●● ●●
● ●●● ●●● ●
● ● ● ●●●●● ●●●●●●
●
●●
●●●
●●
●●●
●
●
●●
●
● ● ●●●
●●
●
●● ●
● ●
● ● ● ● ●●●●● ●● ●●●●●●●●●● ●●● ●
● ● ●●
● ● ● ●
●
●●●● ●
●● ●● ●● ● ● ●●●●●
● ● ● ● ●● ●
●● ● ●● ●●
●● ● ●
● ●● ● ● ●● ●● ● ●
●
● ● ●●●
●
●●●●
●● ●●●●● ●● ●●
●●●●●●● ● ●
● ● ● ●●●●● ●● ● ● ●
●
●
● ● ●●●●● ●● ● ● ●● ● ●
● ●
● ● ● ● ●● ● ●
● ● ● ●● ●● ●● ●
● ● ● ● ● ●
●
● ● ● ● ●●
● ●
● ● ● ● ●
●
● ● ●
● ●
● ● ●
● ● ●
● ● ● ●
●
● ●
● ● ●
● ●
●
●
●
● ●
−0.5
● ●
● ●
●●
● ●
−5
● ●
●
●● ●
●● ●
●
● ●
●●
●● ●
●
● ●
●
●● ●
● ●●
● ●
●
●●● ●●
●
●
●
●
● ●●
●
●
● ●●
●● ●
●●●
●
●
●
●
●●● ●●
●
●● ●
●
●● ●
●
●● ●
●●
●● ●●
●●●
●
●
●●● ●●●
−1.0
●●
●● ●
●
● ● ●
●●●●
●●●
●
●●● ● ●● ●
●●●●● ●● −10
●●● ● ●
●● ● ●●
S1 Y1
10
●●
●
●
● ● ● ●
● ●
●
● ● ●
● ●● ● ●
● ●
●●
● ● ●●● ●
● ● ● ●● ● ●●
●● ● ● ● ●● ●
● ●●● ● ● ●
● ● ● ● ●●● ● ●
●● ● ● ●
● ● ● ●● ●
●
● ● ● ● ●
●
● ● ●●●●
● ●●●●●●●●●● ●●● ●●
● ●●●● ● ●●● ● ●●●● ●●
● ● ●
● ●● ●● ●● ●
●
●
● ●
●●●●● ●● ●●● ●● ●
● ● ●●●●●● ●●● ●●
● ●●
● ●●
●●●
●● ●● ●
●
●●
● ●● ●
●●●
● ●●●
● ●●●
●●●
● ● ● ●● ●●●● ●● ●●●
●● ● ● ●●●● ● ●
●
● ● ● ● ●●● ●●
●●●
●
●
●●
●● ●●●
●● ●●●●
●● ● ●
● ●
●●●●●
●●
●● ● ●
●●
● ●
●●●
●
●
●●
●
●
●
●●● ● ●
●●● ● ●
●●●●●
●●
●●
●
●●● ●
●●
●●●● ●
● ●
5
5
● ● ●●
●●
●●●
● ● ● ● ●● ● ●
● ●● ● ● ● ●
●● ● ● ●
●● ●● ●●●● ●
●●
●● ●●●● ● ●
●● ● ●
● ●● ●● ●● ●
● ●
●●●●
●●
●
● ● ●●● ●●●●● ● ● ● ●
●● ●● ●● ●● ●●
●● ● ● ● ●●
● ● ● ● ● ●●
●●●● ●
●●● ● ● ● ●
● ● ● ●●
● ● ● ●●●● ●
● ●
●●●
● ● ●
● ● ● ● ●
● ● ●
●
● ● ● ● ●● ● ●
● ●● ● ●
● ●
●● ●
● ● ●●● ●
● ● ● ●● ● ●●
●● ● ● ● ●● ●
● ●●● ● ● ● ●
● ● ● ● ●●● ● ●
● ●
●● ● ●●●●● ●●● ● ● ●
● ● ● ●
●●●● ●●● ●●
● ●
●●●●
●●●●●●
● ●●●●●●
● ● ●● ●
● ●● ●●●
●● ●● ● ● ●
●
●
●●●
●●
●●●● ● ●● ●● ●●
●
●
●
●
● ●●● ●● ●●● ●●●● ● ●●
●● ●● ●
● ● ● ●●
●
●●●●●●●●
●●
●● ● ●●
X2
X2
● ●●● ● ●
● ● ● ● ●● ●●●● ●● ●●●
●● ● ● ●●●● ● ●
●
● ● ● ●●●
0
0
●●●●● ●● ● ●● ● ●
● ● ●
●●●●●
●●
●● ●●
●●
●●●●
●
●●
●
●
●
●
●
●
●
●
●●
●
●● ●●●●● ●
● ●● ●
●●
●
●●
●●
● ●
● ● ●● ●
●
● ●● ●● ●
●●
●
●●●
●●●
●
●●●●● ●● ●
●
●
●● ● ●● ●● ●
●
●
●●●●●
●
●
●
●●●●●●●●●
●●●● ●●
● ●
●●● ● ●●
● ● ●
●●●●● ● ● ●
● ●
●●●● ● ● ● ●
●● ●● ● ● ● ●
●● ●● ●●
● ● ● ● ●●
● ● ● ● ● ●●
●●●● ●
●●● ● ● ● ● ●
● ● ●●
●
● ● ●●● ●
● ●
●●●
● ● ●
● ●
● ●● ● ●
●
●
● ●
● ●
●
−5
−5
●
−10
●
−10
−10 −5 0 5 10 −10 −5 0 5 10
X1 X1
Univariate regression
Consider the (univariate) time series regression model
Xt,j = aj + b′j Ft + εt,j , t ∈ {1, . . . , n}.