logo

랜덤벡터 이차형식의 기대값 📂수리통계학

랜덤벡터 이차형식의 기대값

공식

모평균 벡터 μRn\mu \in \mathbb{R}^{n}공분산행렬 ΣRn×n\Sigma \in \mathbb{R}^{n \times n} 에 대해 랜덤벡터 X\mathbf{X}X(μ,Σ)\mathbf{X} \sim \left( \mu , \Sigma \right) 라 하자. 대칭행렬 AA 에 대해 랜덤벡터 이차형식 Q=XTAXQ = \mathbf{X}^{T} A \mathbf{X}기대값은 다음과 같다. E(Q)=trAΣ+μTAμ E (Q) = \operatorname{tr} A \Sigma + \mu^{T} A \mu 여기서 μT\mu^{T}μ\mu전치행렬이고, tr\operatorname{tr}트레이스다.

따름정리

XN(μ,Σ)\mathbf{X} \sim N(\mu, \Sigma)이면,

E((Xμ)TA(Xμ))=tr(AΣ)E((Xμ1)TA(Xμ1))=tr(AΣ)+(μμ1)TA(μμ1) \begin{align*} E((\mathbf{X} - \mu)^{\mathsf{T}}A (\mathbf{X} - \mu)) &= \tr (A \Sigma) \\ E((\mathbf{X} - \mu_{1})^{\mathsf{T}}A (\mathbf{X} - \mu_{1})) &= \tr (A \Sigma) + (\mu - \mu_{1})^{\mathsf{T}} A (\mu - \mu_{1}) \end{align*}

증명

Z=Xμ\mathbf{Z} = \mathbf{X} - \mu 라고 하면 ZN(0,Σ)\mathbf{Z} \sim N(\mathbf{0}, \Sigma) 이므로,

E(ZTAZ)=tr(AΣ) E(\mathbf{Z}^{\mathsf{T}} A \mathbf{Z}) = \tr(A\Sigma)

Z=Xμ1\mathbf{Z} = \mathbf{X} - \mu_{1} 라고 하면, ZN(μμ1,Σ)\mathbf{Z} \sim N(\mu - \mu_{1}, \Sigma) 이므로,

E(ZTAZ)=tr(AΣ)+(μμ1)TA(μμ1) E(\mathbf{Z}^{\mathsf{T}} A \mathbf{Z}) = \tr(A\Sigma) + (\mu - \mu_{1})^{\mathsf{T}} A (\mu - \mu_{1})

유도 1

대각합의 순환 성질: tr(ABC)=tr(BCA)=tr(CAB) \operatorname{tr} (ABC) = \operatorname{tr} (BCA) = \operatorname{tr} (CAB)

랜덤벡터의 기대값과 트레이스: E(tr(X))=tr(E(X))E(\tr(\mathbf{X})) = \tr(E(\mathbf{X}))

공분산행렬: μRp\mathbf{\mu} \in \mathbb{R}^{p}μ:=(EX1,,EXp)\mathbf{\mu} := \left( EX_{1} , \cdots , EX_{p} \right) 와 같이 주어져 있다고 하면 Cov(X)=E[XXT]μμT \operatorname{Cov} (\mathbf{X}) = E \left[ \mathbf{X} \mathbf{X}^{T} \right] - \mathbf{\mu} \mathbf{\mu}^{T}

E(Q)=E(trXTAX)=E(trAXXT)=trAE(XXT)=trA(Σ+μμT)=trAΣ+trAμμT=trAΣ+trμTAμ=trAΣ+μTAμ \begin{align*} & E (Q) \\ =& E \left( \operatorname{tr} \mathbf{X}^{T} A \mathbf{X} \right) \\ =& E \left( \operatorname{tr} A \mathbf{X} \mathbf{X}^{T} \right) \\ =& \operatorname{tr} A E \left( \mathbf{X} \mathbf{X}^{T} \right) \\ =& \operatorname{tr} A \left( \Sigma + \mu \mu^{T} \right) \\ =& \operatorname{tr} A \Sigma + \operatorname{tr} A \mu \mu^{T} \\ =& \operatorname{tr} A \Sigma + \operatorname{tr} \mu^{T} A \mu \\ =& \operatorname{tr} A \Sigma + \mu^{T} A \mu \end{align*}


  1. Hogg et al. (2018). Introduction to Mathematical Statistcs(8th Edition): p556. ↩︎