Tut4 5
Tut4 5
Review
Multivariate Statistics
• Random vector
○ The elements of a vector are random variables
= ⋮
• Mean vector
○ The elements of the vector are the expected value of the elements of the random
vector
( )
= ( )= ⋮
( )
• Covariance matrix
○ X an m 1 random vector and Y an n 1 random vector
○ Variance-covariance matrix / covariance matrix is defined as
Cov X, Y E X μ X Y μ Y T
○ Properties
Cov(X,Y) = (ij) where ij = Cov(xi,yj)
Dimensions of Cov(X,Y) = m n
Cov X, Y E XY T E X E Y T
Σ Var Y E Y μ Y Y μ Y T E YY T E Y E Y T
○ Correlation
, = =
○ Correlation matrix
=
⬚ 0
= diag( ) = ⬚ ⬚
0 ⬚
• Linear combinations
○ Let A and B are constant matrices, b and c are constant vectors and X and Y are
two n 1 random vectors with Var(Y) = , then
E(AY + B X + b) = AE(Y) + BE(X) + b
Cov(AY + b,BY + c)=AB’
1
Var(AY + b)= AA’
2 Σ
p/2 1/ 2
2
Then
E(Y) = , Var(Y) = = (ij)
E(yi) = i, Var(yi) = ii, Cov(yi yj) = ij
yi ~ N (i,ii)
○ For a sample of size n, y1, …, yn, the likelihood function and log-likelihood
function are
1 1
Lμ, Σ | y1 ,, y n i 1 y i i 1 exp y i μ Σ 1 y i μ
n n T
2 Σ
p/2 1/ 2
2
n n 1 n
log L log 2 log Σ y i μ Σ 1 y i μ
p T
2 2 2 i 1
○ MLE
= = ∑ = ~ ,
= ∑ ( − )( − ) = ( − ) ( − )
−2 ~
(asymptotically)
v = number of free parameters under H1
v0 = number of free parameters under H0
2
Exercises
• Consider variables YrHgt (inches), FtFrBody (pounds), PrctFFB (percentages),
BkFat (inches), SaleHt (inches) and SalseWt (pounds) in "TUT402.csv".
a) Calculate the sample means of the variables.
b) Calculate the sample covariance matrix.
c) Calculate the sample correlation matrix.
d) Assuming multivariate normality of the 6 variables, find MLE for
([ , , ]|[ , , ] = [45,1000,80])
e) Calculate the sample means and covariance matrix of the variables when they are
measured in (feet, pounds, percentages, feet, inches, pounds)
f) Calculate the maximum log-likelihood under multivariate normality.
g) Calculate the maximum log-likelihood under multivariate normality and the mean
= [55,1200,68,0.20,50,1200].
h) Conduct a likelihood ratio test for : = [55,1200,68,0.20,50,1200] at
5%.