Chapter 2

Multivariate Normal Distribution

[37/97%]
Show
LaTeX
Problem 2.1

Suppose that ε∼N3(0,σ2I3)\varepsilon \sim N_{3}\left(0, \sigma^{2} \mathrm{I}_{3}\right) and that Y0Y_{0} is N(0,σ02)N\left(0, \sigma_{0}^{2}\right), independently of the εi\varepsilon_{i} 's. Define

Yi=ρYi−1+εi(i=1,2,3) Y_{i} = \rho Y_{i-1}+\varepsilon _{i} \quad (i = 1,2,3)
?
(a)

Find the variance-covariance matrix of Y=(Y1,Y2,Y3)T\mathbf{Y} = \left(Y_{1}, Y_{2}, Y_{3}\right)^{T}.

(b)

What is the distribution of Y\mathbf{Y}?

Problem 2.2

Let Y∼Nn(0,In)\mathbf{Y} \sim N_{n}\left(\mathbf{0}, \mathbf{I}_{n}\right), and put X=AY,U=BY\mathbf{X} = \mathbf{A Y}, \mathbf{U} = \mathbf{B Y} and V=CY\mathbf{V} = \mathbf{C Y}. Suppose that Cov⁡[X,U]=0\operatorname {Cov}\left[ \mathbf{X}, \mathbf{U} \right] = \mathbf{0} and Cov⁡[X,V]=0\operatorname {Cov}\left[ \mathbf{X}, \mathbf{V} \right] = \mathbf{0}. Show that X\mathbf{X} is independent of U+V\mathbf{U}+\mathrm{V}.

?
Problem 2.3

If Y1,Y2,…,YnY_{1}, Y_{2}, \ldots , Y_{n} is a random sample from N(μ,σ2)N\left(\mu , \sigma^{2}\right), prove that Yˉ\bar{Y} is independent of ∑i=1n−1(Yi−Yi+1)2\sum_{i = 1}^{n-1}\left(Y_{i}-Y_{i+1}\right)^{2}.

?
Problem 2.4

If X\mathbf{X} and Y\mathbf{Y} are nn-dimensional vectors with independent multivariate normal distributions, prove that aX+bYa \mathbf{X}+b \mathbf{Y} is also multivariate normal.

?
Problem 2.5

If Y∼Nn(0,In)\mathbf{Y} \sim N_{n}\left(\mathbf{0}, \mathbf{I}_{n}\right) and a\mathbf{a} is a nonzero vector, show that the conditional distribution of YTY\mathbf{Y}^{T} \mathbf{Y} given aTY=0\mathbf{a}^{T} \mathbf{Y} = 0 is χn−12\chi_{n-1}^{2}.

?
Problem 2.6

Let Y∼Nn(μ1n,Σ)\mathbf{Y} \sim N_{n}\left(\mu \mathbf{1}_{n}, \mathbf{\Sigma }\right), where Σ=(1−ρ)In+ρ1n1nT\mathbf{\Sigma } = (1-\rho ) \mathbf{I}_{n}+\rho \mathbf{1}_{n} \mathbf{1}_{n}^{T} and ρ>−1/(n−1)\rho > -1 /(n-1). Show that ∑i(Yi−Yˉ)2/(1−ρ)\sum_{i}\left(Y_{i}-\bar{Y}\right)^{2} /(1-\rho ) is χn−12\chi_{n-1}^{2}.

?
Problem 2.7

Let Yi,i=1,…,n\mathbf{Y}_{i}, i = 1, \ldots , n, be independent Np(μ,Σ)N_{p}(\mathbf{\mu }, \mathbf{\Sigma }) random vectors. Show that

1n−1∑i=1n(Yi−Y‾)(Yi−Y‾)T \frac{1}{n-1} \sum _{i = 1}^{n}\left(\mathbf{Y}_{i}-\overline{\mathbf{Y}}\right)\left(\mathbf{Y}_{i}-\overline{\mathbf{Y}}\right)^{T}

is an unbiased estimate of Σ\mathbf{\Sigma }.

?
Problem 2.8

Let Y∼Nn(0,In)\mathbf{Y} \sim N_{n}\left(\mathbf{0}, \mathbf{I}_{n}\right) and let A\mathbf{A} and B\mathbf{B} be symmetric idempotent matrices with AB=BA=0\mathbf{A B} = \mathbf{B A} = \mathbf{0}. Show that YTAY,YTBY\mathbf{Y}^{T} \mathbf{A Y}, \mathbf{Y}^{T} \mathbf{B Y} and YT(In−A−B)Y\mathbf{Y}^{T}\left(\mathbf{I}_{n}-\mathbf{A}-\mathbf{B}\right) \mathbf{Y} have independent chi-square distributions.

?
Problem 2.9

Let (Xi,Yi),i=1,2,…,n\left(X_{i}, Y_{i}\right), i = 1,2, \ldots , n, be a random sample from a bivariate normal distribution, with means μ1\mu_{1} and μ2\mu_{2}, variances σ12\sigma_{1}^{2} and σ22\sigma_{2}^{2}, and correlation ρ\rho, and let

W=(XT,YT)T=(X1,X2,…,Xn,Y1,Y2,…,Yn)T \mathbf{W} = \left(\mathbf{X}^{T}, \mathbf{Y}^{T}\right)^{T} = \left(X_{1}, X_{2}, \ldots , X_{n}, Y_{1}, Y_{2}, \ldots , Y_{n}\right)^{T}
?
(a)

Show that W\mathbf{W} has a N2n(μ,Σ)N_{2 n}(\mathbf{\mu }, \mathbf{\Sigma }) distribution, where

μ=(μ11nT,μ21nT)T and Σ=[σ12Inρσ1σ2Inρσ1σ2Inσ22In] \mathbf{\mu } = \left(\mu _{1} \mathbf{1}_{n}^{T}, \mu _{2} \mathbf{1}_{n}^{T}\right)^{T} \quad \text{ and } \quad \mathbf{\Sigma } =\left[\begin{smallmatrix} \sigma _{1}^{2} \mathbf{I}_{n} & \rho \sigma _{1} \sigma _{2} \mathbf{I}_{n} \\ \rho \sigma _{1} \sigma _{2} \mathbf{I}_{n} & \sigma _{2}^{2} \mathbf{I}_{n} \end{smallmatrix}\right]
(b)

Find the conditional distribution of X\mathbf{X} given Y\mathbf{Y}.

Problem 2.10

If Y∼N2(0,Σ)\mathbf{Y} \sim N_{2}(\mathbf{0}, \mathbf{\Sigma }), where Σ=(σij)\mathbf{\Sigma } = \left(\sigma_{i j}\right), prove that

(YTΣ−1Y−Y12σ11)∼χ12 \left(\mathbf{Y}^{T} \mathbf{\Sigma }^{-1} \mathbf{Y}-\frac{Y_{1}^{2}}{\sigma _{11}}\right) \sim \chi _{1}^{2}
?
Problem 2.11

Let a0,a1,…,ana_{0}, a_{1}, \ldots , a_{n} be independent N(0,σ2)N\left(0, \sigma^{2}\right) random variables and define

Yi=ai+ϕai−1(i=1,2,…,n) Y_{i} = a_{i}+\phi a_{i-1} \quad (i = 1,2, \ldots , n)

Show that Y=(Y1,Y2,…,Yn)T\mathbf{Y} = \left(Y_{1}, Y_{2}, \ldots , Y_{n}\right)^{T} has a multivariate normal distribution and find its variance-covariance matrix. (The sequence Y1,Y2,…Y_{1}, Y_{2}, \ldots is called a moving average process of order one and is a commonly used model in time series analysis.)

?
Problem 2.12

Suppose that Y∼N3(0,In)\mathbf{Y} \sim N_{3}\left(\mathbf{0}, \mathbf{I}_{n}\right). Find the m.g.f. of 2(Y1Y2−Y2Y3−Y3Y1)2\left(Y_{1} Y_{2}-Y_{2} Y_{3}-Y_{3} Y_{1}\right). Hence show that this random variable has the same distribution as that of 2U1−U2−U32 U_{1}-U_{2}-U_{3}, where the UiU_{i} 's are independent χ12\chi_{1}^{2} random variables.

?
Problem 2.13

Theorem 2.3 can be used as a definition of the multivariate normal distribution. If so, deduce from this definition the following results:

?
(a)

If Z1,Z2,…,ZnZ_{1}, Z_{2}, \ldots , Z_{n} are i.i.d. N(0,1)N(0,1), then Z∼Nn(0,In)\mathbf{Z} \sim N_{n}\left(\mathbf{0}, \mathbf{I}_{n}\right).

(b)

If Y∼Nn(μ,Σ)\mathbf{Y} \sim N_{n}(\mathbf{\mu }, \mathbf{\Sigma }), then Y\mathbf{Y} has m.g.f. (2.5).

(c)

If Σ\mathbf{\Sigma } is positive-definite, prove that Y\mathbf{Y} has density function (2.1).

Problem 2.14

Let Y=(Y1,Y2,…,Yn)T\mathbf{Y} = \left(Y_{1}, Y_{2}, \ldots , Y_{n}\right)^{T} be a vector of nn random variables (n>3)(n > 3) with density function

f(y)=(2π)−n/2exp⁡(−12∑i=1nyi2){1+∏i=1n[yiexp⁡(−12yi2)]}−∞<yi<∞(i=1,2,…,n) \begin{aligned} f(\mathbf{y}) = & (2 \pi )^{-n / 2} \exp \left(-\frac{1}{2} \sum _{i = 1}^{n} y_{i}^{2}\right)\left\{ 1+\prod _{i = 1}^{n}\left[y_{i} \exp \left(-\frac{1}{2} y_{i}^{2}\right)\right]\right\} \\ & -\infty < y_{i} < \infty \quad (i = 1,2, \ldots , n) \end{aligned}

Prove that any subset of n−1n-1 random variables are mutually independent N(0,1)N(0,1) variables.

?
Problem 2.15

Suppose that Y=(Y1,Y2,Y3,Y4)T∼N4(0,I4)\mathbf{Y} = \left(Y_{1}, Y_{2}, Y_{3}, Y_{4}\right)^{T} \sim N_{4}\left(\mathbf{0}, \mathrm{I}_{4}\right), and let Q=Y1Y2−Y3Y4Q = Y_{1} Y_{2}-Y_{3} Y_{4}.

?
(a)

Prove that QQ does not have a chi-square distribution.

(b)

Find the m.g.f. of QQ.

Problem 2.16

If Y∼Nn(0,In)\mathbf{Y} \sim N_{n}\left(\mathbf{0}, \mathbf{I}_{n}\right), find the variance of

(Y1−Y2)2+(Y2−Y3)2+⋯+(Yn−1−Yn)2 \left(Y_{1}-Y_{2}\right)^{2}+\left(Y_{2}-Y_{3}\right)^{2}+\cdots +\left(Y_{n-1}-Y_{n}\right)^{2}
?
Problem 2.17

Given Y∼Nn(μ,Σ)\mathbf{Y} \sim N_{n}(\mathbf{\mu }, \mathbf{\Sigma }), prove that

Var⁡[YTAY]=2tr⁡(AΣAΣ)+4μTAΣAμ. \operatorname {Var}\left[\mathbf{Y}^{T} \mathbf{A} \mathbf{Y}\right] = 2 \operatorname {tr}(\mathbf{A} \mathbf{\Sigma } \mathbf{A} \mathbf{\Sigma })+4 \mathbf{\mu }^{T} \mathbf{A} \mathbf{\Sigma } \mathbf{A} \mathbf{\mu } .
?
Problem 2a.1

Show that

f(y1,y2)=k−1exp⁡[−12(2y12+y22+2y1y2−22y1−14y2+65)] f\left(y_{1}, y_{2}\right) = k^{-1} \exp \left[-\frac{1}{2}\left(2 y_{1}^{2}+y_{2}^{2}+2 y_{1} y_{2}-22 y_{1}-14 y_{2}+65\right)\right]

is the density of a bivariate normal random vector Y=(Y1,Y2)T\mathbf{Y} = \left(Y_{1}, Y_{2}\right)^{T}.

?
(a)

Find kk.

(b)

Find E[Y]\mathbb {E}\left[\mathbf{Y}\right] and Var⁡[Y]\operatorname {Var}\left[\mathbf{Y}\right].

Problem 2a.2

Let U\mathbf{U} have density gg and let Y=A(U+c)\mathbf{Y} = \mathbf{A}(\mathbf{U}+\mathbf{c}), where A\mathbf{A} is nonsingular. Show that the density ff of Y\mathbf{Y} satisfies

f(y)=g(u)/∣det⁡(A)∣ f(\mathbf{y}) = g(\mathbf{u}) /\left|\operatorname {det}\left(\mathbf{A}\right)\right|

where y=A(u+c)\mathbf{y} = \mathbf{A}(\mathbf{u}+\mathbf{c}).

?
Problem 2a.3
?
(a)

Show that the 3×33 \times 3 matrix

Σ=[1ρρρ1ρρρ1] \mathbf{\Sigma } =\left[\begin{smallmatrix} 1 & \rho & \rho \\ \rho & 1 & \rho \\ \rho & \rho & 1 \end{smallmatrix}\right]

is positive-definite for ρ>−12\rho > -\frac{1}{2}.

(b)

Find Σ1/2\mathbf{\Sigma }^{1 / 2} when

Σ=[1ρρ1] \mathbf{\Sigma } =\left[\begin{smallmatrix} 1 & \rho \\ \rho & 1 \end{smallmatrix}\right]
Problem 2b.1

Find the moment generating function of the bivariate normal distribution given in Example 2.2.

?
Problem 2b.2

If Y∼Nn(μ,Σ)\mathbf{Y} \sim N_{n}(\mathbf{\mu }, \mathbf{\Sigma }), show that Yi∼N(μi,σii)Y_{i} \sim N\left(\mu_{i}, \sigma_{i i}\right).

?
Problem 2b.3

Suppose that Y∼N3(μ,Σ)\mathbf{Y} \sim N_{3}(\mathbf{\mu }, \mathbf{\Sigma }), where

μ=[212] and Σ=[211130101] \mathbf{\mu } =\left[\begin{smallmatrix} 2 \\ 1 \\ 2 \end{smallmatrix}\right] \quad \text{ and } \quad \mathbf{\Sigma } =\left[\begin{smallmatrix} 2 & 1 & 1 \\ 1 & 3 & 0 \\ 1 & 0 & 1 \end{smallmatrix}\right]

Find the joint distribution of Z1=Y1+Y2+Y3Z_{1} = Y_{1}+Y_{2}+Y_{3} and Z2=Y1−Y2Z_{2} = Y_{1}-Y_{2}.

?
Problem 2b.4

Given Y∼Nn(μ,In)\mathbf{Y} \sim N_{n}\left(\mathbf{\mu }, \mathbf{I}_{n}\right), find the joint density of aTY\mathbf{a}^{T} \mathbf{Y} and bTY\mathbf{b}^{T} \mathbf{Y}, where aTb=0\mathbf{a}^{T} \mathbf{b} = 0, and hence show that aTY\mathbf{a}^{T} \mathbf{Y} and bTY\mathbf{b}^{T} \mathbf{Y} are independent.

?
Problem 2b.5

Let (Xi,Yi),i=1,2,…,n\left(X_{i}, Y_{i}\right), i = 1,2, \ldots , n, be a random sample from a bivariate normal distribution. Find the joint distribution of (Xˉ,Yˉ)(\bar{X}, \bar{Y}).

?
Problem 2b.6

If Y1Y_{1} and Y2Y_{2} are random variables such that Y1+Y2Y_{1}+Y_{2} and Y1−Y2Y_{1}-Y_{2} are independent N(0,1)N(0,1) random variables, show that Y1Y_{1} and Y2Y_{2} have a bivariate normal distribution. Find the mean and variance matrix of Y=(Y1,Y2)T\mathbf{Y} = \left(Y_{1}, Y_{2}\right)^{T}.

?
Problem 2b.7

Let X1X_{1} and X2X_{2} have joint density

f(x1,x2)=12πexp⁡[−12(x12+x22)][1−x1x2(1+x12)(1+x22)]−∞<x1,x2<∞. \begin{array}{r} f\left(x_{1}, x_{2}\right) = \frac{1}{2 \pi } \exp \left[-\frac{1}{2}\left(x_{1}^{2}+x_{2}^{2}\right)\right]\left[1-\frac{x_{1} x_{2}}{\left(1+x_{1}^{2}\right)\left(1+x_{2}^{2}\right)}\right] \\ -\infty < x_{1}, x_{2} < \infty . \end{array}

Show that X1X_{1} and X2X_{2} have N(0,1)N(0,1) marginal distributions.

?
Problem 2b.8

Suppose that Y1,Y2,…,YnY_{1}, Y_{2}, \ldots , Y_{n} are independently distributed as N(0,1)N(0,1). Calculate the m.g.f: of the random vector

(Yˉ,Y1−Yˉ,Y2−Yˉ,…,Yn−Yˉ) \left(\bar{Y}, Y_{1}-\bar{Y}, Y_{2}-\bar{Y}, \ldots , Y_{n}-\bar{Y}\right)

and hence show that Yˉ\bar{Y} is independent of ∑i(Yi−Yˉ)2\sum_{i}\left(Y_{i}-\bar{Y}\right)^{2}.

?
Problem 2b.9

Let X1,X2X_{1}, X_{2}, and X3X_{3} be i.i.d. N(0,1)N(0,1). Let

Y1=(X1+X2+X3)/3Y2=(X1−X2)/2Y3=(X1+X2−2X3)/6 \begin{aligned} & Y_{1} = \left(X_{1}+X_{2}+X_{3}\right) / \sqrt{3} \\ & Y_{2} = \left(X_{1}-X_{2}\right) / \sqrt{2} \\ & Y_{3} = \left(X_{1}+X_{2}-2 X_{3}\right) / \sqrt{6} \end{aligned}

Show that Y1,Y2Y_{1}, Y_{2} and Y3Y_{3} are i.i.d. N(0,1)N(0,1). (The transformation above is a special case of the so-called Helmert transformation.)

?
Problem 2c.1

If Y1,Y2,…,YnY_{1}, Y_{2}, \ldots , Y_{n} have a multivariate normal distribution and are pairwise independent, are they mutually independent?

?
Problem 2c.2

Let Y∼Nn(μ1n,Σ)\mathbf{Y} \sim N_{n}\left(\mu \mathbf{1}_{n}, \mathbf{\Sigma }\right), where Σ=(1−ρ)In+ρJn\mathbf{\Sigma } = (1-\rho ) \mathbf{I}_{n}+\rho \mathbf{J}_{n} and ρ>−1/(n−1)\rho > -1 /(n-1). When ρ=0,Yˉ\rho = 0, \bar{Y} and ∑i(Yi−Yˉ)2\sum_{i}\left(Y_{i}-\bar{Y}\right)^{2} are independent, by Example 2.8. Are they independent when ρ≠0\rho \neq 0?

?
Problem 2c.3

Given Y∼N3(μ,Σ)\mathbf{Y} \sim N_{3}(\mathbf{\mu }, \mathbf{\Sigma }), where

Σ=σ2[1ρ0ρ1ρ0ρ1] \mathbf{\Sigma } = \sigma ^{2}\left[\begin{smallmatrix} 1 & \rho & 0 \\ \rho & 1 & \rho \\ 0 & \rho & 1 \end{smallmatrix}\right]

for what value(s) of ρ\rho are Y1+Y2+Y3Y_{1}+Y_{2}+Y_{3} and Y1−Y2−Y3Y_{1}-Y_{2}-Y_{3} statistically independent?

?
Problem 2d.1

Show that the m.g.f. for (2.10) is ∏1n(1−2tdi)−1/2\prod_{1}^{n}\left(1-2 t d_{i}\right)^{-1 / 2}.

?
Problem 2d.2

Let Y∼Nn(0,In)\mathbf{Y} \sim N_{n}\left(\mathbf{0}, \mathrm{I}_{n}\right) and let A\mathbf{A} be symmetric.

?
(a)

Show that the m.g.f. of YTAY\mathbf{Y}^{T} \mathbf{A Y} is [det⁡(In−2tA)]−1/2\left[\operatorname {det}\left(\mathbf{I}_{n}-2 t \mathbf{A}\right)\right]^{-1 / 2}.

(b)

If A\mathbf{A} is idempotent of rank⁡r\operatorname {rank} r, show that the m.g.f. is (1−2t)−r/2(1-2 t)^{-r / 2}.

(c)

Find the m.g.f. if Y∼Nn(0,Σ)\mathbf{Y} \sim N_{n}(\mathbf{0}, \mathbf{\Sigma }).

Problem 2d.3

If Y∼N2(0,I2)\mathbf{Y} \sim N_{2}\left(0, \mathbf{I}_{2}\right), find values of aa and bb such that

a(Y1−Y2)2+b(Y1+Y2)2∼χ22 a\left(Y_{1}-Y_{2}\right)^{2}+b\left(Y_{1}+Y_{2}\right)^{2} \sim \chi _{2}^{2}
?
Problem 2d.4

Suppose that Y∼N3(0,In)\mathbf{Y} \sim N_{3}\left(\mathbf{0}, \mathbf{I}_{n}\right). Show that

13[(Y1−Y2)2+(Y2−Y3)2+(Y3−Y1)2] \frac{1}{3}\left[\left(Y_{1}-Y_{2}\right)^{2}+\left(Y_{2}-Y_{3}\right)^{2}+\left(Y_{3}-Y_{1}\right)^{2}\right]

has a χ22\chi_{2}^{2} distribution. Does some multiple of

(Y1−Y2)2+(Y2−Y3)2+⋯+(Yn−1−Yn)2+(Yn−Y1)2 \left(Y_{1}-Y_{2}\right)^{2}+\left(Y_{2}-Y_{3}\right)^{2}+\cdots +\left(Y_{n-1}-Y_{n}\right)^{2}+\left(Y_{n}-Y_{1}\right)^{2}

have a chi-squared distribution for general nn?

?
Problem 2d.5

Let Y∼Nn(0,In)\mathbf{Y} \sim N_{n}\left(0, \mathrm{I}_{n}\right) and let A\mathbf{A} and B\mathbf{B} be symmetric. Show that the joint m.g.f. of YTAY\mathbf{Y}^{T} \mathbf{A Y} and YTBY\mathbf{Y}^{T} \mathbf{B Y} is [det⁡(In−2sA−2tB)]−1/2\left[\operatorname {det}\left(\mathbf{I}_{n}-2 s \mathbf{A}-2 t B\right)\right]^{-1 / 2}. Hence show that the two quadratic forms are independent if AB=0\mathbf{A B} = 0.

?