一、Wishart分布
1.1 基本概念
Wishart分布是χ2\chi^2χ2分布的多元情况,假设有随机变量ξ∼N(0,Σ)∈Rp\xi\sim N(0,\Sigma)\in\bm R^{p}ξ∼N(0,Σ)∈Rp,有iidiidiid样本X1,⋯ ,XnX_1,\cdots,X_nX1,⋯,Xn,定义X=(X1,⋯ ,Xn)T∈Rn×p\bm X=(X_1,\cdots,X_n)^T\in\bm R^{n\times p}X=(X1,⋯,Xn)T∈Rn×p,那么M=XTX=∑i=1nXiXiT∼Wp(Σ,n)(1)M=\bm X^T\bm X=\sum\limits_{i=1}^nX_iX_i^T\sim W_p(\Sigma, n)\tag1M=XTX=i=1∑nXiXiT∼Wp(Σ,n)(1)
nnn为自由度。
1.2 性质
-
如果M∼Wp(Σ,n)M\sim W_p(\Sigma, n)M∼Wp(Σ,n),而且B∈Rp×qB\in\bm R^{p\times q}B∈Rp×q,那么BTMB∼Wq(BTΣB,n)(2)B^TMB\sim W_q(B^T\Sigma B, n)\tag2BTMB∼Wq(BTΣB,n)(2)
当B=Σ−12B=\Sigma^{-\frac{1}{2}}B=Σ−21时,Σ−12MΣ−12∼Wp(Ip,n)\Sigma^{-\frac{1}{2}}M\Sigma^{-\frac{1}{2}}\sim W_p(I_p, n)Σ−21MΣ−21∼Wp(Ip,n)
-
如果M∼Wp(Σ,n)M\sim W_p(\Sigma, n)M∼Wp(Σ,n),那么E(M)=nΣE(M)=n\SigmaE(M)=nΣ
-
如果Mi∼Wp(Σ,ni),i=1,⋯ ,kM_i\sim W_p(\Sigma, n_i),\quad i=1,\cdots,kMi∼Wp(Σ,ni),i=1,⋯,k,且相互独立,那么M=∑i=1kMi∼Wp(Σ,n)(3)M=\sum_{i=1}^kM_i\sim W_p(\Sigma, n)\tag3M=i=1∑kMi∼Wp(Σ,n)(3)
其中n=∑i=1knin=\sum\limits_{i=1}^kn_in=i=1∑kni -
如果X∈Rn×p\bm X\in\bm R^{n\times p}X∈Rn×p,且Xi∼N(0,Σ)X_i\sim N(0,\Sigma)Xi∼N(0,Σ),对于对称矩阵CCC,当且仅当C2=CC^2=CC2=C时,XTCX∼W(Σ,r)(4)\bm X^TC\bm X\sim W(\Sigma, r)\tag4XTCX∼W(Σ,r)

最低0.47元/天 解锁文章
3629

被折叠的 条评论
为什么被折叠?



