内容来源
应用多元统计分析 北京大学出版社 高惠璇编著
1
设 X ∼ N p ( μ , Σ ) , Σ > 0 X\sim N_p(\mu,\Sigma),\Sigma>0 X∼Np(μ,Σ),Σ>0
则 X ′ Σ − 1 X ∼ χ 2 ( p , δ ) X'\Sigma^{-1}X\sim\chi^2(p,\delta) X′Σ−1X∼χ2(p,δ),其中 δ = μ ′ Σ − 1 μ \delta=\mu'\Sigma^{-1}\mu δ=μ′Σ−1μ
证明
因为 Σ > 0 \Sigma>0 Σ>0,所以存在可逆矩阵 C C C
有 Σ = C C ′ \Sigma=CC' Σ=CC′
令 Y = C − 1 X Y=C^{-1}X Y=C−1X,则
Y ∼ N p ( C − 1 μ , C − 1 Σ ( C − 1 ) ′ ) = N p ( C − 1 μ , I ) Y\sim N_p(C^{-1}\mu,C^{-1}\Sigma(C^{-1})')= N_p(C^{-1}\mu,I) Y∼Np(C−1μ,C−1Σ(C−1)′)=Np(C−1μ,I)
且
X ′ Σ − 1 X = Y ′ C ′ Σ − 1 C Y = Y ′ Y ∼ χ 2 ( p , δ ) X'\Sigma^{-1}X=Y'C'\Sigma^{-1}CY=Y'Y\sim\chi^2(p,\delta) X′Σ−1X=Y′C′Σ−1CY=Y′Y∼χ2(p,δ)
其中 δ = ( C − 1 μ ) ′ ( C − 1 μ ) = μ ′ Σ − 1 μ \delta=(C^{-1}\mu)'(C^{-1}\mu)=\mu'\Sigma^{-1}\mu δ=(C−1μ)′(C−1μ)=μ′Σ−1μ
2
设 X ∼ N p ( μ , Σ ) , Σ > 0 X\sim N_p(\mu,\Sigma),\Sigma>0 X∼Np(μ,Σ),Σ>0
A A A 为对称矩阵, r a n k ( A ) = r rank(A)=r rank(A)=r
则
( X − μ ) ′ A ( X − μ ) ∼ χ 2 ( r ) ⟺ Σ A Σ A Σ = Σ A Σ (X-\mu)'A(X-\mu)\sim\chi^2(r)\Longleftrightarrow \Sigma A\Sigma A\Sigma=\Sigma A\Sigma (X−μ)′A(X−μ)∼χ2(r)⟺ΣAΣAΣ=ΣAΣ
证明
因为 Σ > 0 \Sigma>0 Σ>0
所以它的平方根矩阵 Σ 1 / 2 > 0 \Sigma^{1/2}>0 Σ1/2>0
所以存在正交矩阵 Γ \Gamma Γ 使得
Σ 1 / 2 = Γ d i a g ( λ 1 , ⋯ , λ r ) Γ ′ \Sigma^{1/2}= \Gamma diag\Big(\sqrt{\lambda_1},\cdots,\sqrt{\lambda_r}\Big)\Gamma' Σ1/2=Γdiag(λ1,⋯,λr)Γ′
其中 λ i \lambda_i λi 是 Σ \Sigma Σ 的特征值
令 Y = Σ − 1 / 2 ( X − μ ) ∼ N p ( 0 , I ) Y=\Sigma^{-1/2}(X-\mu)\sim N_p(0,I) Y=Σ−1/2(X−μ)∼Np(0,I),则
( X − μ ) ′ A ( X − μ ) = Y ′ Σ 1 / 2 A Σ 1 / 2 Y (X-\mu)'A(X-\mu)=Y'\Sigma^{1/2} A\Sigma^{1/2}Y (X−μ)′A(X−μ)=Y′Σ1/2AΣ1/2Y
设 C = Σ 1 / 2 A Σ 1 / 2 C=\Sigma^{1/2} A\Sigma^{1/2} C=Σ1/2AΣ1/2,结合上篇笔记-分量独立的n维正态随机向量的二次型的结论有
Y ′ C Y ∼ χ 2 ( p ) ⟺ C 2 = C Y'CY\sim\chi^2(p)\Longleftrightarrow C^2=C Y′CY∼χ2(p)⟺C2=C
即
Σ 1 / 2 A Σ 1 / 2 ⋅ Σ 1 / 2 A Σ 1 / 2 = Σ 1 / 2 A Σ 1 / 2 \Sigma^{1/2} A\Sigma^{1/2}\cdot\Sigma^{1/2} A\Sigma^{1/2}=\Sigma^{1/2} A\Sigma^{1/2} Σ1/2AΣ1/2⋅Σ1/2AΣ1/2=Σ1/2AΣ1/2
3
设 X ∼ N p ( μ , Σ ) , Σ > 0 设 X\sim N_p(\mu,\Sigma),\Sigma>0 设X∼Np(μ,Σ),Σ>0
A , B A,B A,B 为 p p p 阶非对称矩阵,则
( X − μ ) ′ A ( X − μ ) 与 ( X − μ ) ′ B ( X − μ ) 独立 ⟺ Σ A Σ B Σ = 0 (X-\mu)'A(X-\mu)与(X-\mu)'B(X-\mu)独立\\\Longleftrightarrow \Sigma A\Sigma B\Sigma=0 (X−μ)′A(X−μ)与(X−μ)′B(X−μ)独立⟺ΣAΣBΣ=0