나는 Normal-Wishart 후부의 파생을 연구하고 있지만 매개 변수 중 하나 (스케일 매트릭스의 후부, 아래쪽 참조)에 붙어 있습니다.
맥락과 완전성을 위해 다음은 모델과 나머지 파생물입니다.
xiμΛ∼N(μ,Λ)∼N(μ0,(κ0Λ)−1)∼W(υ0,W0)
세 가지 요소 각각의 확장 된 형태는 다음과 같습니다 (최대 비례 상수).
가능성 :
N(xi|μ,Λ)∝|Λ|N/2exp(−12∑i=1N(xTiΛxi−2μTΛxi+μTΛμ))
일반적인 이전 :
N(μ|(μ0,κ0Λ)−1)∝|Λ|1/2exp(−12(μTκ0Λμ−2μTκ0Λμ0+μT0κ0Λμ0))
Wishart 이전 :
W(Λ|υ0,W0)∝|Λ|υ0−D−12exp(−12tr(W−10Λ))
우리는 다음 과 같이 분해 될 수있는 사후 Normal-Wishart ( ) 뿐만 아니라 :N ( μ | μ , κ " Λ ) W ( Λ | υ ' , W ' )μ,Λ|μ′,κ′,υ′,W′N(μ|μ,κ′Λ)W(Λ|υ′,W′)
자유의υ′
가능성의 첫 번째 요소와 Wishart를 병합함으로써 우리는 뒤에서 Wishart 요소의 첫 번째 요소를 얻습니다.
따라서 후부의 첫 번째 매개 변수가 있습니다 :
υ'=υ0+N
|Λ|υ0+N−D−12
υ′=υ0+N
스케일 팩터κ′
및 로 둘러 싸인 요소를 식별하여 이전 이 다음 가능성에 의해 업데이트되는 사람을 찾습니다 .
따라서
μ κ 0 Λ μ TμTμκ0Λ κ ' = κ 0 +N
μT((κ0+N)Λ)μ
κ′=κ0+N
평균μ′
세 번째 매개 변수는 내부의 내용을 식별하는 데 있습니다 :
따라서 세 번째 매개 변수가 있습니다 :
2 μ T ( Λ N ¯ x + κ 0 Λ μ 0 )2μT...μ ' = 1
2μT(ΛNx¯¯¯+κ0Λμ0)(ΛNx¯¯¯+κ0Λμ0)(Nx¯¯¯+κ0μ0)=2μTκ′Λμ′=κ′Λμ′=κ′μ′
μ′=1k′(Nx¯¯¯+κ0μ0)
스케일 매트릭스W′
그리고 네 번째 매개 변수는 나머지 매개 변수에 대한 작업에서 비롯됩니다.
t r ( 승' - 1Λ )= t r ( W− 10Λ ) + ∑나는 = 1엔엑스티나는Λ x나는+ μ티0κ0Λ μ0= t r ( W− 10Λ ) + ∑나는 = 1엔t r ( x티나는Λ x나는) + t r ( μ티0κ0Λ μ0)= t r ( W− 10Λ + ∑나는 = 1엔엑스티나는Λ x나는+ μ티0κ0Λ μ0)
여기에서 계속 진행하고 (지금까지 실수하지 않은 경우) 대한 표준 솔루션을 얻는 방법은 무엇입니까?여'
편집 1 :
이제 우리는 표준 솔루션에서와 같이 두 개의 제곱을 얻기 위해 항을 재정렬하고 몇 가지 요소를 더하고 뺍니다.
티 r ( 승' - 1Λ ) ==t r ( 승− 1Λ+ ∑나는 = 1엔( x티나는Λ x나는+ x¯¯¯티Λ x¯¯¯− 2 x티나는Λ x¯¯¯)+ κ0( μ티0Λ μ0+ x¯¯¯티Λ x¯¯¯− 2 x¯¯¯티Λ μ0)− ∑나는 = 1엔엑스¯¯¯티Λ x¯¯¯+ 2 ∑나는 = 1엔엑스티나는Λ x¯¯¯− κ0엑스¯¯¯티Λ x¯¯¯+ 2 κ0엑스¯¯¯티Λ μ0)t r ( 승− 1Λ + ∑나는 = 1엔( x나는− x¯¯¯) Λ ( x나는− x¯¯¯)티+ κ0( x¯¯¯−μ0) Λ (x¯¯¯−μ0)티− N엑스¯¯¯Λ x¯¯¯티+ 2 N엑스¯¯¯Λ x¯¯¯티−κ0엑스¯¯¯Λ x¯¯¯티+ 2κ0엑스¯¯¯Λμ티0)
우리는 사각형에 남아있는 요소들을 단순화합니다 :
t r ( 승' - 1Λ ) =t r ( 승− 1Λ + ∑나는 = 1엔( x나는− x¯¯¯)티Λ ( x나는− x¯¯¯) + κ0( x¯¯¯−μ0)티Λ ( x¯¯¯−μ0)+ ( N−κ0) x¯¯¯티Λ x¯¯¯+ 2κ0엑스¯¯¯티Λμ0)
편집 2 ( @bdeonovic의 답변 덕분에 후속 조치 )
추적은 주기적이므로 입니다. 그런 다음
다음 :
t r ( AB C) = t r ( B CA ) = t r ( CB )
t r ( 승' - 1Λ ) =t r ( 승− 1Λ + ∑나는 = 1엔( x나는− x¯¯¯) ( x나는− x¯¯¯)티Λ + κ0( x¯¯¯−μ0) ( x¯¯¯−μ0)티Λ+ ( N−κ0) x¯¯¯엑스¯¯¯티Λ + 2 κ0엑스¯¯¯μ티0Λ )
t r ( 승' - 1) =t r ( 승− 1+ ∑나는 = 1엔( x나는− x¯¯¯) ( x나는− x¯¯¯)티+κ0( x¯¯¯−μ0) ( x¯¯¯−μ0)티+ ( N−κ0) x¯¯¯엑스¯¯¯티+ 2κ0엑스¯¯¯μ티0)
거의! 그러나 여전히 거기에 없습니다. 목표는 다음과 같습니다.
여− 1+ ∑나는 = 1엔( x나는− x¯¯¯) ( x나는− x¯¯¯)티+κ0엔κ0+ N( x¯¯¯−μ0) ( x¯¯¯−μ0)티