Wishart distribution

deejayosamu·2026년 1월 29일

다변량통계

목록 보기
10/11

: Multivariate extension of χ2\chi^2 distribution

  • pdf

    f(w)=det(w)(rp1)/22rp/2Γp(r/2)det(Σ)r/2exp(12tr(Σ1w))f(w)=\frac{det(w)^{(r-p-1)/2}}{2^{rp/2} \Gamma_p(r/2) det(\Sigma)^{r/2}} exp(-\frac{1}{2} tr(\Sigma^{-1}w))
    for wRp×p:positive definite,r>p1w \in \mathbb{R}^{p \times p}: \text{positive definite}, r>p-1 and ΣRp×p:positive definite\Sigma \in \mathbb{R}^{p \times p}: \text{positive definite}

  • 동치명제

    WWp(r,Σ),WRp×pW \sim W_p(r,\Sigma),W \in \mathbb{R}^{p \times p}
    <=> Wdj=1rZjZjT where ZjiidNp(0,Σ)(j=1,,r)W \overset{d}{\equiv} \sum_{j=1}^{r} Z_j Z_j^T \text{ where } Z_j \overset{iid}{\sim} N_p(0,\Sigma)(j=1,\cdots,r)
    <=> For any dRp,dTWddTΣdχr2d \in \mathbb{R}^p, d^T W d \sim d^T \Sigma d \chi^2_r

  • Properties

    ① Suppose X1,,XniidNp(0,Σ),X=(X1,,Xn)TRn×p and ARn×n:symmtricX_1,\cdots,X_n \overset{iid}{\sim} N_p(0,\Sigma), X=(X_1,\cdots,X_n)^T \in \mathbb{R}^{n \times p} \text{ and } A \in \mathbb{R}^{n \times n}: \text{symmtric}
    If A2=AA^2=A and r=tr(A)r=tr(A), we have XTAXWp(r,Σ)X^T A X \sim W_p(r,\Sigma)

    ② If W1Wp(r1,Σ),W2Wp(r2,Σ) and W1W2W_1 \sim W_p(r_1,\Sigma),W_2 \sim W_p(r_2,\Sigma) \text{ and } W_1 \perp W_2, then W1+W2Wp(r1+r2,Σ)W_1 + W_2 \sim W_p(r_1+r_2, \Sigma)

    ③ If WWp(r,Σ)W \sim W_p(r, \Sigma), then CWCTWq(r,CΣCT)CWC^T \sim W_q(r, C \Sigma C^T) for any CRq×pC \in \mathbb{R}^{q \times p} with rank(C)=qprank(C)=q \leq p

    ④ If W=(wij)Wp(r,Σ) where Σ=(σij)W=(w_{ij}) \sim W_p(r,\Sigma) \text{ where } \Sigma=(\sigma_{ij}), then E(w)=rΣE(w)=r \Sigma and Var(wij)=r(σij2+σiiσjj)Var(w_{ij})=r(\sigma^2_{ij}+\sigma_{ii}\sigma_{jj})
    pf)

    pf_1

    pf_2

    pf_3

    pf_4

0개의 댓글