توزيع غاما المعكوس

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

لعلّ الاستخدام الرئيسي لتوزيع غاما المعكوس يكمن في الإحصاء البايزي ، حيث يظهر هذا التوزيع كتوزيع احتمالي خلفي هامشي للتباين المجهول لتوزيع طبيعي ، في حال استخدام توزيع احتمالي مسبق غير مُعلِم ، وكتوزيع احتمالي مسبق مترافق قابل للتحليل ، في حال اشتراط استخدام توزيع احتمالي مسبق مُعلِم. [ 1 ] يشيع بين بعض البايزيين النظر في تمثيل بديل للتوزيع الطبيعي بدلالة الدقة ، المُعرَّفة بأنها مقلوب التباين، مما يسمح باستخدام توزيع غاما مباشرةً كتوزيع احتمالي مسبق مترافق. بينما يُفضِّل بايزيون آخرون تمثيل توزيع غاما المعكوس بشكل مختلف، كتوزيع كاي تربيع معكوس مُقاس .

توصيف

دالة كثافة الاحتمال

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

و(x؛α،β)=βαΓ(α)(1/x)α+1خبرة(-β/x){\displaystyle f(x;\alpha ,\beta )={\frac {\beta ^{\alpha }}{\Gamma (\alpha )}}(1/x)^{\alpha +1}\exp \left(-\beta /x\right)}

مع معلمة الشكلα{\displaystyle \alpha }ومعامل المقياسβ{\displaystyle \beta }[ 2 ] هناΓ(){\displaystyle \Gamma (\cdot )}يرمز إلى دالة جاما .

بخلاف توزيع جاما ، الذي يحتوي على حد أسي مشابه إلى حد ما،β{\displaystyle \beta }هو معامل مقياس حيث أن دالة الكثافة تحقق ما يلي:

و(x؛α،β)=و(x/β؛α،1)β{\displaystyle f(x;\alpha ,\beta )={\frac {f(x/\beta ;\alpha ,1)}{\beta }}}

دالة التوزيع التراكمي

دالة التوزيع التراكمي هي دالة غاما المنتظمة

F(x؛α،β)=Γ(α،βx)Γ(α)=سؤال(α،βx){\displaystyle F(x;\alpha ,\beta )={\frac {\Gamma \left(\alpha ,{\frac {\beta }{x}}\right)}{\Gamma (\alpha )}}=Q\left(\alpha ,{\frac {\beta }{x}}\right)\!}

حيث يمثل البسط دالة غاما غير الكاملة العليا، ويمثل المقام دالة غاما . تسمح العديد من حزم البرامج الرياضية بالحساب المباشر لـسؤال{\displaystyle Q}، دالة غاما المنتظمة.

لحظات

بشرط أنα>ن{\displaystyle \alpha >n}، الن{\displaystyle n}يتم إعطاء العزم رقم -th لتوزيع جاما العكسي بواسطة [ 3 ]

هـ[Xن]=βنΓ(α-ن)Γ(α)=βن(α-1)(α-ن).{\displaystyle \mathrm {E} [X^{n}]=\beta ^{n}{\frac {\Gamma (\alpha -n)}{\Gamma (\alpha )}}={\frac {\beta ^{n}}{(\alpha -1)\cdots (\alpha -n)}}.}

الدالة المميزة

يتميز توزيع غاما العكسي بدالة مميزة2(-أناβت)α2Γ(α)كα(-4أناβت){\displaystyle {\frac {2\left(-i\beta t\right)^{\!\!{\frac {\alpha }{2}}}}{\Gamma (\alpha )}}K_{\alpha }\left({\sqrt {-4i\beta t}}\right)}أينكα{\displaystyle K_{\alpha }}هي دالة بيسل المعدلة من النوع الثاني.

ملكيات

لα>0{\displaystyle \alpha >0}وβ>0{\displaystyle \beta >0}،

هـ[ln(X)]=ln(β)-ψ(α){\displaystyle \mathbb {E} [\ln(X)]=\ln(\beta )-\psi (\alpha )\,}

و

هـ[X-1]=αβ،{\displaystyle \mathbb {E} [X^{-1}]={\frac {\alpha }{\beta }},\,}

إنتروبيا المعلومات هي

ح(X)=هـ[-ln(ص(X))]=هـ[-αln(β)+ln(Γ(α))+(α+1)ln(X)+βX]=-αln(β)+ln(Γ(α))+(α+1)ln(β)-(α+1)ψ(α)+α=α+ln(βΓ(α))-(α+1)ψ(α).{\displaystyle {\begin{aligned}\operatorname {H} (X)&=\operatorname {E} [-\ln(p(X))]\\&=\operatorname {E} \left[-\alpha \ln(\beta )+\ln(\Gamma (\alpha ))+(\alpha +1)\ln(X)+{\frac {\beta }{X}}\right]\\&=-\alpha \ln(\beta )+\ln(\Gamma (\alpha ))+(\alpha +1)\ln(\beta )-(\alpha +1)\psi (\alpha )+\alpha \\&=\alpha +\ln(\beta \Gamma (\alpha ))-(\alpha +1)\psi (\alpha ).\end{aligned}}}

أين ψ(α){\displaystyle \psi (\alpha )}هي دالة ديغاما .

إن انحراف Kullback-Leibler لـ Inverse-Gamma( α p , β p ) من Inverse-Gamma( α q , β q ) هو نفس انحراف KL لـ Gamma ( α p , β p ) من Gamma ( α q , β q ) :

دكل(αص،βص؛αq،βq)=هـ[سجلρ(X)π(X)]=هـ[سجلρ(1/Y)π(1/Y)]=هـ[سجلρجي(Y)πجي(Y)]،{\displaystyle D_{\mathrm {KL} }(\alpha _{p},\beta _{p};\alpha _{q},\beta _{q})=\mathbb {E} \left[\log {\frac {\rho (X)}{\pi (X)}}\right]=\mathbb {E} \left[\log {\frac {\rho (1/Y)}{\pi (1/Y)}}\right]=\mathbb {E} \left[\log {\frac {\rho _{G}(Y)}{\pi _{G}(Y)}}\right],}

أين ρ،π{\displaystyle \rho ,\pi }هي دوال كثافة الاحتمال لتوزيعات غاما المعكوسة و ρجي،πجي{\displaystyle \rho _{G},\pi _{G}}هي دوال كثافة الاحتمال لتوزيعات جاما،Y{\displaystyle Y}يتم توزيع جاما( α p , β p ) .

دكل(αص،βص؛αq،βq)=(αص-αq)ψ(αص)-سجلΓ(αص)+سجلΓ(αq)+αq(سجلβص-سجلβq)+αصβq-βصβص.{\displaystyle {\begin{aligned}D_{\mathrm {KL} }(\alpha _{p},\beta _{p};\alpha _{q},\beta _{q})={}&(\alpha _{p}-\alpha _{q})\psi (\alpha _{p})-\log \Gamma (\alpha _{p})+\log \Gamma (\alpha _{q})+\alpha _{q}(\log \beta _{p}-\log \beta _{q})+\alpha _{p}{\frac {\beta _{q}-\beta _{p}}{\beta _{p}}}.\end{aligned}}}
  • لوXإنف-غاما(α،β){\displaystyle X\sim {\mbox{Inv-Gamma}}(\alpha ,\beta )}ثمكXإنف-غاما(α،كβ){\displaystyle kX\sim {\mbox{Inv-Gamma}}(\alpha ,k\beta )\,}، لك>0{\displaystyle k>0}
  • لوXإنف-غاما(α،12){\displaystyle X\sim {\mbox{Inv-Gamma}}(\alpha ,{\tfrac {1}{2}})}ثمXInv-χ2(2α){\displaystyle X\sim {\mbox{Inv-}}\chi ^{2}(2\alpha )\,}( توزيع مربع كاي العكسي )
  • لوXإنف-غاما(α2،12){\displaystyle X\sim {\mbox{Inv-Gamma}}({\tfrac {\alpha }{2}},{\tfrac {1}{2}})}ثمXانعكاس مُقاسχ2(α،1α){\displaystyle X\sim {\mbox{Scaled Inv-}}\chi ^{2}(\alpha ,{\tfrac {1}{\alpha }})\,}( توزيع مربع كاي المعكوس المقياس )
  • لوXإنف-غاما(12،ج2){\displaystyle X\sim {\textrm {Inv-Gamma}}({\tfrac {1}{2}},{\tfrac {c}{2}})}ثمXليفي(0،ج){\displaystyle X\sim {\textrm {Levy}}(0,c)\,}( توزيع ليفي )
  • لوXإنف-غاما(1،ج){\displaystyle X\sim {\textrm {Inv-Gamma}}(1,c)}ثم1Xتجربة(ج){\displaystyle {\tfrac {1}{X}}\sim {\textrm {Exp}}(c)\,}( التوزيع الأسي )
  • لوXجاما(α،β){\displaystyle X\sim {\mbox{Gamma}}(\alpha ,\beta )\,}( توزيع جاما مع معامل المعدل )β{\displaystyle \beta }) ثم1Xإنف-غاما(α،β){\displaystyle {\tfrac {1}{X}}\sim {\mbox{Inv-Gamma}}(\alpha ,\beta )\,}(انظر الاشتقاق في الفقرة التالية لمزيد من التفاصيل)
  • لاحظ أنه إذاXجاما(ك،θ){\displaystyle X\sim {\mbox{Gamma}}(k,\theta )}(توزيع جاما مع معامل المقياس)θ{\displaystyle \theta }) ثم1/Xإنف-غاما(ك،1/θ){\displaystyle 1/X\sim {\mbox{Inv-Gamma}}(k,1/\theta )}
  • يُعد توزيع غاما المعكوس حالة خاصة من توزيع بيرسون من النوع 5
  • يُعد توزيع ويشارت العكسي تعميمًا متعدد المتغيرات لتوزيع جاما العكسي .
  • للاطلاع على توزيع مجموع متغيرات جاما المعكوسة المستقلة، انظر ويتكوفسكي (2001).

الاشتقاق من توزيع جاما

يتركXجاما(α،β){\displaystyle X\sim {\mbox{Gamma}}(\alpha ,\beta )}وتذكر أن دالة كثافة الاحتمال لتوزيع جاما هي

وX(x)=βαΓ(α)xα-1هـ-βx{\displaystyle f_{X}(x)={\frac {\beta ^{\alpha }}{\Gamma (\alpha )}}x^{\alpha -1}e^{-\beta x}}،x>0{\displaystyle x>0}.

لاحظ أنβ{\displaystyle \beta }يمثل معامل المعدل من منظور توزيع جاما.

حدد التحويلY=ز(X)=1X{\displaystyle Y=g(X)={\tfrac {1}{X}}}ثم، ملف PDF الخاص بـY{\displaystyle Y}يكون

وY(y)=وX(ز-1(y))|ددyز-1(y)|=βαΓ(α)(1y)α-1خبرة(-βy)1y2=βαΓ(α)(1y)α+1خبرة(-βy)=βαΓ(α)y-α-1خبرة(-βy){\displaystyle {\begin{aligned}f_{Y}(y)&=f_{X}\left(g^{-1}(y)\right)\left|{\frac {d}{dy}}g^{-1}(y)\right|\\[6pt]&={\frac {\beta ^{\alpha }}{\Gamma (\alpha )}}\left({\frac {1}{y}}\right)^{\alpha -1}\exp \left({\frac {-\beta }{y}}\right){\frac {1}{y^{2}}}\\[6pt]&={\frac {\beta ^{\alpha }}{\Gamma (\alpha )}}\left({\frac {1}{y}}\right)^{\alpha +1}\exp \left({\frac {-\beta }{y}}\right)\\[6pt]&={\frac {\beta ^{\alpha }}{\Gamma (\alpha )}}y^{-\alpha -1}\exp \left({\frac {-\beta }{y}}\right)\\[6pt]\end{aligned}}}

لاحظ أنβ{\displaystyle {\beta }}يمثل معامل المقياس من منظور توزيع غاما العكسي. ويمكن إثبات ذلك بسهولة من خلال ملاحظة أنβ{\displaystyle {\beta }}يستوفي شروط كونه معامل مقياس .

و(y/β؛α،1)β=1β1Γ(α)(yβ)-α-1خبرة(-1y/β)=βαΓ(α)y-α-1خبرة(-βy)=و(y؛α،β){\displaystyle {\begin{aligned}{\frac {f(y/\beta ;\alpha ,1)}{\beta }}&={\frac {1}{\beta }}{\frac {1}{\Gamma (\alpha )}}\left({\frac {y}{\beta }}\right)^{-\alpha -1}\exp \left(-{\frac {1}{y/\beta }}\right)\\[6pt]&={\frac {\beta ^{\alpha }}{\Gamma (\alpha )}}y^{-\alpha -1}\exp \left(-{\frac {\beta }{y}}\right)\\[6pt]&=f(y;\alpha ,\beta )\end{aligned}}}

حدوث

انظر أيضاً

مراجع

  1. هوف، ب. (2009). "النموذج الطبيعي". دورة تمهيدية في الأساليب الإحصائية البايزية . سبرينغر. ص 67-88 . ISBN  978-0-387-92299-7.
  2. "توزيع غاما العكسي - وثائق لغة وولفرام" . reference.wolfram.com . تم الاطلاع عليه بتاريخ 9 أبريل 2018 .
  3. جون د. كوك (3 أكتوبر 2008). "توزيع غاما العكسي" (ملف PDF) . تم الاطلاع عليه في 3 ديسمبر 2018 .
  4. لودكوفسكي، مايك (2007). "ملاحظات حول الحركة البراونية في الرياضيات 526" (ملف PDF) . جامعة كاليفورنيا، سانتا باربرا. الصفحات 5-6 . مؤرشف من النسخة الأصلية (ملف PDF) بتاريخ 26 يناير 2022. تاريخ الاطلاع: 13 أبريل 2021 . 
  • ويتكوفسكي، ف. (2001). " حساب توزيع توليفة خطية من متغيرات غاما المعكوسة". كيبرنيتيكا . 37 (1): 79-90 . MR 1825758. Zbl 1263.62022 .