स्केल्ड व्युत्क्रम ची-वर्ग वितरण

Scaled inverse chi-squared

Probability density function
Cumulative distribution function
Parameters	$\nu >0\,$ $\tau ^{2}>0\,$
Support	$x\in (0,\infty )$
PDF	${\frac {(\tau ^{2}\nu /2)^{\nu /2}}{\Gamma (\nu /2)}}~{\frac {\exp \left[{\frac {-\nu \tau ^{2}}{2x}}\right]}{x^{1+\nu /2}}}$
CDF	$\Gamma \left({\frac {\nu }{2}},{\frac {\tau ^{2}\nu }{2x}}\right)\left/\Gamma \left({\frac {\nu }{2}}\right)\right.$
Mean	${\frac {\nu \tau ^{2}}{\nu -2}}$ for $\nu >2\,$
Mode	${\frac {\nu \tau ^{2}}{\nu +2}}$
Variance	${\frac {2\nu ^{2}\tau ^{4}}{(\nu -2)^{2}(\nu -4)}}$for $\nu >4\,$
Skewness	${\frac {4}{\nu -6}}{\sqrt {2(\nu -4)}}$for $\nu >6\,$
Ex. kurtosis	${\frac {12(5\nu -22)}{(\nu -6)(\nu -8)}}$for $\nu >8\,$
Entropy	${\frac {\nu }{2}}\!+\!\ln \left({\frac {\tau ^{2}\nu }{2}}\Gamma \left({\frac {\nu }{2}}\right)\right)$ $\!-\!\left(1\!+\!{\frac {\nu }{2}}\right)\psi \left({\frac {\nu }{2}}\right)$
MGF	${\frac {2}{\Gamma ({\frac {\nu }{2}})}}\left({\frac {-\tau ^{2}\nu t}{2}}\right)^{\!\!{\frac {\nu }{4}}}\!\!K_{\frac {\nu }{2}}\left({\sqrt {-2\tau ^{2}\nu t}}\right)$
CF	${\frac {2}{\Gamma ({\frac {\nu }{2}})}}\left({\frac {-i\tau ^{2}\nu t}{2}}\right)^{\!\!{\frac {\nu }{4}}}\!\!K_{\frac {\nu }{2}}\left({\sqrt {-2i\tau ^{2}\nu t}}\right)$

स्केल्ड व्युत्क्रम ची-वर्ग वितरण x = 1/s² के लिए वितरण है, जहां s², v स्वतंत्र सामान्य वितरण यादृच्छिक चर के वर्गों का प्रतिरूप माध्य है जिसका माध्य 0 एवं व्युत्क्रम विचरण 1/σ² = τ² है। इसलिए वितरण दो मात्राओं ν एवं τ² द्वारा परिचालित है, जिसे क्रमशः स्वतंत्रता की ची-वर्ग डिग्री की संख्या एवं स्केलिंग पैरामीटर के रूप में जाना जाता है।

स्केल किए गए व्युत्क्रम ची-वर्ग वितरण का यह सदस्य दो अन्य वितरण सदस्यों व्युत्क्रम-ची-वर्ग वितरण एवं व्युत्क्रम-गामा वितरण के निकटता से संबंधित है। व्युत्क्रम ची-वर्ग वितरण की अपेक्षा में, स्केल किए गए वितरण में अतिरिक्त पैरामीटर τ² होता है, जो वितरण को क्षैतिज एवं लंबवत रूप से मापता है, जो मूल अंतर्निहित प्रक्रिया के व्युत्क्रम-विचरण का प्रतिनिधित्व करता है। साथ ही, स्केल किए गए व्युत्क्रम ची-वर्ग वितरण को उनके योग के व्युत्क्रम के अतिरिक्त ν वर्ग विचलन के माध्य के व्युत्क्रम के वितरण के रूप में प्रस्तुत किया जाता है। इस प्रकार दोनों वितरणों में यह संबंध है कि यदि;

$X\sim {\mbox{Scale-inv-}}\chi ^{2}(\nu ,\tau ^{2})$ तब ${\frac {X}{\tau ^{2}\nu }}\sim {\mbox{inv-}}\chi ^{2}(\nu )$ होता है।

व्युत्क्रम गामा वितरण की अपेक्षा में, स्केल्ड व्युत्क्रम ची-वर्ग वितरण समान डेटा वितरण का वर्णन करता है, परन्तु भिन्न सांख्यिकीय पैरामीटर का उपयोग करता है, जो कुछ परिस्थितियों में अधिक सुविधाजनक हो सकता है। विशेष रूप से, यदि;

$X\sim {\mbox{Scale-inv-}}\chi ^{2}(\nu ,\tau ^{2})$ तब $X\sim {\textrm {Inv-Gamma}}\left({\frac {\nu }{2}},{\frac {\nu \tau ^{2}}{2}}\right)$ होता है।

किसी भी रूप का उपयोग निश्चित प्रथम व्युत्क्रम क्षण (गणित) के लिए अधिकतम एन्ट्रापी संभाव्यता वितरण, वितरण का प्रतिनिधित्व करने के लिए किया जा सकता है। $(E(1/X))$ एवं प्रथम लघुगणक क्षण $(E(\ln(X))$ है।

स्केल किए गए व्युत्क्रम ची-वर्ग वितरण का बायेसियन सांख्यिकी में भी विशेष उपयोग होता है, जो x = 1/s² के लिए पूर्वानुमानित वितरण के रूप में इसके उपयोग से कुछ सीमा तक असंबंधित है। विशेष रूप से, स्केल किए गए व्युत्क्रम ची-वर्ग वितरण का उपयोग सामान्य वितरण के विचरण पैरामीटर के लिए संयुग्म पूर्व के रूप में किया जा सकता है। इस संदर्भ में स्केलिंग पैरामीटर को τ² के अतिरिक्त σ₀² द्वारा प्रदर्शित किया गया है, एवं इसकी भिन्न व्याख्या है। इसके अतिरिक्त प्रोग्राम को सामान्यतः व्युत्क्रम-गामा वितरण सूत्रीकरण का उपयोग करके प्रस्तुत किया गया है; चूँकि, कुछ लेखक, विशेष रूप से गेलमैन एट अल (1995/2004) का अनुसरण कर रहे हैं जिसका तर्क है कि व्युत्क्रम ची-स्क्वेर्ड पैरामीट्रिज़ेशन अधिक सहज है।

विशेषता

स्केल किए गए व्युत्क्रम ची-वर्ग वितरण की संभाव्यता घनत्व फलन डोमेन $x>0$ पर विस्तृत है एवं

$f(x;\nu ,\tau ^{2})={\frac {(\tau ^{2}\nu /2)^{\nu /2}}{\Gamma (\nu /2)}}~{\frac {\exp \left[{\frac {-\nu \tau ^{2}}{2x}}\right]}{x^{1+\nu /2}}}$ है,

जहाँ $\nu$ स्वतंत्रता की डिग्री (सांख्यिकी) पैरामीटर है एवं $\tau ^{2}$ स्केल पैरामीटर है, संचयी वितरण फलन

$F(x;\nu ,\tau ^{2})=\Gamma \left({\frac {\nu }{2}},{\frac {\tau ^{2}\nu }{2x}}\right)\left/\Gamma \left({\frac {\nu }{2}}\right)\right.$

$=Q\left({\frac {\nu }{2}},{\frac {\tau ^{2}\nu }{2x}}\right)$ है,

जहाँ $\Gamma (a,x)$ अधूरा गामा फलन है, $\Gamma (x)$ गामा फलन है एवं $Q(a,x)$ नियमित गामा फलन है। विशेषता फलन है:

$\varphi (t;\nu ,\tau ^{2})=$

${\frac {2}{\Gamma ({\frac {\nu }{2}})}}\left({\frac {-i\tau ^{2}\nu t}{2}}\right)^{\!\!{\frac {\nu }{4}}}\!\!K_{\frac {\nu }{2}}\left({\sqrt {-2i\tau ^{2}\nu t}}\right),$

जहाँ $K_{\frac {\nu }{2}}(z)$ दूसरे प्रकार का संशोधित बेसेल फलन है।

पैरामीटर अनुमान

$\tau ^{2}$ की अधिकतम संभावना अनुमान

$\tau ^{2}=n/\sum _{i=1}^{n}{\frac {1}{x_{i}}}$ है,

${\frac {\nu }{2}}$ की अधिकतम संभावना अनुमान न्यूटन की विधि का उपयोग करके पाया जा सकता है:

$\ln \left({\frac {\nu }{2}}\right)-\psi \left({\frac {\nu }{2}}\right)={\frac {1}{n}}\sum _{i=1}^{n}\ln \left(x_{i}\right)-\ln \left(\tau ^{2}\right),$

जहाँ $\psi (x)$ डिगामा फलन है। माध्य का सूत्र लेकर एवं इसका निवारण करके प्रारंभिक अनुमान $\nu$ प्राप्त किया जा सकता है। ${\bar {x}}={\frac {1}{n}}\sum _{i=1}^{n}x_{i}$ प्रतिरूप माध्य हो, फिर प्रारंभिक अनुमान $\nu$ द्वारा दिया गया है:

${\frac {\nu }{2}}={\frac {\bar {x}}{{\bar {x}}-\tau ^{2}}}$ है।

सामान्य वितरण के विचरण का बायेसियन अनुमान

सामान्य वितरण के विचरण के बायेसियन अनुमान में स्केल्ड व्युत्क्रम ची-वर्ग वितरण का दूसरा महत्वपूर्ण अनुप्रयोग है।

बेयस प्रमेय के अनुसार, ब्याज की मात्राओं के लिए पश्च संभाव्यता वितरण, मात्राओं एवं संभावना फलन के लिए पूर्व वितरण के उत्पाद के समानुपाती होता है:

$p(\sigma ^{2}|D,I)\propto p(\sigma ^{2}|I)\;p(D|\sigma ^{2})$ जहां D डेटा का प्रतिनिधित्व करता है एवं I, σ² के विषय में किसी प्रारंभिक जानकारी का प्रतिनिधित्व करता है, जो हमारे पास पूर्व से ही हो सकता है।

सबसे सरल परिदृश्य तब उत्पन्न होता है जब माध्य μ पूर्व से ही ज्ञात हो; या, वैकल्पिक रूप से, यदि यह σ² की सशर्त संभावना है, जो कि μ के विशेष कल्पित मान के लिए लिया गया है।

तब संभाव्यता पद L(σ)²|D) = p(D|p²) का परिचित रूप

${\mathcal {L}}(\sigma ^{2}|D,\mu )={\frac {1}{\left({\sqrt {2\pi }}\sigma \right)^{n}}}\;\exp \left[-{\frac {\sum _{i}^{n}(x_{i}-\mu )^{2}}{2\sigma ^{2}}}\right]$ है।

इसे पुनर्स्केलिंग-अपरिवर्तनीय पूर्व p(σ)²|I) = 1/s² के साथ संयोजित करना, जिसके विषय में तर्क दिया जा सकता है (उदाहरण के लिए मानक विचलन पैरामीटर के साथ जेफरीज़ का अनुसरण करते हुए) कि यह समस्या में σ² के लिए पूर्व संभव सबसे कम सूचनापूर्ण है, संयुक्त पश्चवर्ती संभावना देता है:

$p(\sigma ^{2}|D,I,\mu )\propto {\frac {1}{\sigma ^{n+2}}}\;\exp \left[-{\frac {\sum _{i}^{n}(x_{i}-\mu )^{2}}{2\sigma ^{2}}}\right]$,

इस फॉर्म को स्केल किए गए व्युत्क्रम ची-वर्ग वितरण के रूप में पैरामीटर ν = n एवं के साथ τ² = s² = (1/n) Σ (x_i-μ)² के साथ पहचाना जा सकता है।

गेलमैन एट अल की टिप्पणी है कि इस वितरण की पुन: उपस्थिति, जिसे पनिवारणे प्रतिरूप संदर्भ में देखा गया था, उल्लेखनीय लग सकता है; परन्तु पनिवारणे के विकल्प को देखते हुए परिणाम आश्चर्यजनक नहीं है। विशेष रूप से, σ² के लिए पनिवारणे पुनर्स्केलिंग-अपरिवर्तनीय का विकल्प का परिणाम यह है कि σ ²/s² के अनुपात की संभावना रूप कंडीशनिंग चर से स्वतंत्र) समान होता है जब s² पर अनुकूलित किया जाता है, जैसे कि जब σ² पर अनुकूलित किया जाता है:

$p({\tfrac {\sigma ^{2}}{s^{2}}}|s^{2})=p({\tfrac {\sigma ^{2}}{s^{2}}}|\sigma ^{2})$,

प्रतिरूप-सिद्धांत विषय में, σ² पर अनुकूलित, (1/s²) के लिए संभाव्यता वितरण) स्केल्ड व्युत्क्रम ची-वर्ग वितरण है; एवं इसलिए s² पर अनुकूलित σ² के लिए संभाव्यता वितरण, स्केल-अज्ञेयवादी पूर्व दिया गया स्केल्ड व्युत्क्रम ची-वर्ग वितरण भी है।

पूर्व सूचनात्मक के रूप में उपयोग करें

यदि σ² के संभावित मूल्यों के विषय में अधिक सूचना है, स्केल्ड व्युत्क्रम ची-स्क्वायर सदस्य से वितरण, जैसे स्केल-इनव-χ²(n₀, s₀²), σ² के लिए अधिक सूचनापूर्ण पूर्व का प्रतिनिधित्व करने के लिए सुविधाजनक रूप हो सकता है, n₀ के परिणाम से पूर्व अवलोकन के परिणाम से (चूँकि n₀ आवश्यक नहीं कि पूर्ण संख्या हो):

$p(\sigma ^{2}|I^{\prime },\mu )\propto {\frac {1}{\sigma ^{n_{0}+2}}}\;\exp \left[-{\frac {n_{0}s_{0}^{2}}{2\sigma ^{2}}}\right]$,

इस प्रकार के पूर्व से पश्चवर्ती वितरण को बढ़ावा मिलता है,

$p(\sigma ^{2}|D,I^{\prime },\mu )\propto {\frac {1}{\sigma ^{n+n_{0}+2}}}\;\exp \left[-{\frac {ns^{2}+n_{0}s_{0}^{2}}{2\sigma ^{2}}}\right]$,

जो स्वयं स्केल्ड व्युत्क्रम ची-वर्ग वितरण है। इस प्रकार स्केल किए गए व्युत्क्रम ची-वर्ग वितरण σ² अनुमान के लिए सुविधाजनक संयुग्मित पूर्व सदस्य हैं।

माध्य अज्ञात होने पर विचरण का अनुमान

यदि माध्य ज्ञात नहीं है, तो इसके लिए जो सबसे असूचनात्मक पूर्व लिया जा सकता है, वह संभवतः अनुवाद-अपरिवर्तनीय पूर्व p(μ|I) ∝ स्थिरांक है, जो μ एवं σ² के लिए निम्नलिखित संयुक्त पश्च वितरण देता है,

${\begin{aligned}p(\mu ,\sigma ^{2}\mid D,I)&\propto {\frac {1}{\sigma ^{n+2}}}\exp \left[-{\frac {\sum _{i}^{n}(x_{i}-\mu )^{2}}{2\sigma ^{2}}}\right]\\&={\frac {1}{\sigma ^{n+2}}}\exp \left[-{\frac {\sum _{i}^{n}(x_{i}-{\bar {x}})^{2}}{2\sigma ^{2}}}\right]\exp \left[-{\frac {n(\mu -{\bar {x}})^{2}}{2\sigma ^{2}}}\right]\end{aligned}}$

σ² के लिए सीमांत पश्च μ वितरण पर एकीकृत करके संयुक्त पश्च वितरण से प्राप्त किया जाता है,

${\begin{aligned}p(\sigma ^{2}|D,I)\;\propto \;&{\frac {1}{\sigma ^{n+2}}}\;\exp \left[-{\frac {\sum _{i}^{n}(x_{i}-{\bar {x}})^{2}}{2\sigma ^{2}}}\right]\;\int _{-\infty }^{\infty }\exp \left[-{\frac {n(\mu -{\bar {x}})^{2}}{2\sigma ^{2}}}\right]d\mu \\=\;&{\frac {1}{\sigma ^{n+2}}}\;\exp \left[-{\frac {\sum _{i}^{n}(x_{i}-{\bar {x}})^{2}}{2\sigma ^{2}}}\right]\;{\sqrt {2\pi \sigma ^{2}/n}}\\\propto \;&(\sigma ^{2})^{-(n+1)/2}\;\exp \left[-{\frac {(n-1)s^{2}}{2\sigma ^{2}}}\right]\end{aligned}}$

यह पुनः $\scriptstyle {n-1}\;$ एवं $\scriptstyle {s^{2}=\sum (x_{i}-{\bar {x}})^{2}/(n-1)}$ मापदंडों के साथ स्केल्ड व्युत्क्रम ची-वर्ग वितरण है।

संबंधित वितरण

• यदि $X\sim {\mbox{Scale-inv-}}\chi ^{2}(\nu ,\tau ^{2})$ तब $kX\sim {\mbox{Scale-inv-}}\chi ^{2}(\nu ,k\tau ^{2})\,$होता है।
• यदि $X\sim {\mbox{inv-}}\chi ^{2}(\nu )\,$ (विपरीत-ची-वर्ग वितरण) तो $X\sim {\mbox{Scale-inv-}}\chi ^{2}(\nu ,1/\nu )\,$होता है।
• यदि $X\sim {\mbox{Scale-inv-}}\chi ^{2}(\nu ,\tau ^{2})$ तब ${\frac {X}{\tau ^{2}\nu }}\sim {\mbox{inv-}}\chi ^{2}(\nu )\,$ (व्युत्क्रम-ची-वर्ग वितरण) होता है।
• यदि $X\sim {\mbox{Scale-inv-}}\chi ^{2}(\nu ,\tau ^{2})$ तब $X\sim {\textrm {Inv-Gamma}}\left({\frac {\nu }{2}},{\frac {\nu \tau ^{2}}{2}}\right)$ (विपरीत-गामा वितरण) होता है।
• स्केल्ड व्युत्क्रम ची वर्ग वितरण प्रकार 5 पियर्सन वितरण का विशेष विषय है।

संदर्भ

• Gelman A. et al (1995), Bayesian Data Analysis, pp 474–475; also pp 47, 480