توزيع جاما

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

  1. مع معلمة الشكلα{\displaystyle \alpha }ومعامل المقياس θ
  2. مع معلمة الشكلα{\displaystyle \alpha }ومعامل معدلβ=1/θ{\displaystyle \beta =1/\theta }

في كل من هذه الأشكال، يكون كلا المعاملين أعدادًا حقيقية موجبة.

لهذا التوزيع تطبيقات مهمة في مجالات متنوعة، منها الاقتصاد القياسي ، والإحصاء البايزي ، واختبارات العمر. [ 3 ] في الاقتصاد القياسي، يشيع استخدام معلمات ( α ، θ ) لنمذجة أوقات الانتظار، مثل الوقت حتى الوفاة، حيث غالبًا ما يأخذ شكل توزيع إيرلانغ لقيم α الصحيحة . أما الإحصائيون البايزيون فيفضلون معلمات ( α ، β )، مستخدمين توزيع غاما كتوزيع مسبق مترافق لعدة معلمات ذات مقياس معكوس، مما يسهل التحليل في حسابات التوزيع اللاحق.

يُعد توزيع جاما توزيع الاحتمالية ذو الإنتروبيا القصوى (سواء بالنسبة إلى مقياس أساسي موحد أو1/x{\displaystyle 1/x}القياس الأساسي) لمتغير عشوائي X حيث E [ X ] = αθ = α / β ثابت وأكبر من الصفر، و E [ln X ] = ψ ( α ) + ln θ = ψ ( α ) − ln β ثابت ( ψ هي دالة digamma ). [ 4 ]

التعريفات

يبدو أن استخدام المعاملات α و θ أكثر شيوعًا في الاقتصاد القياسي والمجالات التطبيقية الأخرى، حيث يُستخدم توزيع غاما بكثرة لنمذجة أوقات الانتظار. على سبيل المثال، في اختبارات العمر ، يُعد وقت الانتظار حتى الوفاة متغيرًا عشوائيًا يُنمذج غالبًا باستخدام توزيع غاما. انظر هوغ وكريغ [ 5 ] للاطلاع على دوافع واضحة.

يُعدّ استخدام المعاملات α و β أكثر شيوعًا في الإحصاء البايزي ، [ 6 ] حيث يُستخدم توزيع غاما كتوزيع احتمالي مسبق مترافق لأنواع مختلفة من معاملات المقياس العكسي (المعدل)، مثل λ للتوزيع الأسي أو توزيع بواسون [ 7 ] ، أو حتى β لتوزيع غاما نفسه. ويُستخدم توزيع غاما العكسي، ذو الصلة الوثيقة ، كتوزيع احتمالي مسبق مترافق لمعاملات المقياس، مثل تباين التوزيع الطبيعي .

إذا كان α عددًا صحيحًا موجبًا ، فإن التوزيع يمثل توزيع إيرلانغ ؛ أي مجموع α من المتغيرات العشوائية المستقلة الموزعة أُسّيًا ، ولكل منها متوسط ​​θ .

التوصيف باستخدام الشكل α والمعدل β

يمكن تمثيل توزيع غاما بدلالة معامل الشكل α ومعامل المقياس العكسي β = 1/ θ ، والذي يُسمى معامل المعدل . ويُرمز للمتغير العشوائي X الذي يتبع توزيع غاما بمعامل شكل α ومعدل β بالرمز التالي:

XΓ(α،β)جاما(α،β){\displaystyle X\sim \Gamma (\alpha ,\beta )\equiv \operatorname {Gamma} (\alpha ,\beta )}

دالة كثافة الاحتمال المقابلة في معلمات الشكل والمعدل هي

و(x؛α،β)=xα-1هـ-βxβαΓ(α) ل x>0α،β>0،{\displaystyle {\begin{aligned}f(x;\alpha ,\beta )&={\frac {x^{\alpha -1}e^{-\beta x}\beta ^{\alpha }}{\Gamma (\alpha )}}\quad {\text{ لـ }}x>0\quad \alpha ,\beta >0,\\[6pt]\end{aligned}}}

أينΓ(α){\displaystyle \Gamma (\alpha )}دالة غاما . لجميع الأعداد الصحيحة الموجبة،Γ(α)=(α-1)!{\displaystyle \Gamma (\alpha )=(\alpha -1)!}.

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

F(x؛α،β)=0xو(u؛α،β)دu=γ(α،βx)Γ(α)،{\displaystyle F(x;\alpha ,\beta )=\int _{0}^{x}f(u;\alpha ,\beta )\,du={\frac {\gamma (\alpha ,\beta x)}{\Gamma (\alpha )}},}

أينγ(α،βx){\displaystyle \gamma (\alpha ,\beta x)}هي دالة غاما غير المكتملة الدنيا .

إذا كان α عددًا صحيحًا موجبًا (أي أن التوزيع هو توزيع إيرلانغ )، فإن دالة التوزيع التراكمي لها التوسع المتسلسل التالي: [ 8 ]

F(x؛α،β)=1-أنا=0α-1(βx)أناأنا!هـ-βx=هـ-βxأنا=α(βx)أناأنا!.{\displaystyle {\begin{aligned}F(x;\alpha ,\beta )&=1-\sum _{i=0}^{\alpha -1}{\frac {\left(\beta x\right)^{i}}{i!}}e^{-\beta x}\\[1ex]&=e^{-\beta x}\sum _{i=\alpha }^{\infty }{\frac {\left(\beta x\right)^{i}}{i!}}.\end{aligned}}}

التوصيف باستخدام الشكل α والمقياس θ

يُرمز إلى المتغير العشوائي X الذي يتبع توزيع جاما بشكل α ومقياس θ بالرمز التالي:

XΓ(α،θ)جاما(α،θ){\displaystyle X\sim \Gamma (\alpha ,\theta )\equiv \operatorname {Gamma} (\alpha ,\theta )}

توضيح لدالة كثافة الاحتمال (PDF) لغاما لقيم المعاملات α و x مع ضبط θ على 1، 2، 3، 4، 5، و 6. يمكن رؤية كل طبقة θ على حدة هنا . وكذلك بواسطة α و x . .

دالة كثافة الاحتمال باستخدام معلمات الشكل والمقياس هي

و(x؛α،θ)=xα-1هـ-x/θθαΓ(α) ل x>0 و α،θ>0.{\displaystyle f(x;\alpha ,\theta )={\frac {x^{\alpha -1}e^{-x/\theta }}{\theta ^{\alpha }\Gamma (\alpha )}}\quad {\text{ for }}x>0{\text{ and }}\alpha ,\theta >0.}

هنا Γ( α ) هي دالة جاما المحسوبة عند α .

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

F(x؛α،θ)=0xو(u؛α،θ)دu=γ(α،xθ)Γ(α)،{\displaystyle F(x;\alpha ,\theta )=\int _{0}^{x}f(u;\alpha ,\theta )\,du={\frac {\gamma {\left(\alpha ,{\frac {x}{\theta }}\right)}}{\Gamma (\alpha )}},}

أينγ(α،xθ){\textstyle \gamma {\left(\alpha ,{\frac {x}{\theta }}\right)}}هي دالة غاما غير المكتملة الدنيا .

ويمكن التعبير عنها أيضًا على النحو التالي، إذا كان α عددًا صحيحًا موجبًا (أي أن التوزيع هو توزيع إيرلانغ ): [ 8 ]

F(x؛α،θ)=1-أنا=0α-11أنا!(xθ)أناهـ-x/θ=هـ-x/θأنا=α1أنا!(xθ)أنا.{\displaystyle F(x;\alpha ,\theta )=1-\sum _{i=0}^{\alpha -1}{\frac {1}{i!}}\left({\frac {x}{\theta }}\right)^{i}e^{-x/\theta }=e^{-x/\theta }\sum _{i=\alpha }^{\infty }{\frac {1}{i!}}\left({\frac {x}{\theta }}\right)^{i}.}

كلا النموذجين شائعان لأن أياً منهما قد يكون أكثر ملاءمة حسب الموقف.

ملكيات

المتوسط ​​والتباين

يُعطى متوسط ​​توزيع جاما بضرب معاملي الشكل والمقياس: μ=αθ=α/β{\displaystyle \mu =\alpha \theta =\alpha /\beta } التباين هو: σ2=αθ2=α/β2{\displaystyle \sigma ^{2}=\alpha \theta ^{2}=\alpha /\beta ^{2}} الجذر التربيعي لمعامل الشكل العكسي يعطي معامل التباين : σ/μ=α-0.5=1/α{\displaystyle \sigma /\mu =\alpha ^{-0.5}=1/{\sqrt {\alpha }}}

الانحراف

يعتمد انحراف توزيع جاما فقط على معامل شكله، α ، وهو يساوي2/α.{\displaystyle 2/{\sqrt {\alpha }}.}

لحظات أسمى

العزم الخام رقم r يُعطى بالصيغة التالية:

هـ[Xر]=θرΓ(α+ر)Γ(α)=θرαر¯{\displaystyle \mathrm {E} [X^{r}]=\theta ^{r}{\frac {\Gamma (\alpha +r)}{\Gamma (\alpha )}}=\theta ^{r}\alpha ^{\overline {r}}}

معαر¯{\displaystyle \alpha ^{\overline {r}}}العامل التصاعدي .

التراكمات

يُعطى التراكم من الرتبة r بالصيغة التالية :

κر=θرα(ر-1)!=θرαΓ(ر){\displaystyle \kappa _{r}=\theta ^{r}\alpha (r-1)!=\theta ^{r}\alpha \Gamma (r)}.

التقريبات والحدود الوسيطة

حدود وتقريبات تقاربية لوسيط توزيع جاما. تشير المنطقة ذات اللون السماوي إلى الفجوة الكبيرة بين الحدود الدنيا والعليا المنشورة قبل عام 2021.

بخلاف المنوال والمتوسط، اللذين لهما صيغ حسابية سهلة تعتمد على المعلمات، فإن الوسيط لا يملك معادلة مغلقة. الوسيط لهذا التوزيع هو القيمةν{\displaystyle \nu }بحيث 1Γ(α)θα0νxα-1هـ-x/θدx=12.{\displaystyle {\frac {1}{\Gamma (\alpha )\theta ^{\alpha }}}\int _{0}^{\nu}x^{\alpha -1}e^{-x/\theta }dx={\frac {1}{2}}.}

قام تشين وروبين بمعالجة دقيقة لمشكلة تحديد التوسع التقاربي وحدود الوسيط لتوزيع جاما، حيث أثبتا أن (لـθ=1{\displaystyle \theta =1}) α-13<ν(α)<α،{\displaystyle \alpha -{\tfrac {1}{3}}<\nu (\alpha )<\alpha ,} أينμ(α)=α{\displaystyle \mu (\alpha )=\alpha }هو المتوسط ​​وν(α){\displaystyle \nu (\alpha )}هو الوسيط لـجاما(α،1){\displaystyle {\text{Gamma}}(\alpha ,1)}التوزيع. [ 9 ] بالنسبة للقيم الأخرى لمعامل المقياس، يتناسب المتوسط ​​معμ=αθ{\displaystyle \mu =\alpha \theta }، وسيتم قياس الحدود والتقريبات الوسيطة بشكل مماثل بواسطة θ .

وجد كي بي تشوي الحدود الخمسة الأولى في تقريب لوران التقاربي للوسيط من خلال مقارنة الوسيط بـ رامانوجانθ{\displaystyle \theta }دالة . [ 10 ] وجد بيرج وبيدرسن المزيد من المصطلحات: [ 11 ]ν(α)=α-13+8405α-1+18425515α-2+22483444525α-3-1900640815345358875α-4-يا(α-5)+{\displaystyle {\begin{aligned}\nu (\alpha )=\alpha &-{\frac {1}{3}}+{\frac {8}{405}}\alpha ^{-1}+{\frac {184}{25\,515}}\alpha ^{-2}+{\frac {2248}{3\,444\,525}}\alpha ^{-3}\\[1ex]&-{\frac {19\,006\,408}{15\,345\,358\,875}}\alpha ^{-4}-{\mathcal {O}}{\left(\alpha ^{-5}\right)}+\cdots \end{aligned}}}

خطان تقاربيان للوسيط لتوزيع جاما، تم إثبات أنهما حدود (الخط الأحمر المتصل العلوي والخط الأحمر المتقطع السفلي)، منν(α)2-1/α(أ+α){\displaystyle \nu (\alpha )\approx 2^{-1/\alpha }(A+\alpha )}، وتقريب بينهما يُعطي تقريبًا (باللون الأحمر المنقط) دقيقًا عند α = 1 ، ويبلغ أقصى خطأ نسبي له حوالي 0.6%. المنطقة المظللة باللون السماوي هي الفجوة المتبقية بين الحدود العليا والسفلى (أو الحدود المفترضة)، بما في ذلك هذه الحدود الجديدة والحدود الواردة في الشكل السابق.
رسم بياني لوغاريتمي للحدود العليا (خط متصل) والسفلى (خط متقطع) لوسيط توزيع غاما والفجوات بينهما. تمثل المناطق الخضراء والصفراء والسماوية الفجوة قبل ورقة ليون البحثية لعام ٢٠٢١. تعمل المناطق الخضراء والصفراء على تضييق هذه الفجوة بالحدود الدنيا التي أثبتها ليون. كما تعمل حدود ليون التي تم إثباتها في عام ٢٠٢٣ على تضييق المنطقة الصفراء بشكل أكبر. ضمن المنطقة الصفراء في الغالب، تم رسم الحدود التخمينية المُستنتجة باستخدام دالة كسرية مغلقة الصيغة، إلى جانب قيمة الوسيط المحسوبة عدديًا (خط منقط). توجد حدود مُستنتجة أدق، لكنها غير مرسومة، لأنها لن تُحل على هذا النطاق.

تعتبر المجاميع الجزئية لهذه السلاسل تقريبات جيدة لقيم α العالية بما فيه الكفاية ؛ فهي غير مرسومة في الشكل، الذي يركز على منطقة α المنخفضة التي يتم تقريبها بشكل أقل جودة.

أثبت بيرغ وبيدرسن أيضًا العديد من خصائص الوسيط، موضحين أنه دالة محدبة لـ α ، [ 12 ] وأن السلوك التقاربي بالقرب منα=0{\displaystyle \alpha =0}يكونν(α)هـ-γ2-1/α{\displaystyle \nu (\alpha )\approx e^{-\gamma }2^{-1/\alpha }}(حيث γ هو ثابت أويلر-ماسكيروني )، وذلك لجميعα>0{\displaystyle \alpha >0}الوسيط محدود بـα2-1/α<ν(α)<كهـ-1/3ك{\displaystyle \alpha 2^{-1/\alpha }<\nu (\alpha )<ke^{-1/3k}}[ 11 ]

حد أعلى خطي أدق، لـα1{\displaystyle \alpha \geq 1}تم تقديمها فقط في عام 2021 بواسطة غونت وميركل، [ 13 ] بالاعتماد على نتيجة بيرغ وبيدرسن التي مفادها أن ميلν(α){\displaystyle \nu (\alpha )}يكون في كل مكان أقل من 1: ν(α)α-1+سجل2  {\displaystyle \nu (\alpha )\leq \alpha -1+\log 2~~}لα1{\displaystyle \alpha \geq 1}(مع المساواة في)α=1{\displaystyle \alpha =1}) والتي يمكن توسيعها لتشمل حدًا لجميعα>0{\displaystyle \alpha >0}بأخذ القيمة القصوى مع الوتر الموضح في الشكل، حيث ثبت أن الوسيط محدب. [ 12 ]

تقريب للوسيط يكون دقيقًا تقاربيًا عند قيم α العالية ومعقولًا حتىα=0.5{\displaystyle \alpha =0.5}أو أقل قليلاً، ويتبع ذلك من تحويل ويلسون-هيلفرتي : ν(α)=α(1-19α)3{\displaystyle \nu (\alpha )=\alpha \left(1-{\frac {1}{9\alpha }}\right)^{3}} والتي تصبح سالبة لـα<1/9{\displaystyle \alpha <1/9}.

في عام 2021، اقترح ليون عدة تقريبات من الشكلν(α)2-1/α(أ+بα){\displaystyle \nu (\alpha )\approx 2^{-1/\alpha }(A+B\alpha )}لقد افترض قيمًا لـ A و B التي يكون عندها هذا التقريب حدًا أعلى أو أدنى محكمًا تقاربيًا لجميعα>0{\displaystyle \alpha >0}[ 14 ] على وجه الخصوص ، اقترح هذه الحدود المغلقة الشكل، والتي أثبتها في عام 2023: [ 15 ]

νل(α)=2-1/α(سجل2-13+α){\displaystyle \nu _{L\infty }(\alpha )=2^{-1/\alpha }\left(\log 2-{\tfrac {1}{3}}+\alpha \right)}يمثل حدًا أدنى، محكمًا تقاربًا كماα{\displaystyle \alpha \to \infty }νيو(α)=2-1/α(هـ-γ+α){\displaystyle \nu _{U}(\alpha )=2^{-1/\alpha }(e^{-\gamma }+\alpha )\quad }يمثل حدًا أعلى، محكمًا تقاربًا كماα0{\displaystyle \alpha \to 0}

أظهر ليون أيضًا (بشكل غير رسمي في عام 2021، وبشكل دقيق في عام 2023) حدين أدنى آخرين ليسا تعبيرين مغلقين ، بما في ذلك هذا الحد الذي يتضمن دالة غاما ، استنادًا إلى حل التعبير التكاملي باستبدال 1 بـهـ-x{\displaystyle e^{-x}}: ν(α)>(2Γ(α+1))-1/α{\displaystyle \nu (\alpha )>\left({\frac {2}{\Gamma (\alpha +1)}}\right)^{-1/\alpha }}(الاقتراب من المساواة كـك0{\displaystyle k\to 0}) وخط المماس عندα=1{\displaystyle \alpha =1}حيث تم العثور على المشتقةν(1)0.9680448{\displaystyle \nu ^{\prime }(1)\approx 0.9680448}: ν(α)ν(1)+(α-1)ν(1){\displaystyle \nu (\alpha )\geq \nu (1)+(\alpha -1)\nu ^{\prime }(1)\quad }(مع المساواة في)ك=1{\displaystyle k=1}) ν(α)سجل2+(α-1)[γ-2إي(-سجل2)-سجلسجل2]{\displaystyle \nu (\alpha )\geq \log 2+(\alpha -1)\left[\gamma -2\operatorname {Ei} (-\log 2)-\log \log 2\right]} حيث Ei هو التكامل الأسي . [ 14 ] [ 15 ]

بالإضافة إلى ذلك، أظهر أن عمليات الاستيفاء بين الحدود يمكن أن توفر تقريبات ممتازة أو حدودًا أدق للوسيط، بما في ذلك تقريب دقيق عندα=1{\displaystyle \alpha =1}(أينν(1)=سجل2{\displaystyle \nu (1)=\log 2}ويبلغ الحد الأقصى للخطأ النسبي أقل من 0.6%. جميع التقريبات والحدود المُستكملة هي من الشكل التالي: ν(α)ز~(α)νل(α)+(1-ز~(α))νيو(α){\displaystyle \nu (\alpha )\approx {\tilde {g}}(\alpha )\nu _{L\infty }(\alpha )+(1-{\tilde {g}}(\alpha ))\nu _{U}(\alpha )} أينز~{\displaystyle {\tilde {g}}}هي دالة استيفاء تعمل بشكل رتيب من 0 عند قيم α المنخفضة إلى 1 عند قيم α العالية ، وهي تقارب دالة استيفاء مثالية أو دقيقة.ز(α){\displaystyle g(\alpha )}: ز(α)=νيو(α)-ν(α)νيو(α)-νل(α){\displaystyle g(\alpha )={\frac {\nu _{U}(\alpha )-\nu (\alpha )}{\nu _{U}(\alpha )-\nu _{L\infty }(\alpha )}}} بالنسبة لأبسط دالة استيفاء تم النظر فيها، دالة كسرية من الدرجة الأولى ز~1(α)=αب0+α{\displaystyle {\tilde {g}}_{1}(\alpha )={\frac {\alpha }{b_{0}+\alpha }}} أضيق حد أدنى لديه ب0=8405+هـ-γسجل2-سجل222هـ-γ-سجل2+13-سجل20.143472{\displaystyle b_{0}={\frac {{\frac {8}{405}}+e^{-\gamma }\log 2-{\frac {\log ^{2}2}{2}}}{e^{-\gamma }-\log 2+{\frac {1}{3}}}}-\log 2\approx 0.143472} وأضيق حد أعلى لديه ب0=هـ-γ-سجل2+131-هـ-γπ2120.374654{\displaystyle b_{0}={\frac {e^{-\gamma }-\log 2+{\frac {1}{3}}}{1-{\frac {e^{-\gamma }\pi ^{2}}{12}}}}\approx 0.374654} تم رسم الحدود المُستكمَلة (معظمها داخل المنطقة الصفراء) في الرسم البياني اللوغاريتمي الموضح. تتوفر حدود أدق باستخدام دوال استيفاء مختلفة، ولكن ليس عادةً بمعاملات مغلقة الشكل كهذه. [ 14 ]

الخلاصة

إذا كان للمتغير العشوائي Xᵢ توزيع جاما ( αᵢ , θ ) حيث i = 1, 2, ..., N (أي أن جميع التوزيعات لها نفس معامل المقياس θ )، فإن

أنا=1شمالXأناجيأممأ(أنا=1شمالαأنا،θ){\displaystyle \sum _{i=1}^{N}X_{i}\sim \mathrm {Gamma} \left(\sum _{i=1}^{N}\alpha _{i},\theta \right)}

بشرط أن تكون جميع المتغيرات X i مستقلة .

بالنسبة للحالات التي تكون فيها المتغيرات X i مستقلة ولكن لها معلمات قياس مختلفة، انظر ماثاي [ 16 ] أو موشوبولوس [ 17 ] .

يُظهر توزيع جاما قابلية قسمة لا نهائية .

التوسع

لو Xجيأممأ(α،θ)،{\displaystyle X\sim \mathrm {Gamma} (\alpha ,\theta ),}

ثم، لأي قيمة لـ c > 0 ،

جXجيأممأ(α،جθ)،{\displaystyle cX\sim \mathrm {Gamma} (\alpha ,c\,\theta ),}باستخدام دوال توليد العزوم،

أو ما يعادل ذلك، إذا

Xجيأممأ(α،β){\displaystyle X\sim \mathrm {Gamma} \left(\alpha ,\beta \right)}(تحديد معلمات الشكل والمعدل)

جXجيأممأ(α،βج)،{\displaystyle cX\sim \mathrm {Gamma} \left(\alpha ,{\frac {\beta }{c}}\right),}

في الواقع، نعلم أنه إذا كان X متغيرًا عشوائيًا أسيًا بمعدل β ، فإن cX هو متغير عشوائي أسي بمعدل β / c ؛ وينطبق الشيء نفسه على متغيرات جاما (ويمكن التحقق من ذلك باستخدام دالة توليد العزوم ، انظر، على سبيل المثال، هذه الملاحظات ، 10.4-(ii)): الضرب في ثابت موجب c يقسم المعدل (أو، بشكل مكافئ، يضرب المقياس).

عائلة الدوال الأسية

توزيع جاما هو عائلة أسية ذات معلمتين مع معلمات طبيعية α − 1 و −1/ θ (بشكل مكافئ، α1 و βوإحصائيات طبيعية X و ln X.

إذا تم تثبيت معلمة الشكل α ، فإن عائلة التوزيعات ذات المعلمة الواحدة الناتجة هي عائلة أسية طبيعية .

التوقع والتباين اللوغاريتميان

يمكن للمرء أن يثبت ذلك

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

أو ما يعادل ذلك،

هـ[lnX]=ψ(α)+lnθ{\displaystyle \operatorname {E} [\ln X]=\psi (\alpha )+\ln \theta }

حيث ψ هي دالة ديغاما . وبالمثل،

متغير[lnX]=ψ(1)(α){\displaystyle \operatorname {var} [\ln X]=\psi ^{(1)}(\alpha )}

أينψ(1){\displaystyle \psi ^{(1)}}هي دالة تريغاما .

يمكن اشتقاق ذلك باستخدام صيغة العائلة الأسية لدالة توليد العزوم للإحصائية الكافية ، لأن إحدى الإحصائيات الكافية لتوزيع جاما هي ln x .

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

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

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

في تمثيل المعاملات α و θ ، تُعطى إنتروبيا المعلومات بواسطة

ح(X)=α+lnθ+lnΓ(α)+(1-α)ψ(α).{\displaystyle \operatorname {H} (X)=\alpha +\ln \theta +\ln \Gamma (\alpha )+(1-\alpha )\psi (\alpha ).}

تباعد كولباك-لايبير

توضيح لتباعد كولباك -لايبير (KL) لدالتين احتماليتين من نوع جاما. هنا ، β = β₀ + 1، والتي تم ضبطها على 1، 2، 3، 4، 5، و 6 . يظهر بوضوح عدم التناظر النموذجي لتباعد كولباك-لايبير. 

يتم إعطاء تباعد كولباك-لايبير (KL-divergence) لتوزيع جاما ( α p , β p ) ("التوزيع الحقيقي") من توزيع جاما ( α q , β q ) ("التوزيع التقريبي") بواسطة [ 18 ]

دكل(αص،βص؛αq،βq)=(αص-αq)ψ(αص)-سجلΓ(αص)Γ(αq)+αqسجلβصβq+αص(βqβص-1).{\displaystyle {\begin{aligned}D_{\mathrm {KL} }(\alpha _{p},\beta _{p};\alpha _{q},\beta _{q})={}&(\alpha _{p}-\alpha _{q})\psi (\alpha _{p})-\log {\frac {\Gamma (\alpha _{p})}{\Gamma (\alpha _{q})}}\\&{}+\alpha _{q}\log {\frac {\beta _{p}}{\beta _{q}}}+\alpha _{p}\left({\frac {\beta _{q}}{\beta _{p}}}-1\right).\end{aligned}}}

مكتوبًا باستخدام المعلمات α , θ ، يتم إعطاء اختلاف KL لجاما ( α p , θ p ) من Gamma ( α q , θ q ) بواسطة

دكل(αص،θص؛αq،θq)=(αص-αq)ψ(αص)-سجلΓ(αص)Γ(αq)+αqسجلθqθص+αص(θصθq-1).{\displaystyle {\begin{aligned}D_{\mathrm {KL} }(\alpha _{p},\theta _{p};\alpha _{q},\theta _{q})={}&(\alpha _{p}-\alpha _{q})\psi (\alpha _{p})-\log {\frac {\Gamma (\alpha _{p})}{\Gamma (\alpha _{q})}}\\&{}+\alpha _{q}\log {\frac {\theta _{q}}{\theta _{p}}}+\alpha _{p}\left({\frac {\theta _{p}}{\theta _{q}}}-1\right).\end{aligned}}}

تحويل لابلاس

تحويل لابلاس لدالة كثافة الاحتمال لتوزيع جاما، وهي دالة توليد العزوم لتوزيع جاما، هو

F(s)=هـ[هـ-sX]=1(1+θs)α=(ββ+s)α{\displaystyle F(s)=\operatorname {E} \left[e^{-sX}\right]={\frac {1}{\left(1+\theta s\right)^{\alpha }}}=\left({\frac {\beta }{\beta +s}}\right)^{\alpha }}

(أينX{\textstyle X}(متغير عشوائي ذو هذا التوزيع).

العرض الكامل عند نصف الحد الأقصى

يُعرَّف العرض الكامل عند نصف القيمة القصوى لتوزيع احتمالي أحادي النمط بأنه الفرق بين قيمتي الدعم اللتين تُنتجان قيمة تساوي نصف القيمة القصوى للدالة عند تقييمها بواسطة دالة كثافة الاحتمال للتوزيع . بالنسبة لتوزيع جاما ذي معامل الشكلα>1{\displaystyle \alpha >1}، قيمة FWHM هي [ 19 ]

Fدبليوحم=(α-1)θ[دبليو0(-2-1/(α-1)هـ)-دبليو-1(-2-1/(α-1)هـ)]{\displaystyle \mathrm {FWHM} =(\alpha -1)\theta \left[W_{0}\left(-{\frac {2^{-1/(\alpha -1)}}{e}}\right)-W_{-1}\left(-{\frac {2^{-1/(\alpha -1)}}{e}}\right)\right]}

متىα=1{\displaystyle \alpha =1}يتحول توزيع جاما إلى توزيع أسي ، ويكون عرضه عند نصف الارتفاع الأقصى (FWHM) هو

Fدبليوحم=θln2{\displaystyle \mathrm {FWHM} =\theta \ln {2}}

متىα<1{\displaystyle \alpha <1}، توزيع جاما ليس له نمط نهائي وبالتالي ليس له عرض نصف القيمة القصوى (FWHM) محدد جيدًا.

عام

  • يتركX1،X2،...،Xن{\displaystyle X_{1},X_{2},\ldots ,X_{n}}يكونن{\displaystyle n}المتغيرات العشوائية المستقلة والمتطابقة التوزيع التي تتبع التوزيع الأسي بمعامل معدل β ، إذنأناXأناجاما(ن،β){\textstyle \sum _{i}X_{i}\sim \operatorname {Gamma} (n,\beta )}حيث n هو معامل الشكل و β هو المعدل، وX¯=1نأناXأناجاما(ن،نβ){\textstyle {\bar {X}}={\frac {1}{n}}\sum _{i}X_{i}\sim \operatorname {Gamma} (n,n\beta )}.
  • إذا كان X يتبع توزيع غاما (1، β ) (في نموذج الشكل-المعدل)، فإن X يتبع توزيعًا أسيًا بمعامل معدل β . أما في نموذج الشكل-المقياس، فإن X يتبع توزيع غاما (1، θ ) بمعامل معدل 1/ θ .
  • إذا كان X يتبع توزيع غاما ( ν /2, 2) (في نموذج الشكل والمقياس)، فإن X يكون مطابقًا لتوزيع كاي تربيع ( χ² ( ν )) ذي درجات الحرية ν . وبالعكس، إذا كان Q يتبع توزيع كاي تربيع (χ² ( ν )) وكان c ثابتًا موجبًا، فإن cQ يتبع توزيع غاما ( ν /2, 2c ) .
  • إذا كانت θ = 1/ α ، فإن المرء يحصل على توزيع شولز-زيم ، والذي يستخدم بشكل بارز لنمذجة أطوال سلسلة البوليمر.
  • إذا كان α عددًا صحيحًا ، فإن توزيع جاما هو توزيع إيرلانغ ، وهو التوزيع الاحتمالي لوقت الانتظار حتى "وصول" الحدث رقم α في عملية بواسون أحادية البعد بكثافة 1/ θ .
XΓ(αZ،θ)،Yبويس(xθ)،{\displaystyle X\sim \Gamma (\alpha \in \mathbb {Z} ,\theta ),\qquad Y\sim \operatorname {Pois} \left({\frac {x}{\theta }}\right),}
ثم
برو(X>x)=برو(Y<α).{\displaystyle \Pr(X>x)=\Pr(Y<\alpha ).}
X2Γ(32،2أ2).{\displaystyle X^{2}\sim \Gamma {\left({\tfrac {3}{2}},2a^{2}\right)}.}
  • إذا كان X ~ Gamma( α , θ ) ، فإنخبرةX{\textstyle \exp X}يتبع توزيع لوغاريتمي غاما. [ 20 ]
  • إذا كان X ~ Gamma( α , θ ) ، فإنسجلX{\textstyle \log X}يتبع توزيعًا أُسّيًا غاما (يُختصر إلى exp-gamma). [ 21 ] ويُشار إليه أحيانًا بشكل خاطئ باسم توزيع لوغاريتم غاما. [ 22 ] توجد صيغ حساب متوسطه وتباينه في قسم #التوقع والتباين اللوغاريتمي .
  • إذا كان X ~ Gamma( α , θ ) ، فإنX{\displaystyle {\sqrt {X}}}يتبع توزيع غاما المعمم بمعاملات p = 2 و d = وأ=θ{\displaystyle a={\sqrt {\theta }}}.
  • وبشكل أعم، إذا كان X ~ Gamma( α , θ ) ، فإنXq{\displaystyle X^{q}}لq>0{\displaystyle q>0}يتبع توزيع غاما المعمم بمعاملات p = 1/ q و d = α / q وأ=θq{\displaystyle a=\theta ^{q}}.
  • إذا كان X ~ Gamma( α , θ ) بالشكل α والمقياس θ , ثم 1/ X ~ Inv-Gamma( α , θ −1 ) (انظر توزيع جاما العكسي للاشتقاق).
  • المعامل 1: إذاXكΓ(αك،θك){\displaystyle X_{k}\sim \Gamma (\alpha _{k},\theta _{k})\,}إذا كانوا مستقلين، فإنهم كذلك.α2θ2X1α1θ1X2F(2α1،2α2){\displaystyle {\frac {\alpha _{2}\theta _{2}X_{1}}{\alpha _{1}\theta _{1}X_{2}}}\sim \mathrm {F} (2\alpha _{1},2\alpha _{2})}أو ما يعادل ذلك،X1X2β(α1،α2،1،θ1θ2){\displaystyle {\frac {X_{1}}{X_{2}}}\sim \beta '\left(\alpha _{1},\alpha _{2},1,{\frac {\theta _{1}}{\theta _{2}}}\right)}
  • المعامل الثاني: إذاXكΓ(αك،βك){\displaystyle X_{k}\sim \Gamma (\alpha _{k},\beta _{k})\,}إذا كانوا مستقلين، فإنهم كذلك.α2β1X1α1β2X2F(2α1،2α2){\displaystyle {\frac {\alpha _{2}\beta _{1}X_{1}}{\alpha _{1}\beta _{2}X_{2}}}\sim \mathrm {F} (2\alpha _{1},2\alpha _{2})}أو ما يعادل ذلك،X1X2β(α1،α2،1،β2β1){\displaystyle {\frac {X_{1}}{X_{2}}}\sim \beta '\left(\alpha _{1},\alpha _{2},1,{\frac {\beta _{2}}{\beta _{1}}}\right)}
  • إذا كان X ~ Gamma( α , θ ) و Y ~ Gamma( β , θ ) موزعين بشكل مستقل، فإن X /( X + Y ) له توزيع بيتا بمعاملات α و β ، و X /( X + Y ) مستقل عن X + Y ، الذي هو توزيع Gamma( α + β , θ ) .
  • لوXنبيتا(α،نβ){\displaystyle X_{n}\sim {\text{Beta}}(\alpha ,n\beta )\,}وYن=نXن{\displaystyle Y_{n}=nX_{n}}، ثمYن{\displaystyle Y_{n}}يتقارب في التوزيع إلىجاما(α،β){\displaystyle {\text{Gamma}}(\alpha ,\beta )}تم تحديدها بموجب المعلمة 2.
  • إذا كانت X i ~ Gamma( α i , 1) موزعة بشكل مستقل، فإن المتجه ( X 1 / S , ..., X n / S ) ، حيث S = X 1 + ... + X n ، يتبع توزيع Dirichlet مع المعلمات α 1 , ..., α n .
  • بالنسبة إلى α الكبيرة ، يتقارب توزيع جاما مع التوزيع الطبيعي بمتوسط ​​μ = αθ والتباين σ 2 = αθ 2 .
  • يُعد توزيع جاما التوزيع المسبق المترافق لدقة التوزيع الطبيعي ذي المتوسط ​​المعروف .
  • يُعد توزيع جاما المصفوفي وتوزيع ويشارت تعميمات متعددة المتغيرات لتوزيع جاما (العينات عبارة عن مصفوفات موجبة محددة بدلاً من أعداد حقيقية موجبة).
  • يُعد توزيع جاما حالة خاصة من توزيع جاما المعمم ، وتوزيع جاما الصحيح المعمم ، وتوزيع جاوس المعكوس المعمم .
  • من بين التوزيعات المنفصلة، ​​يعتبر التوزيع ذو الحدين السالب أحيانًا النظير المنفصل لتوزيع جاما.
  • توزيعات تويدي - توزيع جاما هو أحد أفراد عائلة نماذج التشتت الأسي لتويدي .
  • يُعد توزيع جاما أحد أفراد عائلة التوزيعات نصف الطبيعية المعدلة .
  • لمعاملات الشكل والمقياسx|θΓ(α،θ){\displaystyle x|\theta \sim \Gamma (\alpha ,\theta )}، إذا كان معامل المقياسθأناجي(ب،1){\displaystyle \theta \sim IG(b,1)}أينأناجي{\displaystyle IG}يشير إلى توزيع غاما المعكوس ، ثم التوزيع الهامشيxβ(α،ب){\displaystyle x\sim \beta '(\alpha ,b)}أينβ{\displaystyle \beta '}يشير إلى توزيع بيتا الأولي .

مركب جاما

إذا كانت معلمة شكل توزيع غاما معروفة، ولكن معلمة المقياس العكسي غير معروفة، فإن توزيع غاما للمقياس العكسي يشكل توزيعًا مسبقًا مترافقًا. أما التوزيع المركب ، الناتج عن تكامل المقياس العكسي، فله حل مغلق يُعرف باسم توزيع غاما المركب . [ 23 ]

أما إذا كانت معلمة الشكل معروفة ولكن المتوسط ​​غير معروف، مع إعطاء التوزيع المسبق للمتوسط ​​بواسطة توزيع جاما آخر، فإن ذلك يؤدي إلى توزيع K.

الاستدلال الإحصائي

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

تقدير الاحتمال الأقصى

دالة الاحتمال لـ N مشاهدة مستقلة ومتطابقة التوزيع ( x1 , ..., xN ) هي

ل(α،θ)=أنا=1شمالو(xأنا؛α،θ){\displaystyle L(\alpha ,\theta )=\prod _{i=1}^{N}f(x_{i};\alpha ,\theta )}

ومنها نحسب دالة الاحتمال اللوغاريتمي

(α،θ)=(α-1)أنا=1شمالlnxأنا-أنا=1شمالxأناθ-شمالαlnθ-شمالlnΓ(α){\displaystyle \ell (\alpha ,\theta )=(\alpha -1)\sum _{i=1}^{N}\ln x_{i}-\sum _{i=1}^{N}{\frac {x_{i}}{\theta }}-N\alpha \ln \theta -N\ln \Gamma (\alpha )}

إيجاد القيمة القصوى بالنسبة إلى θ عن طريق اشتقاق الدالة ومساواتها بالصفر يُعطي مُقدِّر الاحتمال الأقصى لمعامل θ ، والذي يساوي متوسط ​​العينة.x¯{\displaystyle {\bar {x}}}مقسومًا على معامل الشكل α :

θ^=1αشمالأنا=1شمالxأنا=x¯α{\displaystyle {\hat {\theta }}={\frac {1}{\alpha N}}\sum _{i=1}^{N}x_{i}={\frac {\bar {x}}{\alpha }}}

وبتعويض هذا في دالة الاحتمال اللوغاريتمي نحصل على

(α)=(α-1)أنا=1شمالlnxأنا-شمالα-شمالαlnأناxأناαشمال-شمالlnΓ(α){\displaystyle \ell (\alpha )=(\alpha -1)\sum _{i=1}^{N}\ln x_{i}-N\alpha -N\alpha \ln {\frac {\sum _{i}x_{i}}{\alpha N}}-N\ln \Gamma (\alpha )}

نحتاج إلى عينتين على الأقل:شمال2{\displaystyle N\geq 2}لأن لـشمال=1{\displaystyle N=1}، الوظيفة(α){\displaystyle \ell (\alpha )}يزداد بلا حدود معα{\displaystyle \alpha \to \infty }. لα>0{\displaystyle \alpha >0}، ويمكن التحقق من ذلك(α){\displaystyle \ell (\alpha )}تكون الدالة مقعرة تمامًا ، وذلك باستخدام خصائص المتباينات لدالة بوليغاما . إيجاد القيمة العظمى بالنسبة إلى α عن طريق اشتقاق الدالة ومساواتها بالصفر يُعطي

lnα-ψ(α)=ln(1شمالأنا=1شمالxأنا)-1شمالأنا=1شمالlnxأنا=lnx¯-lnx¯{\displaystyle {\begin{aligned}\ln \alpha -\psi (\alpha )&=\ln \left({\frac {1}{N}}\sum _{i=1}^{N}x_{i}\right)-{\frac {1}{N}}\sum _{i=1}^{N}\ln x_{i}\\[1ex]&=\ln {\bar {x}}-{\overline {\ln x}}\end{aligned}}}

حيث ψ هي دالة ديغاما وlnx¯{\displaystyle {\overline {\ln x}}}يمثل المتوسط ​​الحسابي لـ ln x . لا يوجد حل مغلق لـ α . الدالة ذات سلوك عددي جيد للغاية، لذا إذا رُغِبَ في إيجاد حل عددي، فيمكن إيجاده باستخدام، على سبيل المثال، طريقة نيوتن . يمكن إيجاد قيمة ابتدائية لـ k إما باستخدام طريقة العزوم ، أو باستخدام التقريب.

lnα-ψ(α)12α(1+16α+1){\displaystyle \ln \alpha -\psi (\alpha )\approx {\frac {1}{2\alpha }}\left(1+{\frac {1}{6\alpha +1}}\right)}

إذا سمحنا

s=ln(1شمالأنا=1شمالxأنا)-1شمالأنا=1شمالlnxأنا=lnx¯-lnx¯{\displaystyle {\begin{aligned}s&=\ln \left({\frac {1}{N}}\sum _{i=1}^{N}x_{i}\right)-{\frac {1}{N}}\sum _{i=1}^{N}\ln x_{i}\\[1ex]&=\ln {\bar {x}}-{\overline {\ln x}}\end{aligned}}}

إذن فإن قيمة α تساوي تقريبًا

ك3-s+(s-3)2+24s12s{\displaystyle k\approx {\frac {3-s+{\sqrt {\left(s-3\right)^{2}+24s}}}{12s}}}

وهو ما يقع ضمن نطاق 1.5% من القيمة الصحيحة. [ 24 ] الصيغة الصريحة لتحديث نيوتن-رافسون لهذا التخمين الأولي هي: [ 25 ]

αα-lnα-ψ(ك)-s1α-ψ(α).{\displaystyle \alpha \leftarrow \alpha -{\frac {\ln \alpha -\psi (k)-s}{{\frac {1}{\alpha }}-\psi \prime (\alpha )}}.}

عند تقدير الاحتمالية القصوى(α^،θ^){\displaystyle ({\hat {\alpha }},{\hat {\theta }})}، القيم المتوقعة لـ x وlnx{\displaystyle \ln x}أتفق مع المتوسطات التجريبية: α^θ^=x¯وψ(α^)+lnθ^=lnx¯.{\displaystyle {\begin{aligned}{\hat {\alpha }}{\hat {\theta }}&={\bar {x}}&&{\text{and}}&\psi ({\hat {\alpha }})+\ln {\hat {\theta }}&={\overline {\ln x}}.\end{aligned}}}

تحذير بشأن معلمات الشكل الصغيرة

للحصول على البيانات،(x1،...،xشمال){\displaystyle (x_{1},\ldots ,x_{N})}، أي يتم تمثيلها بتنسيق الفاصلة العائمة الذي يتحول إلى 0 للقيم الأصغر منε{\displaystyle \varepsilon }ستؤدي اللوغاريتمات اللازمة لتقدير الاحتمال الأقصى إلى فشل العملية في حال وجود أي تدفقات سفلية. إذا افترضنا أن البيانات تم توليدها بواسطة توزيع جاما ذي دالة توزيع تراكميF(x؛α،θ){\displaystyle F(x;\alpha ,\theta )}إذن، فإن احتمال وجود تدفق سفلي واحد على الأقل هو: برو(التدفق السفلي)=1-(1-F(ε؛α،θ))شمال{\displaystyle \Pr({\text{underflow}})=1-(1-F(\varepsilon ;\alpha ,\theta ))^{N}} سيقترب هذا الاحتمال من 1 عندما تكون قيمة α صغيرة وقيمة N كبيرة . على سبيل المثال، عندα=10-2{\displaystyle \alpha =10^{-2}}، شمال=104{\displaystyle N=10^{4}}وε=2.25×10-308{\displaystyle \varepsilon =2.25\times 10^{-308}}،برو(التدفق السفلي)0.9998{\displaystyle \Pr({\text{underflow}})\approx 0.9998}يتمثل الحل البديل في استخدام البيانات بتنسيق لوغاريتمي.

لاختبار تطبيق مُقدِّر الاحتمال الأقصى الذي يأخذ بيانات لوغاريتمية كمدخلات، من المفيد أن تكون قادرًا على توليد لوغاريتمات غير متدفقة سفليًا لمتغيرات جاما العشوائية، عندماα<1{\displaystyle \alpha <1}بعد تطبيق ذلك في [26] scipy.stats.loggamma، يمكن القيام بذلك على النحو التالي: [ 26 ] مثالYجاما(α+1،θ){\displaystyle Y\sim {\text{Gamma}}(\alpha +1,\theta )}ويوزي مُوحد{\displaystyle U\sim {\text{Uniform}}}بشكل مستقل. ثم تكون العينة اللوغاريتمية المطلوبة هيZ=ln(Y)+ln(يو)/α{\displaystyle Z=\ln(Y)+\ln(U)/\alpha }، لهذا السببخبرة(Z)جاما(ك،θ){\displaystyle \exp(Z)\sim {\text{Gamma}}(k,\theta )}.

المقدرات ذات الصيغة المغلقة

توجد مقدرات متسقة مغلقة الشكل لـ α و θ مشتقة من احتمالية توزيع جاما المعمم . [ 27 ]

التقدير للشكل α هو

α^=شمالأنا=1شمالxأناشمالأنا=1شمالxأناlnxأنا-أنا=1شمالxأناأنا=1شمالlnxأنا{\displaystyle {\hat {\alpha }}={\frac {N\sum \limits _{i=1}^{N}x_{i}}{N\sum \limits _{i=1}^{N}x_{i}\ln x_{i}-\sum \limits _{i=1}^{N}x_{i}\sum \limits _{i=1}^{N}\ln x_{i}}}}

والتقدير للمقياس θ هو

θ^=1شمال2(شمالأنا=1شمالxأناlnxأنا-أنا=1شمالxأناأنا=1شمالlnxأنا){\displaystyle {\hat {\theta }}={\frac {1}{N^{2}}}\left(N\sum _{i=1}^{N}x_{i}\ln x_{i}-\sum _{i=1}^{N}x_{i}\sum _{i=1}^{N}\ln x_{i}\right)}

باستخدام المتوسط ​​الحسابي لـ x ، والمتوسط ​​الحسابي لـ ln x ، والمتوسط ​​الحسابي لحاصل ضرب x ·ln تتبسط التعبيرات إلى:

α^=x¯θ^{\displaystyle {\hat {\alpha }}={\frac {\bar {x}}{\hat {\theta }}}}θ^=xlnx¯-x¯lnx¯.{\displaystyle {\hat {\theta }}={\overline {x\ln x}}-{\bar {x}}{\overline {\ln x}}.}

إذا تم استخدام معلمات المعدل، فإن تقديرβ^=1/θ^{\displaystyle {\hat {\beta }}=1/{\hat {\theta }}}.

لا تُعدّ هذه المُقدِّرات مُقدِّرات احتمالية قصوى بالمعنى الدقيق، بل تُعرف باسم مُقدِّرات العزوم اللوغاريتمية المختلطة. ومع ذلك، فهي تتمتع بكفاءة مماثلة لمُقدِّرات الاحتمالية القصوى.

على الرغم من اتساق هذه التقديرات، إلا أنها تنطوي على انحياز طفيف. ويُعدّ أحد المتغيرات المصححة للانحياز لتقدير المقياس θ هو

θ~=شمالشمال-1θ^{\displaystyle {\tilde {\theta }}={\frac {N}{N-1}}{\hat {\theta }}}

يتم إعطاء تصحيح الانحياز لمعامل الشكل α على النحو التالي [ 28 ]

α~=α^-1شمال(3α^-23(α^1+α^)-45α^(1+α^)2){\displaystyle {\tilde {\alpha }}={\hat {\alpha }}-{\frac {1}{N}}\left(3{\hat {\alpha }}-{\frac {2}{3}}\left({\frac {\hat {\alpha }}{1+{\hat {\alpha }}}}\right)-{\frac {4}{5}}{\frac {\hat {\alpha }}{(1+{\hat {\alpha }})^{2}}}\right)}

متوسط ​​الخطأ التربيعي الأدنى البايزي

بمعرفة قيمة α وعدم معرفة قيمة θ ، تكون دالة الكثافة الاحتمالية اللاحقة لـ θ (باستخدام التوزيع الاحتمالي المسبق القياسي غير المتغير للمقياس لـ θ ) هي

برو(θ|α،x1،...،xشمال)1θأنا=1شمالو(xأنا؛α،θ){\displaystyle \Pr(\theta \mid \alpha ,x_{1},\dots ,x_{N})\propto {\frac {1}{\theta }}\prod _{i=1}^{N}f(x_{i};\alpha ,\theta )}

يدل على

yأنا=1شمالxأنا،برو(θ|α،x1،...،xشمال)=ج(xأنا)θ-شمالα-1هـ-y/θ{\displaystyle y\equiv \sum _{i=1}^{N}x_{i},\qquad \Pr(\theta \mid \alpha ,x_{1},\dots ,x_{N})=C(x_{i})\theta ^{-N\alpha -1}e^{-y/\theta }}

حيث لا يعتمد ثابت التكامل ( C ) على θ . يُظهر شكل دالة الكثافة الاحتمالية اللاحقة أن 1/ θ يتبع توزيع غاما بمعامل شكل + 2 ومعامل معدل y . يمكن إجراء التكامل بالنسبة إلى θ باستخدام تغيير المتغيرات لإيجاد ثابت التكامل.

0θ-شمالα-1+مهـ-y/θدθ=0xشمالα-1-مهـ-xyدx=y-(شمالα-م)Γ(شمالα-م){\displaystyle {\begin{aligned}\int _{0}^{\infty }\theta ^{-N\alpha -1+m}e^{-y/\theta }\,d\theta &=\int _{0}^{\infty }x^{N\alpha -1-m}e^{-xy}\,dx\\&=y^{-(N\alpha -m)}\Gamma (N\alpha -m)\!\end{aligned}}}

يمكن حساب العزوم عن طريق أخذ النسبة ( م / م = 0 )

هـ[xم]=Γ(شمالα-م)Γ(شمالα)yم{\displaystyle \operatorname {E} [x^{m}]={\frac {\Gamma (N\alpha -m)}{\Gamma (N\alpha )}}y^{m}}

وهذا يُظهر أن متوسط ​​± الانحراف المعياري لتقدير التوزيع الاحتمالي اللاحق لـ θ هو

yشمالα-1±y2(شمالα-1)2(شمالα-2).{\displaystyle {\frac {y}{N\alpha -1}}\pm {\sqrt {\frac {y^{2}}{\left(N\alpha -1\right)^{2}(N\alpha -2)}}}.}

الاستدلال البايزي

المرافق السابق

في الاستدلال البايزي ، يعتبر توزيع جاما هو التوزيع المسبق المترافق للعديد من توزيعات الاحتمالية: بواسون ، والأسي ، والطبيعي (مع متوسط ​​معروف)، وباريتو ، وجاما مع شكل معروف σ ، وجاما المعكوس مع معلمة شكل معروفة، وجومبيرتز مع معلمة مقياس معروفة.

التوزيع المسبق المترافق لتوزيع جاما هو: [ 29 ]

ص(α،θ|ص،q،ر،s)=1Zصα-1هـ-θ-1qΓ(α)رθαs،{\displaystyle p(\alpha ,\theta \mid p,q,r,s)={\frac {1}{Z}}{\frac {p^{\alpha -1}e^{-\theta ^{-1}q}}{\Gamma (\alpha )^{r}\theta ^{\alpha s}}},}

حيث Z هو ثابت التوحيد الذي لا يوجد له حل مغلق. يمكن إيجاد التوزيع الاحتمالي اللاحق بتحديث المعاملات كما يلي:

ص=صأناxأنا،q=q+أناxأنا،ر=ر+ن،s=s+ن،{\displaystyle {\begin{aligned}p'&=p\prod \nolimits _{i}x_{i},\\q'&=q+\sum \nolimits _{i}x_{i},\\r'&=r+n,\\s'&=s+n,\end{aligned}}}

حيث n هو عدد المشاهدات، و x i هي المشاهدة رقم i من توزيع جاما.

التواجد والتطبيقات

لنفترض سلسلة من الأحداث، حيث يكون وقت الانتظار لكل حدث توزيعًا أُسّيًا بمعدل β . عندئذٍ، يكون وقت الانتظار لحدوث الحدث رقم n هو توزيع جاما ذو شكل عددي صحيح.α=ن{\displaystyle \alpha =n}يُتيح هذا البناء لتوزيع غاما نمذجة مجموعة واسعة من الظواهر التي تتطلب حدوث عدة أحداث فرعية، يستغرق كل منها وقتًا بتوزيع أسي، بالتتابع حتى يقع حدث رئيسي. [ 30 ] تشمل الأمثلة وقت انتظار أحداث انقسام الخلايا ، [ 31 ] وعدد الطفرات التعويضية لطفرة معينة، [ 32 ] ووقت الانتظار حتى يصبح إصلاح النظام الهيدروليكي ضروريًا، [ 33 ] وما إلى ذلك.

في الفيزياء الحيوية، يكون زمن التوقف بين خطوات محرك جزيئي مثل إنزيم ATP synthase أسيًا تقريبًا عند تركيز ثابت من ATP، مما يدل على أن كل خطوة من خطوات المحرك تتطلب تحللًا واحدًا لجزيء ATP. إذا كان هناك n حدثًا لتحلل ATP، فسيكون التوزيع غاما من الدرجة n. [ 34 ]

استُخدم توزيع جاما لنمذجة حجم مطالبات التأمين [ 35 ] وكميات الأمطار. [ 36 ] وهذا يعني أن إجمالي مطالبات التأمين وكمية الأمطار المتراكمة في خزان يتم نمذجتها بواسطة عملية جاما - تمامًا كما يُولّد التوزيع الأسي عملية بواسون .

يُستخدم توزيع جاما أيضًا لنمذجة الأخطاء في نماذج انحدار بواسون متعددة المستويات لأن مزيجًا من توزيعات بواسون مع معدلات موزعة وفقًا لتوزيع جاما له توزيع مغلق معروف يسمى التوزيع ذو الحدين السالب .

في الاتصالات اللاسلكية، يتم استخدام توزيع جاما لنمذجة تلاشي طاقة الإشارة متعدد المسارات ؛ انظر أيضًا توزيع رايلي وتوزيع رايسيان .

في علم الأورام ، غالباً ما يتبع التوزيع العمري لحالات الإصابة بالسرطان توزيع جاما، حيث تتنبأ معلمات الشكل والمقياس، على التوالي، بعدد الأحداث المحركة والفاصل الزمني بينها. [ 37 ] [ 38 ]

في علم الأعصاب ، يُستخدم توزيع جاما غالبًا لوصف توزيع الفترات الزمنية بين النبضات العصبية . [ 39 ] [ 40 ]

في التعبير الجيني البكتيري حيث يمكن أن يحدث إنتاج البروتين على شكل دفعات، غالبًا ما يتبع عدد نسخ بروتين معين توزيع جاما، حيث يمثل معامل الشكل ومعامل المقياس، على التوالي، متوسط ​​عدد الدفعات لكل دورة خلوية ومتوسط ​​عدد جزيئات البروتين المنتجة لكل دفعة. [ 41 ]

في علم الجينوم ، تم تطبيق توزيع جاما في خطوة تحديد الذروة (أي في التعرف على الإشارة) في تحليل بيانات ChIP-chip [ 42 ] و ChIP-seq [ 43 ] .

في الإحصاء البايزي، يُستخدم توزيع غاما على نطاق واسع كتوزيع احتمالي مسبق مترافق . وهو التوزيع الاحتمالي المسبق المترافق لدقة ( أي معكوس التباين) التوزيع الطبيعي . كما أنه التوزيع الاحتمالي المسبق المترافق للتوزيع الأسي .

في علم الوراثة العرقي ، يُعد توزيع غاما النهج الأكثر شيوعًا لنمذجة تباين معدل التطور بين المواقع [ 44 ] عند استخدام طرق الاحتمال الأقصى ، أو الطرق البايزية ، أو طرق مصفوفة المسافة لتقدير الأشجار التطورية. تُقدّر التحليلات التطورية العرقية التي تستخدم توزيع غاما لنمذجة تباين معدل التطور مُعاملًا واحدًا من البيانات، لأنها تُقيّد الدراسة بالتوزيعات التي يكون فيها α = β . يعني هذا التحديد أن متوسط ​​هذا التوزيع هو 1، والتباين هو 1/ α . تستخدم طرق الاحتمال الأقصى والطرق البايزية عادةً تقريبًا منفصلًا لتوزيع غاما المستمر. [ 45 ] [ 46 ]

توليد المتغيرات العشوائية

بالنظر إلى خاصية القياس المذكورة أعلاه، يكفي توليد متغيرات جاما مع θ = 1 ، حيث يمكننا لاحقًا التحويل إلى أي قيمة لـ β من خلال عملية قسمة بسيطة.

لنفترض أننا نرغب في توليد متغيرات عشوائية من توزيع غاما ( n + δ , 1) ، حيث n عدد صحيح غير سالب و 0 < δ < 1. باستخدام حقيقة أن توزيع غاما (1, 1) هو نفسه توزيع Exp(1) ، وبالنظر إلى طريقة توليد المتغيرات الأسية ، نستنتج أنه إذا كان U موزعًا توزيعًا منتظمًا على الفترة (0, 1]، فإن −ln U موزعًا وفقًا لتوزيع غاما (1, 1) (أي أخذ العينات بالتحويل العكسي ). الآن، باستخدام خاصية " الجمع α " لتوزيع غاما، نوسع هذه النتيجة:

-ك=1نlnيوكΓ(ن،1){\displaystyle -\sum _{k=1}^{n}\ln U_{k}\sim \Gamma (n,1)}

حيث أن جميع المتغيرات U k موزعة توزيعًا منتظمًا على الفترة (0، 1] ومستقلة . كل ما تبقى الآن هو توليد متغير موزع وفقًا لتوزيع جاما ( δ ، 1) حيث 0 < δ < 1 ، وتطبيق خاصية " الجمع α " مرة أخرى. هذا هو الجزء الأصعب.

ناقش ديفروي [47] بالتفصيل توليد متغيرات جاما عشوائيًا، [ 47 ] : 401-428، مشيرًا إلى أن أيًا منها ليس سريعًا بشكل موحد لجميع معلمات الشكل. بالنسبة للقيم الصغيرة لمعلمة الشكل، غالبًا ما تكون الخوارزميات غير صالحة. [ 47 ] : 406. بالنسبة للقيم العشوائية لمعلمة الشكل، يمكن تطبيق خوارزمية القبول والرفض المعدلة لأهرنز وديتر [ 48 ] ، خوارزمية GD (الشكل α ≥ 1 )، أو طريقة التحويل [ 49 ] عندما 0 < α < 1. انظر أيضًا خوارزمية تشنغ وفيست GKM 3 [ 50 ] أو طريقة مارساجليا للضغط. [ 51 ]

فيما يلي نسخة من طريقة أهرنز-ديتر للقبول والرفض : [ 48 ]

  1. قم بتوليد U و V و W كمتغيرات مستقلة ومتطابقة التوزيع (0، 1).
  2. لويوهـهـ+دلتا{\displaystyle U\leq {\frac {e}{e+\delta }}}ثمξ=V1/دلتا{\displaystyle \xi =V^{1/\delta }}وη=دبليوξدلتا-1{\displaystyle \eta =W\xi ^{\delta -1}}. خلاف ذلك،ξ=1-lnV{\displaystyle \xi =1-\ln V}وη=دبليوهـ-ξ{\displaystyle \eta =We^{-\xi }}.
  3. لوη>ξدلتا-1هـ-ξ{\displaystyle \eta >\xi ^{\delta -1}e^{-\xi }}ثم انتقل إلى الخطوة 1.
  4. ξ يتم توزيعها كـ Γ( δ , 1) .

خلاصة هذا هي θ(ξ-أنا=1αlnيوأنا)Γ(α،θ){\displaystyle \theta \left(\xi -\sum _{i=1}^{\lfloor \alpha \rfloor }\ln U_{i}\right)\sim \Gamma (\alpha ,\theta )} أينα{\displaystyle \scriptstyle \lfloor \alpha \rfloor }يتم توليد الجزء الصحيح من α ، ويتم توليد ξ عبر الخوارزمية أعلاه مع δ = { α } (الجزء الكسري من α ) و U k كلها مستقلة.

على الرغم من صحة المنهج المذكور أعلاه من الناحية الفنية، يشير ديفروي إلى أنه يعتمد خطيًا على قيمة α، وبالتالي فهو ليس خيارًا جيدًا عمومًا. وبدلًا من ذلك، يوصي باستخدام طرق تعتمد على الرفض أو على الجداول، حسب السياق. [ 47 ] : 401-428

على سبيل المثال، طريقة التحويل والرفض البسيطة لمارساجليا التي تعتمد على متغير طبيعي واحد X ومتغير موحد واحد U : [ 26 ]

  1. تعييند=أ-13{\displaystyle d=a-{\frac {1}{3}}}وج=19د{\displaystyle c={\frac {1}{\sqrt {9d}}}}.
  2. تعيينv=(1+جX)3{\displaystyle v=(1+cX)^{3}}.
  3. لوv>0{\displaystyle v>0}وlnيو<X22+د-دv+دlnv{\displaystyle \ln U<{\frac {X^{2}}{2}}+d-dv+d\ln v}يعوددv{\displaystyle dv}وإلا فارجع إلى الخطوة 2.

مع1أ=α{\displaystyle 1\leq a=\alpha }يُولّد هذا الأسلوب عددًا عشوائيًا يتبع توزيع جاما مع مرور الوقت، ويكون هذا العدد ثابتًا تقريبًا مع قيمة α . يعتمد معدل القبول على قيمة α ، حيث يبلغ 0.95 و0.98 و0.99 عندما تكون α  تساوي  1 و  2  و  4 على التوالي. أما عندما تكون α أقل من 1 ، فيمكن استخدامγα=γ1+αيو1/α{\displaystyle \gamma _{\alpha }=\gamma _{1+\alpha }U^{1/\alpha }}لزيادة قيمة k لتكون قابلة للاستخدام مع هذه الطريقة.

في برنامج Matlab، يمكن توليد الأرقام باستخدام الدالة gamrnd()التي تستخدم تمثيل α و θ .

مراجع

  1. "توزيع جاما" . بريتانيكا . 6 يونيو 2026.
  2. وايسشتاين، إريك و. "توزيع جاما" . mathworld.wolfram.com . مؤرشف من الأصل بتاريخ 28-05-2024 . تم الاطلاع عليه بتاريخ 09-10-2024 .
  3. "توزيع غاما | دالة غاما | الخصائص | دالة كثافة الاحتمال" . مؤرشف من الأصل بتاريخ 13-06-2024 . تم الاطلاع عليه بتاريخ 09-10-2024 .
  4. بارك، سونغ ي.؛ بيرا، أنيل ك. (يونيو 2009). "نموذج الانحدار الذاتي الشرطي ذو التباين غير المتجانس ذو الإنتروبيا القصوى". مجلة الاقتصاد القياسي . 150 (2): 219-230 . Bibcode : 2009JEcon.150..219P . doi : 10.1016/j.jeconom.2008.12.014 .
  5. هوغ، آر في ؛ كريغ، إيه تي (1978). مقدمة في الإحصاء الرياضي ( الطبعة الرابعة). نيويورك: ماكميلان. ص. ملاحظة 3.3.1. ISBN   0-02-355710-9.
  6. جيلمان، أندرو؛ كارلين، جون ب.؛ ستيرن، هال س.؛ دونسون، ديفيد ب.؛ فيهتاري، آكي؛ روبين، دونالد ب. (2013). تحليل البيانات البايزي . الصفحات 576-577 . doi : 10.1201/b16018 . ISBN  978-0-429-11307-9.
  7. جوبالان، بريم؛ هوفمان، جيك م.؛ بلي، ديفيد م. (2013). "التوصية القابلة للتوسع باستخدام تحليل بواسون". arXiv : 1311.1704 [ cs.IR ].
  8. 1 2 بابوليس، أثاناسيوس؛ بيلاي، إس. أونيكريشنا (2002). الاحتمالات، والمتغيرات العشوائية، والعمليات العشوائية . ماكجرو هيل. ISBN 978-0-07-112256-6.
  9. تشين، جيسن؛ روبين، هيرمان (أكتوبر 1986). "حدود الفرق بين الوسيط والمتوسط ​​لتوزيعات جاما وبواسون". رسائل الإحصاء والاحتمالات . 4 (6): 281-283 . doi : 10.1016/0167-7152(86)90044-1 .
  10. تشوي، كيه بي (مايو 1994). "حول وسائط توزيعات جاما ومعادلة رامانوجان". وقائع الجمعية الرياضية الأمريكية . 121 (1): 245-251 . doi : 10.1090/S0002-9939-1994-1195477-8 .
  11. 1 2 بيرغ، كريستيان؛ بيدرسن، هنريك ل. (2006). "تخمين تشين-روبين في سياق متصل" . أساليب وتطبيقات التحليل . 13 (1): 63-88 . doi : 10.4310/MAA.2006.v13.n1.a4 .
  12. 1 2 بيرغ، كريستيان؛ بيدرسن، هنريك ل. (2006). "تحدب الوسيط في توزيع غاما". arXiv : math/0609442 .
  13. غونت، روبرت إي.؛ ميركل، ميلان (يناير 2021). "حول حدود المنوال والوسيط للتوزيع الزائدي المعمم والتوزيعات ذات الصلة". مجلة التحليل الرياضي والتطبيقات . 493 (1) 124508. arXiv : 2002.01884 . doi : 10.1016/j.jmaa.2020.124508 .
  14. 1 2 3 ليون، ريتشارد ف. (13 مايو 2021). "حول الحدود والتقريبات المغلقة المحكمة لوسيط توزيع جاما" . PLOS One . 16 (5) e0251626. arXiv : 2011.04060 . Bibcode : 2021PLoSO..1651626L . doi : 10.1371/journal.pone.0251626 . PMC 8118309. PMID 33984053 .  
  15. 1 2 ليون، ريتشارد ف . (13 مايو 2021). "حدود دقيقة لوسيط توزيع جاما" . PLOS One . 18 (9) e0288601. doi : 10.1371/journal.pone.0288601 . PMC 10490949. PMID 37682854 .  
  16. ماثاي، أ.م. (ديسمبر 1982). "سعة تخزين سد مع مدخلات من نوع جاما". حوليات معهد الإحصاء الرياضي . 34 (3): 591-597 . doi : 10.1007/BF02481056 .
  17. موشوبولوس، بي جي (1985). "توزيع مجموع متغيرات جاما العشوائية المستقلة". حوليات معهد الإحصاء الرياضي . 37 (3): 541-544 . Bibcode : 1985AnISM..37..541M . doi : 10.1007/BF02481123 .
  18. بيني، دبليو دي "انحرافات KL للكثافات الطبيعية، وجاما، وديريشليت، وويشارت" .
  19. لوبريت، أنتوني؛ بورج، يوهانس (أغسطس 2026). "عرض نطاق الدوال ذات الشكل التوزيعي غاما عبر دالة لامبرت W" . رسائل الإحصاء والاحتمالات . 235 110707. doi : 10.1016/j.spl.2026.110707 . PMC 13246144. PMID 42267171 .  
  20. "LogGammaDistribution—Wolfram Language Documentation" .
  21. "ExpGammaDistribution—Wolfram Language Documentation" .
  22. "scipy.stats.loggamma — SciPy v1.8.0 Manual" . docs.scipy.org .
  23. دوبي، ساتيا د. (ديسمبر 1970). "توزيعات جاما وبيتا وF المركبة". ميتريكا . 16 (1): 27-31 . doi : 10.1007/BF02613934 .
  24. مينكا، توماس ب. (2002). "تقدير توزيع جاما" (PDF) .
  25. تشوي، إس سي؛ ويت، آر. (1969). "تقدير الاحتمال الأقصى لمعلمات توزيع جاما وانحيازها". تكنومتركس . 11 (4): 683-690 . doi : 10.1080/00401706.1969.10490731 .
  26. 1 2 مارساجليا، جي.؛ تسانغ، دبليو دبليو (2000). "طريقة بسيطة لتوليد متغيرات جاما". معاملات ACM في البرمجيات الرياضية . 26 (3): 363-372 . doi : 10.1145/358407.358414 .
  27. يي، تشي-شنغ؛ تشين، نان (3 أبريل 2017). "مُقدِّرات الصيغة المغلقة لتوزيع غاما المُستمدة من معادلات الاحتمال". الإحصائي الأمريكي . 71 (2): 177-181 . doi : 10.1080/00031305.2016.1209129 .
  28. لوزادا، فرانسيسكو؛ راموس، بيدرو ل.؛ راموس، إدواردو (3 أبريل 2019). "ملاحظة حول انحياز مُقدِّرات الصيغة المغلقة لتوزيع جاما المُستمدة من معادلات الاحتمال". الإحصائي الأمريكي . 73 (2): 195-199 . doi : 10.1080/00031305.2018.1513376 . hdl : 11449/184543 .
  29. فينك، دانيال (مايو 1995). مجموعة من الأسبقيات المترافقة (PDF) (تقرير).
  30. جيسيكا، شاينر، صموئيل م.، 1956- غوريفيتش (2001). "13. تحليل وقت الفشل" . تصميم وتحليل التجارب البيئية . مطبعة جامعة أكسفورد. ISBN 0-19-513187-8. OCLC 43694448 . {{cite book}}: صيانة CS1: أسماء متعددة: قائمة المؤلفين ( رابط ) صيانة CS1: أسماء رقمية: قائمة المؤلفين ( رابط )
  31. غولوبيف، أ. (مارس 2016). "تطبيقات وآثار توزيع غاما المعدل أُسّيًا كنموذج للتغيرات الزمنية المتعلقة بتكاثر الخلايا والتعبير الجيني". مجلة البيولوجيا النظرية . 393 : 203-217 . Bibcode : 2016JThBi.393..203G . doi : 10.1016/j.jtbi.2015.12.027 . PMID 26780652 . 
  32. بون، آرت؛ ديفيس، برادلي هـ؛ تشاو، لين (يوليو 2005). " جامع القسائم وطفرة الكابت" . علم الوراثة . 170 (3): 1323-1332 . doi : 10.1534/genetics.104.037259 . PMC 1451182. PMID 15879511 .  
  33. فينيارد، مايكل؛ أمواكو-جيامبا، كواسي؛ ميريديث، جاك ر. (يوليو 1999). "توزيعات معدل الفشل لأنظمة التصنيع المرنة: دراسة تجريبية". المجلة الأوروبية لبحوث العمليات . 116 (1): 139-155 . doi : 10.1016/s0377-2217(98)00096-4 .
  34. ريف، ماتياس؛ روك، رونالد س.؛ ميهتا، أميت د.؛ موسكر، مارك س.؛ تشيني، ريتشارد إي.؛ سبوديتش، جيمس أ. (15 أغسطس 2000). "حركية حركة الميوسين-V: نموذج جزيئي للاستمرارية" . وقائع الأكاديمية الوطنية للعلوم . 97 (17): 9482-9486 . Bibcode : 2000PNAS...97.9482R . doi : 10.1073 / pnas.97.17.9482 . PMC 16890. PMID 10944217 .  
  35. بولاند، فيليب ج. (2007). الأساليب الإحصائية والاحتمالية في العلوم الاكتوارية . ص 43. doi : 10.1201/9781584886969 . ISBN  978-0-429-14557-5.
  36. ويلكس، دانيال س. (1990). "تقدير الاحتمال الأقصى لتوزيع جاما باستخدام بيانات تحتوي على أصفار" . مجلة المناخ . 3 ( 12): 1495-1501 . Bibcode : 1990JCli....3.1495W . doi : 10.1175/1520-0442(1990)003 < 1495:MLEFTG > 2.0.CO ; 2. JSTOR 26196366 . 
  37. بيليكوف، أليكسي ف. (22 سبتمبر 2017). "يمكن التنبؤ بعدد الأحداث المسرطنة الرئيسية من خلال معدل الإصابة بالسرطان" . التقارير العلمية . 7 (1) 12170. Bibcode : 2017NatSR...712170B . doi : 10.1038/s41598-017-12448-7 . PMC 5610194. PMID 28939880 .  
  38. بيليكوف، أليكسي ف.؛ فياتكين، أليكسي؛ ليونوف، سيرجي ف. (6 أغسطس 2021). "توزيع إيرلانغ يُقارب التوزيع العمري لحالات الإصابة بسرطانات الطفولة والشباب" . PeerJ . 9 e11976. Bibcode : 2021PeerJ...911976B . doi : 10.7717/peerj.11976 . PMC 8351573. PMID 34434669 .  
  39. روبسون، جي جي؛ تروي، جي بي (ديسمبر 1987). "طبيعة التفريغ المستمر لخلايا العقدة الشبكية Q وX وY في القطط". مجلة الجمعية البصرية الأمريكية أ . 4 (12): 2301-2307 . Bibcode : 1987JOSAA...4.2301R . doi : 10.1364/josaa.4.002301 . PMID 3430217 . 
  40. رايت، إم سي إم؛ وينتر، آي إم؛ فورستر، جيه جيه؛ بليك، إس. (نوفمبر 2014). "تخضع استجابة النواة القوقعية البطنية لنبضات النغمات ذات التردد الأمثل لإحصائيات فترات ما بين النبضات المرتبة". بحوث السمع . 317 : 23-32 . doi : 10.1016/j.heares.2014.09.006 . PMID 25261771 . 
  41. فريدمان، نير؛ كاي، لونغ؛ شي، إكس. سوني (19 أكتوبر 2006). "ربط الديناميكيات العشوائية بتوزيع السكان: إطار تحليلي للتعبير الجيني". رسائل المراجعة الفيزيائية . 97 (16) 168302. رمز Bibcode : 2006PhRvL..97p8302F . doi : 10.1103/PhysRevLett.97.168302 . PMID 17155441 . 
  42. ريس، ديفيد جيه؛ فاسيوطي، مارك تي؛ باليغا، نيتين إس. (فبراير 2008). "فك تشفير ارتباط الحمض النووي على مستوى الجينوم باستخدام النماذج". المعلوماتية الحيوية . 24 (3): 396-403 . doi : 10.1093/bioinformatics/btm592 . PMID 18056063 . 
  43. ميندوزا-بارا، ماركو-أنطونيو؛ نوفيكا، مالغورزاتا؛ فان غول، ووتر؛ غرونماير، هينريش (2013). "توصيف أنماط ارتباط ChIP-seq من خلال فك تشفير شكل الذروة القائم على النموذج" . BMC Genomics . 14 (1): 834. doi : 10.1186/1471-2164-14-834 . PMC 4046686. PMID 24279297 .  
  44. يانغ، زيهينغ (سبتمبر 1996). "تباين معدل الطفرات بين المواقع وتأثيره على التحليلات التطورية". اتجاهات في علم البيئة والتطور . 11 (9): 367-372 . Bibcode : 1996TEcoE..11..367Y . doi : 10.1016/0169-5347(96)10041-0 . PMID 21237881 . 
  45. يانغ، زيهينغ (سبتمبر 1994). "تقدير العلاقات التطورية باستخدام طريقة الاحتمال الأقصى من تسلسلات الحمض النووي ذات المعدلات المتغيرة عبر المواقع: طرق تقريبية". مجلة التطور الجزيئي . 39 (3): 306-314 . Bibcode : 1994JMolE..39..306Y . doi : 10.1007/BF00160154 . PMID 7932792 . 
  46. فيلسنشتاين، جوزيف (أكتوبر 2001). "مراعاة تباين معدلات التطور بين المواقع عند استنتاج العلاقات التطورية". مجلة التطور الجزيئي . 53 ( 4-5 ): 447-455 . Bibcode : 2001JMolE..53..447F . doi : 10.1007/s002390010234 . PMID 11675604 . 
  47. 1 2 3 ديفروي، لوك (1986). توليد المتغيرات العشوائية غير المنتظمة . doi : 10.1007/978-1-4613-8643-8 . ISBN 978-1-4613-8645-2.
  48. 1 2 أهرنز، جيه إتش؛ ديتر، يو (يناير 1982). "توليد متغيرات جاما بتقنية رفض معدلة" . اتصالات رابطة آلات الحوسبة . 25 (1): 47-54 . doi : 10.1145/358315.358390 .انظر الخوارزمية GD، صفحة  53.
  49. أهرنز، جيه إتش؛ ديتر، يو. (سبتمبر 1974). "طرق حاسوبية لأخذ العينات من توزيعات جاما، وبيتا، وبواسون، والتوزيعات الثنائية". الحوسبة . 12 (3): 223-246 . Bibcode : 1974Compu..12..223A . doi : 10.1007/BF02293108 .
  50. تشنغ، آر سي إتش؛ فيست، جي إم (1979). "بعض مولدات متغيرات جاما البسيطة". مجلة الجمعية الإحصائية الملكية. السلسلة ج (الإحصاء التطبيقي) . 28 (3): 290-295 . doi : 10.2307/2347200 . JSTOR 2347200 . 
  51. مارساجليا، جورج (1977). "طريقة الضغط لتوليد متغيرات جاما". الحوسبة والرياضيات مع التطبيقات . 3 (4): 321-325 . doi : 10.1016/0898-1221(77)90089-X .