第二章多元正态分布
i1
i
S
n
(
i1
Xi
X)(Xi
X)
则有
S
n
X
j Xj
nXX
i1
1、
~
N
p
(,
1 n
)
2、和S相互独立
3、S ~ Wp (n 1,)
证明:
*
设
*
1
n
*
* 1 n
*
*
ij
nn
1
那么在多元正态的情形下,是否有相同的问题呢?回答
时肯定的。
定义: 设 ~ Wp (n,)和u ~ N p (μ,)相互独立,则
2 nu1u ~ T 2 ( p, n,μ)
称T2服从参数为P和n的非中心霍特林(Hotelling)分布,当。
当 μ 0时,2 nu1u 服从自由度为n的中心霍特林分布, 记为 2 nu1u ~ T 2 ( p, n) 。
a
)
,
x(a 2
)
,,
x(a na
)
1 k na x x
n a1 i1
(a) i
x 1 x a
na
i 1
(a) i
na
W
(x k na
a1 i1
(a) i
x)(x(ia)
x)
E
(x k na
a1 i1
(a) i
xa
)(x(ia)
xa )
B
分布近似。
2、Λ统计量和Λ分布
设k个总体G1,,Gk ,它们服从 N p ( (i) ,) 。分别抽出
如下的样本:
x(1) 1
,
x(1) 2
,,
x(1) n1
x1(
2
)
,
x(2 2
)
,,
x(2 n2
)
x1(
k
)
,
x(k 2
)
,,
x(k nk
)
n
n k
a1
a
xa
x1(
Z
2 p
~
(2 p)
相互独立的标准正态分布的平方和为自由度为 p 的卡方分布。
在一元正态的情形下,我们有样本的统计量
Z x ~ N (0,1) n
当总体的方差未知时,我们必须用样本的方差
S *2
n
1
1
n
(
i 1
xi
x)2
来代替总体的方差,则
t
x S*
n
~
t(n
1)
x2 (x21, x22,, x2 p ) x1 (x11, x12,, x1p )
… xn (xn1, xn2 ,, xnp )
令
1 n
n
i
i 1
样本均值
S
(X n
i 1
i
X)(Xi
X)样本叉积矩阵
则 2 (n 1)n(x )S 1(x ) ~ T 2 ( p, n 1)
定理1:若 ~ Wp (n,) ,且 0 ,n p ,则 的分布密度
为
|
a
1(n p1)
|2
exp( 1 tr1A)
F (a)
np
22
p ( p1 ) 2
n
| |2
2
,a
p
(
n
i
1)
0
i1
2
特别,当 p 1和 1 时, 服从 2 分布。
且 T 2 ( p, n) np F ( p, n p 1) n p1
定理:设 x1,x2,,xn1 是来自多元正态总体 N p (,) 的简单
随机样本,
x2 (x21, x22,, x2 p )
x1 (x11, x12,, x1p )
…
xn1 (xn11, xn1 2 ,, xn1p )
x2 )
~
n1 n2 n1n2
T
2(
p, n
1)
Sp
(n1 1)S1 (n2 1)S2 n1 n2 2
四、基于维斯特(Wishart)分布的统计量
在一元方差分析中,常常遇到基于独立的 2 分布随机变量比值
的F统计量。在多元统计分析中,起到相同作用的是统计量和 分布。
设随机矩阵 X
x21
x22
x2
p
xn1
xn2
xnp
矩阵中的每一个元素均为随机变量,则矩阵X的分布是其列 向量拉长,组成一个长向量
x x11 x1p x21 x2 p xn1 xnp 的分布。
定义 维希特(Wishart)分布的统计量 设 n 个随机向量 Xi ( X i1, X i2 ,, X ip )(i 1,2,3,, n)
2
为一正交矩阵n 令 (1 2 n ) X1 X2 Xn
由于Xi (i 1,2,3,4,n)独立同正态分布,且为正交矩阵,所以
(1 2 n )独立同正态分布
n
1 n
n
i 1
i
E(n)
1 n
x12
x22
xn 2
x21
x22
x2
p
x1
p
x2 p
xnp
xn1
xn 2
xnp
n
X il X lj
l 1
服从自由度为 n 的非中心维斯特分布,记为 ~ Wp (n,,。μ)
特别当 X是 p 阶对称阵,则X 的分布为的下三角部分组 成的长向量
p
(
,
1 n
)
n
S
(X
i1
j
X)(X
j
X)
S
n
X
i1
j Xj
nXX
S
n
i1
X
j
Xj
n n
S nn
S
n1
j 1
j
j
与S相互独立
S
n1
j 1
j
j
~ Wp (n 1,)
当 1,p 1时,由卡方分布的定义可知
为Wm (n,CC) 分布。
三、 抽样分布
定理1:设X1,X2,……Xn是来自多元正态总体Np(,) 的简单随机样本,有
x1 (x11, x12 ,, x1p )
x2 (x21, x22 ,, x2 p )
xn (xn1, xn2 ,, xnp )
令
1 n
n
第2章 抽样分布
Sampling Distributions
§1 样本的联合概率密度函数
设x ~ N p (,), 0, 则总体的密度函数为
f
( x1 ,
x2 ,,
xp
)
(2
)
p
2
1
2
exp[
1 2
(x
) 1 ( x
)]
X1,X2,……,Xn是从总体中抽取的一个简单随机样本,满足X1,
x x11 x1p x22 x2 p xp1, p1 xp, p1 xp, p
在一元正态随机变量中,我们曾经讨论了 2 分布,在多元
正态随机变量也有类似的样本分布。维希特分布(Wishart)相当
于一元统计中的 2 分布。
维希特( Wishart)分布的密度函数
X11 X12 X1p
X
X 21
X 22
X
2
p
X (1) X (2)
X
n1
X n2
X
np
n
p
X
(n)
独立同分布于
Np
(μ, )
,则随机矩阵
n
i
i
i1
A X X
x11 x21 xn1 x11 x12 x1p
证:
由于样本均值
x
~
Np
(
,
1 n
)
令
1
n 2 (X )
1
E() E[ n 2 (X )]
1
D() D[ n 2 (X )] p
1
n 2 (X ) ~ N p (o,I)
所以
2
Z12
Z
2 2
n
i 1
(2 ) p
2
1
2
exp[
1 2
(
xi
)1(xi )]