मैं एक नॉर्मल-विशरट पोस्टीरियर के व्युत्पन्न पर काम कर रहा हूं, लेकिन मैं एक पैरामीटर (स्केल मैट्रिक्स के पीछे, नीचे देखें) में फंस गया हूं।
केवल संदर्भ और पूर्णता के लिए, यहाँ मॉडल और बाकी व्युत्पत्तियाँ हैं:
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Λ))
हम पोस्टीरियर नॉर्मल-विशरट ( ) चाहते हैं, जो अपघटित हो सकता है | साथ ही :μ,Λ|μ′,κ′,υ′,W′N(μ|μ,κ′Λ)W(Λ|υ′,W′)
आजादी के degressυ′
संभावना और Wishart के पहले कारकों को मिलाने से हमें बाद में Wishart कारक का पहला कारक मिलता है:
और इसलिए हमारे पास पीछे के पहले पैरामीटर हैं:
|Λ|υ0+N−D−12
υ′=υ0+N
स्केल फैक्टरκ′
हम पहचान तत्वों से घिरा हुआ और को खोजने के लिए जो पूर्व संभावना से अद्यतन किया जाता है:
और इसलिए हमें दूसरा पैरामीटर मिला:
μTμκ0Λ
μT((κ0+N)Λ)μ
κ′=κ0+N
मीनμ′
तीसरा पैरामीटर यह पहचानने से आता है कि :
अंदर क्या है
और इसलिए हमें तीसरा पैरामीटर मिला:
2μT...
2μT(ΛNx¯¯¯+κ0Λμ0)(ΛNx¯¯¯+κ0Λμ0)(Nx¯¯¯+κ0μ0)=2μTκ′Λμ′=κ′Λμ′=κ′μ′
μ′=1k′(Nx¯¯¯+κ0μ0)
स्केल मैट्रिक्सW′
और चौथा पैरामीटर शेष मापदंडों पर काम करने से आता है:
tr(W′−1Λ)=tr(W−10Λ)+∑i=1NxTiΛxi+μT0κ0Λμ0=tr(W−10Λ)+∑i=1Ntr(xTiΛxi)+tr(μT0κ0Λμ0)=tr(W−10Λ+∑i=1NxTiΛxi+μT0κ0Λμ0)
यहां से कैसे जाना है (यदि मैंने अब तक कोई गलती नहीं की है) और लिए मानक समाधान प्राप्त करें ?W′
1 संपादित करें :
अब हम मानक के रूप में दो वर्गों को प्राप्त करने के लिए कुछ कारकों को फिर से व्यवस्थित करते हैं, जोड़ते हैं और कुछ कारकों को प्रतिस्थापित करते हैं:
tr(W′−1Λ)==tr(W−1Λ+∑i=1N(xTiΛxi+x¯¯¯TΛx¯¯¯−2xTiΛx¯¯¯)+κ0(μT0Λμ0+x¯¯¯TΛx¯¯¯−2x¯¯¯TΛμ0)−∑i=1Nx¯¯¯TΛx¯¯¯+2∑i=1NxTiΛx¯¯¯−κ0x¯¯¯TΛx¯¯¯+2κ0x¯¯¯TΛμ0)tr(W−1Λ+∑i=1N(xi−x¯¯¯)Λ(xi−x¯¯¯)T+κ0(x¯¯¯−μ0)Λ(x¯¯¯−μ0)T−Nx¯¯¯Λx¯¯¯T+2Nx¯¯¯Λx¯¯¯T−κ0x¯¯¯Λx¯¯¯T+2κ0x¯¯¯ΛμT0)
हम उन कारकों को सरल करते हैं जो वर्गों से बाहर रहते हैं:
tr(W′−1Λ)=tr(W−1Λ+∑i=1N(xi−x¯¯¯)TΛ(xi−x¯¯¯)+κ0(x¯¯¯−μ0)TΛ(x¯¯¯−μ0)+(N−κ0)x¯¯¯TΛx¯¯¯+2κ0x¯¯¯TΛμ0)
संपादित 2 ( अनुवर्ती @bdeonovic के उत्तर के लिए धन्यवाद )
ट्रेस चक्रीय है, इसलिए । फिर:
और फिर:
tr(ABC)=tr(BCA)=tr(CAB)
tr(W′−1Λ)=tr(W−1Λ+∑i=1N(xi−x¯¯¯)(xi−x¯¯¯)TΛ+κ0(x¯¯¯−μ0)(x¯¯¯−μ0)TΛ+(N−κ0)x¯¯¯x¯¯¯TΛ+2κ0x¯¯¯μT0Λ)
tr(W′−1)=tr(W−1+∑i=1N(xi−x¯¯¯)(xi−x¯¯¯)T+κ0(x¯¯¯−μ0)(x¯¯¯−μ0)T+(N−κ0)x¯¯¯x¯¯¯T+2κ0x¯¯¯μT0)
लगभग! लेकिन अभी भी वहाँ नहीं है। लक्ष्य है:
W−1+∑i=1N(xi−x¯¯¯)(xi−x¯¯¯)T+κ0Nκ0+N(x¯¯¯−μ0)(x¯¯¯−μ0)T