توزيع كاي تربيع العكسي المُقاس

توزيع كاي تربيع العكسي المُقاسψinv-χ2(ν){\displaystyle \psi \,{\mbox{inv-}}\chi ^{2}(\nu )}، أينψ{\displaystyle \psi }يمثل معامل المقياس، ويساوي توزيع ويشارت العكسي أحادي المتغيردبليو-1(ψ،ν){\displaystyle {\mathcal {W}}^{-1}(\psi ,\nu )}بدرجات من الحريةν{\displaystyle \nu }.

ترتبط هذه العائلة من توزيعات كاي تربيع العكسية المُقاسة بتوزيع كاي تربيع العكسي وبتوزيع كاي تربيع :

لوXψinv-χ2(ν){\displaystyle X\sim \psi \,{\mbox{inv-}}\chi ^{2}(\nu )}ثمX/ψinv-χ2(ν){\displaystyle X/\psi \sim {\mbox{inv-}}\chi ^{2}(\nu )}إلى جانبψ/Xχ2(ν){\displaystyle \psi /X\sim \chi ^{2}(\nu )}و1/Xψ-1χ2(ν){\displaystyle 1/X\sim \psi ^{-1}\chi ^{2}(\nu )}.

بدلاً منψ{\displaystyle \psi }ومع ذلك، فإن توزيع كاي تربيع العكسي المُقاس يتم تحديده في أغلب الأحيان بواسطة مُعامل المقياس.τ2=ψ/ν{\displaystyle \tau ^{2}=\psi /\nu }والتوزيعντ2inv-χ2(ν){\displaystyle \nu \tau ^{2}\,{\mbox{inv-}}\chi ^{2}(\nu )}يُرمز إليه بـالمقياس-الاستثمار-χ2(ν،τ2){\displaystyle {\mbox{Scale-inv-}}\chi ^{2}(\nu ,\tau ^{2})}.

من ناحيةτ2{\displaystyle \tau ^{2}}يمكن كتابة العلاقات المذكورة أعلاه على النحو التالي:

لوXالمقياس-الاستثمار-χ2(ν،τ2){\displaystyle X\sim {\mbox{Scale-inv-}}\chi ^{2}(\nu ,\tau ^{2})}ثمXντ2inv-χ2(ν){\displaystyle {\frac {X}{\nu \tau ^{2}}}\sim {\mbox{inv-}}\chi ^{2}(\nu )}إلى جانبντ2Xχ2(ν){\displaystyle {\frac {\nu \tau ^{2}}{X}}\sim \chi ^{2}(\nu )}و1/X1ντ2χ2(ν){\displaystyle 1/X\sim {\frac {1}{\nu \tau ^{2}}}\chi ^{2}(\nu )}.

تُعد هذه العائلة من توزيعات مربع كاي المعكوسة المُقاسة إعادة تحديد معلمات توزيع جاما المعكوس .

على وجه التحديد، إذا

Xψinv-χ2(ν)=المقياس-الاستثمار-χ2(ν،τ2){\displaystyle X\sim \psi \,{\mbox{inv-}}\chi ^{2}(\nu )={\mbox{Scale-inv-}}\chi ^{2}(\nu ,\tau ^{2})}  ثم  Xإنف-غاما(ν2،ψ2)=إنف-غاما(ν2،ντ22){\displaystyle X\sim {\textrm {Inv-Gamma}}\left({\frac {\nu }{2}},{\frac {\psi }{2}}\right)={\textrm {Inv-Gamma}}\left({\frac {\nu }{2}},{\frac {\nu \tau ^{2}}{2}}\right)}

يمكن استخدام أي من الشكلين لتمثيل توزيع الإنتروبيا القصوى لعزم معكوس أول ثابت(هـ(1/X)){\displaystyle (E(1/X))}والعزم اللوغاريتمي الأول(هـ(ln(X)){\displaystyle (E(\ln(X))}.

يُستخدم توزيع كاي تربيع العكسي المُقاس بشكل خاص في الإحصاء البايزي . تحديدًا، يُمكن استخدام هذا التوزيع كتوزيع احتمالي مُرافق لمعامل التباين في التوزيع الطبيعي . ويُعطى التوزيع الاحتمالي المُرافق نفسه في نماذج بديلة بواسطة توزيع جاما العكسي .

توصيف

تمتد دالة كثافة الاحتمال لتوزيع كاي تربيع العكسي المُقاس على المجالx>0{\displaystyle x>0}وهو

و(x؛ν،τ2)=(τ2ν/2)ν/2Γ(ν/2) خبرة[-ντ22x]x1+ν/2{\displaystyle 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}}}}

أينν{\displaystyle \nu }يمثل معامل درجات الحرية وτ2{\displaystyle \tau ^{2}}يمثل معامل المقياس . دالة التوزيع التراكمي هي

F(x؛ν،τ2)=Γ(ν2،τ2ν2x)/Γ(ν2){\displaystyle F(x;\nu ,\tau ^{2})=\Gamma \left({\frac {\nu }{2}},{\frac {\tau ^{2}\nu }{2x}}\right)\left/\Gamma \left({\frac {\nu }{2}}\right)\right.}
=سؤال(ν2،τ2ν2x){\displaystyle =Q\left({\frac {\nu }{2}},{\frac {\tau ^{2}\nu }{2x}}\right)}

أينΓ(أ،x){\displaystyle \Gamma (a,x)}دالة غاما غير الكاملة ،Γ(x){\displaystyle \Gamma (x)} دالة غاما وسؤال(أ،x){\displaystyle Q(a,x)}هي دالة غاما منتظمة . الدالة المميزة هي

φ(ت؛ν،τ2)={\displaystyle \varphi (t;\nu ,\tau ^{2})=}
2Γ(ν2)(-أناτ2νت2)ν4كν2(-2أناτ2νت)،{\displaystyle {\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}}\يمين)،}

أينكν2(z){\displaystyle K_{\frac {\nu }{2}}(z)}هي دالة بيسل المعدلة من النوع الثاني .

تقدير المعلمات

تقدير الاحتمال الأقصى لـτ2{\displaystyle \tau ^{2}}يكون

τ2=ن/أنا=1ن1xأنا.{\displaystyle \tau ^{2}=n/\sum _{i=1}^{n}{\frac {1}{x_{i}}}.}

تقدير الاحتمال الأقصى لـν2{\displaystyle {\frac {\nu }{2}}}يمكن إيجادها باستخدام طريقة نيوتن على:

ln(ν2)-ψ(ν2)=1نأنا=1نln(xأنا)-ln(τ2)،{\displaystyle \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),}

أينψ(x){\displaystyle \psi (x)}دالة ديغاما . يمكن إيجاد تقدير أولي لها بأخذ صيغة المتوسط ​​وحلها لإيجاد قيمة .ν.{\displaystyle \nu .} يتركx¯=1نأنا=1نxأنا{\displaystyle {\bar {x}}={\frac {1}{n}}\sum _{i=1}^{n}x_{i}}ليكن متوسط ​​العينة. ثم تقدير أولي لـν{\displaystyle \nu }يُعطى بواسطة:

ν2=x¯x¯-τ2.{\displaystyle {\frac {\nu }{2}}={\frac {\bar {x}}{{\bar {x}}-\tau ^{2}}}.}

التقدير البايزي لتباين التوزيع الطبيعي

للتوزيع العكسي المقياسي لتوزيع كاي تربيع تطبيق مهم ثانٍ، وهو التقدير البايزي لتباين التوزيع الطبيعي.

وفقًا لنظرية بايز ، فإن التوزيع الاحتمالي اللاحق للكميات ذات الأهمية يتناسب مع حاصل ضرب التوزيع المسبق للكميات ودالة الاحتمال :

ص(σ2|د،أنا)ص(σ2|أنا)ص(د|σ2){\displaystyle p(\sigma ^{2}|D,I)\propto p(\sigma ^{2}|I)\;p(D|\sigma ^{2})}

حيث يمثل D البيانات ويمثل I أي معلومات أولية حول σ 2 التي قد تكون لدينا بالفعل.

ينشأ السيناريو الأبسط إذا كان المتوسط ​​μ معروفًا بالفعل؛ أو، بدلاً من ذلك، إذا كان المطلوب هو التوزيع الشرطي لـ σ 2 ، لقيمة معينة مفترضة لـ μ .

عندئذٍ، يكون حد الاحتمال L ( σ 2 | D ) = p ( D | σ 2 ) بالشكل المألوف

ل(σ2|د،μ)=1(2πσ)نخبرة[-أنان(xأنا-μ)22σ2]{\displaystyle {\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/ σ² ، والذي يمكن اعتباره (على سبيل المثال، باتباع جيفريز ) أقل توزيع احتمالي مسبق ممكن لـ σ² في هذه المسألة، نحصل على احتمال لاحق مُجمّع

ص(σ2|د،أنا،μ)1σن+2خبرة[-أنان(xأنا-μ)22σ2]{\displaystyle 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 و τ 2 = s 2 = (1/ n ) Σ ( xi - μ ) 2

يشير جيلمان وزملاؤه إلى أن ظهور هذا التوزيع مرة أخرى، والذي شوهد سابقًا في سياق أخذ العينات، قد يبدو أمرًا لافتًا للنظر؛ ولكن بالنظر إلى اختيار التوزيع المسبق "فإن هذه النتيجة ليست مفاجئة". [ 1 ]

على وجه الخصوص، فإن اختيار توزيع احتمالي مسبق ثابت عند إعادة التحجيم لـ σ 2 يؤدي إلى أن احتمال نسبة σ 2 / s 2 له نفس الشكل (مستقل عن متغير التكييف) عند التكييف على s 2 كما هو الحال عند التكييف على σ 2 :

ص(σ2s2|s2)=ص(σ2s2|σ2){\displaystyle p({\tfrac {\sigma ^{2}}{s^{2}}}|s^{2})=p({\tfrac {\sigma ^{2}}{s^{2}}}|\sigma ^{2})}

في حالة نظرية أخذ العينات، المشروطة بـ σ 2 ، يكون توزيع الاحتمال لـ (1/s 2 ) هو توزيع كاي تربيع معكوس مُقاس؛ وبالتالي فإن توزيع الاحتمال لـ σ 2 المشروط بـ s 2 ، بالنظر إلى توزيع مسبق غير متأثر بالمقياس، هو أيضًا توزيع كاي تربيع معكوس مُقاس.

استخدمه كمرجع إعلامي مسبق

إذا كانت هناك معلومات أكثر حول القيم المحتملة لـ σ 2 ، فإن التوزيع من عائلة مربع كاي المعكوس المُقاس، مثل Scale-inv- χ 2 ( n 0 , s 0 2 ) يمكن أن يكون شكلاً مناسباً لتمثيل معلومات مسبقة أكثر دقة لـ σ 2 ، كما لو كان ذلك من نتيجة n 0 من الملاحظات السابقة (على الرغم من أن n 0 ليس بالضرورة أن يكون عددًا صحيحًا):

ص(σ2|أنا،μ)1σن0+2خبرة[-ن0s022σ2]{\displaystyle 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]}

سيؤدي هذا التوزيع الاحتمالي المسبق إلى التوزيع الاحتمالي اللاحق

ص(σ2|د،أنا،μ)1σن+ن0+2خبرة[-نs2+ن0s022σ2]{\displaystyle 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 )   const.، والذي يعطي التوزيع الاحتمالي اللاحق المشترك التالي لـ μ و σ 2 ،

ص(μ،σ2|د،أنا)1σن+2خبرة[-أنان(xأنا-μ)22σ2]=1σن+2خبرة[-أنان(xأنا-x¯)22σ2]خبرة[-ن(μ-x¯)22σ2]{\displaystyle {\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}}}

يتم الحصول على التوزيع الاحتمالي الهامشي اللاحق لـ σ 2 من التوزيع الاحتمالي المشترك اللاحق عن طريق التكامل على μ ،

ص(σ2|د،أنا)1σن+2خبرة[-أنان(xأنا-x¯)22σ2]-خبرة[-ن(μ-x¯)22σ2]دμ=1σن+2خبرة[-أنان(xأنا-x¯)22σ2]2πσ2/ن(σ2)-(ن+1)/2خبرة[-(ن-1)s22σ2]{\displaystyle {\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}}}

هذا مرة أخرى توزيع كاي تربيع معكوس مُقاس، مع معلماتن-1{\displaystyle \scriptstyle {n-1}\;}وs2=(xأنا-x¯)2/(ن-1){\displaystyle \scriptstyle {s^{2}=\sum (x_{i}-{\bar {x}})^{2}/(n-1)}}.

  • لوXالمقياس-الاستثمار-χ2(ν،τ2){\displaystyle X\sim {\mbox{Scale-inv-}}\chi ^{2}(\nu ,\tau ^{2})}ثمكXالمقياس-الاستثمار-χ2(ν،كτ2){\displaystyle kX\sim {\mbox{Scale-inv-}}\chi ^{2}(\nu ,k\tau ^{2})\,}
  • لوXinv-χ2(ν){\displaystyle X\sim {\mbox{inv-}}\chi ^{2}(\nu )\,}( توزيع مربع كاي العكسي ) إذنXالمقياس-الاستثمار-χ2(ν،1/ν){\displaystyle X\sim {\mbox{Scale-inv-}}\chi ^{2}(\nu ,1/\nu )\,}
  • لوXالمقياس-الاستثمار-χ2(ν،τ2){\displaystyle X\sim {\mbox{Scale-inv-}}\chi ^{2}(\nu ,\tau ^{2})}ثمXτ2νinv-χ2(ν){\displaystyle {\frac {X}{\tau ^{2}\nu }}\sim {\mbox{inv-}}\chi ^{2}(\nu )\,}( توزيع مربع كاي العكسي )
  • لوXالمقياس-الاستثمار-χ2(ν،τ2){\displaystyle X\sim {\mbox{Scale-inv-}}\chi ^{2}(\nu ,\tau ^{2})}ثمXإنف-غاما(ν2،ντ22){\displaystyle X\sim {\textrm {Inv-Gamma}}\left({\frac {\nu }{2}},{\frac {\nu \tau ^{2}}{2}}\right)}( توزيع غاما المعكوس )
  • يُعد توزيع مربع كاي العكسي المُقاس حالة خاصة من توزيع بيرسون من النوع 5

مراجع

  • جيلمان، أندرو؛ وآخرون  . (2014). تحليل البيانات البايزي (  الطبعة الثالثة). بوكا راتون: مطبعة سي آر سي. ص  583. ISBN 978-1-4398-4095-5.
  1. جيلمان، أندرو؛ وآخرون . (2014). تحليل البيانات البايزي ( الطبعة الثالثة). بوكا راتون: مطبعة سي آر سي. ص 65. ISBN    978-1-4398-4095-5.