توزيع جاما
في نظرية الاحتمالات والإحصاء ، يُعد توزيع غاما عائلةً متعددة الاستخدامات من التوزيعات الاحتمالية المستمرة ذات معلَمين . [ 1 ] ويُعتبر كلٌّ من التوزيع الأسي ، وتوزيع إيرلانغ ، وتوزيع كاي تربيع حالاتٍ خاصة من توزيع غاما. [ 2 ] وهناك طريقتان متكافئتان شائعتا الاستخدام لتمثيل هذا التوزيع :
- مع معلمة الشكلومعامل المقياس θ
- مع معلمة الشكلومعامل معدل
في كل من هذه الأشكال، يكون كلا المعاملين أعدادًا حقيقية موجبة.
لهذا التوزيع تطبيقات مهمة في مجالات متنوعة، منها الاقتصاد القياسي ، والإحصاء البايزي ، واختبارات العمر. [ 3 ] في الاقتصاد القياسي، يشيع استخدام معلمات ( α ، θ ) لنمذجة أوقات الانتظار، مثل الوقت حتى الوفاة، حيث غالبًا ما يأخذ شكل توزيع إيرلانغ لقيم α الصحيحة . أما الإحصائيون البايزيون فيفضلون معلمات ( α ، β )، مستخدمين توزيع غاما كتوزيع مسبق مترافق لعدة معلمات ذات مقياس معكوس، مما يسهل التحليل في حسابات التوزيع اللاحق.
يُعد توزيع جاما توزيع الاحتمالية ذو الإنتروبيا القصوى (سواء بالنسبة إلى مقياس أساسي موحد أوالقياس الأساسي) لمتغير عشوائي X حيث E [ X ] = αθ = α / β ثابت وأكبر من الصفر، و E [ln X ] = ψ ( α ) + ln θ = ψ ( α ) − ln β ثابت ( ψ هي دالة digamma ). [ 4 ]
التعريفات
يبدو أن استخدام المعاملات α و θ أكثر شيوعًا في الاقتصاد القياسي والمجالات التطبيقية الأخرى، حيث يُستخدم توزيع غاما بكثرة لنمذجة أوقات الانتظار. على سبيل المثال، في اختبارات العمر ، يُعد وقت الانتظار حتى الوفاة متغيرًا عشوائيًا يُنمذج غالبًا باستخدام توزيع غاما. انظر هوغ وكريغ [ 5 ] للاطلاع على دوافع واضحة.
يُعدّ استخدام المعاملات α و β أكثر شيوعًا في الإحصاء البايزي ، [ 6 ] حيث يُستخدم توزيع غاما كتوزيع احتمالي مسبق مترافق لأنواع مختلفة من معاملات المقياس العكسي (المعدل)، مثل λ للتوزيع الأسي أو توزيع بواسون [ 7 ] ، أو حتى β لتوزيع غاما نفسه. ويُستخدم توزيع غاما العكسي، ذو الصلة الوثيقة ، كتوزيع احتمالي مسبق مترافق لمعاملات المقياس، مثل تباين التوزيع الطبيعي .
إذا كان α عددًا صحيحًا موجبًا ، فإن التوزيع يمثل توزيع إيرلانغ ؛ أي مجموع α من المتغيرات العشوائية المستقلة الموزعة أُسّيًا ، ولكل منها متوسط θ .
التوصيف باستخدام الشكل α والمعدل β
يمكن تمثيل توزيع غاما بدلالة معامل الشكل α ومعامل المقياس العكسي β = 1/ θ ، والذي يُسمى معامل المعدل . ويُرمز للمتغير العشوائي X الذي يتبع توزيع غاما بمعامل شكل α ومعدل β بالرمز التالي:
دالة كثافة الاحتمال المقابلة في معلمات الشكل والمعدل هي
أيندالة غاما . لجميع الأعداد الصحيحة الموجبة،.
دالة التوزيع التراكمي هي دالة جاما المنتظمة:
أينهي دالة غاما غير المكتملة الدنيا .
إذا كان α عددًا صحيحًا موجبًا (أي أن التوزيع هو توزيع إيرلانغ )، فإن دالة التوزيع التراكمي لها التوسع المتسلسل التالي: [ 8 ]
التوصيف باستخدام الشكل α والمقياس θ
يُرمز إلى المتغير العشوائي X الذي يتبع توزيع جاما بشكل α ومقياس θ بالرمز التالي:

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

بخلاف المنوال والمتوسط، اللذين لهما صيغ حسابية سهلة تعتمد على المعلمات، فإن الوسيط لا يملك معادلة مغلقة. الوسيط لهذا التوزيع هو القيمةبحيث
قام تشين وروبين بمعالجة دقيقة لمشكلة تحديد التوسع التقاربي وحدود الوسيط لتوزيع جاما، حيث أثبتا أن (لـ) أينهو المتوسط وهو الوسيط لـالتوزيع. [ 9 ] بالنسبة للقيم الأخرى لمعامل المقياس، يتناسب المتوسط مع، وسيتم قياس الحدود والتقريبات الوسيطة بشكل مماثل بواسطة θ .
وجد كي بي تشوي الحدود الخمسة الأولى في تقريب لوران التقاربي للوسيط من خلال مقارنة الوسيط بـ رامانوجاندالة . [ 10 ] وجد بيرج وبيدرسن المزيد من المصطلحات: [ 11 ]


تعتبر المجاميع الجزئية لهذه السلاسل تقريبات جيدة لقيم α العالية بما فيه الكفاية ؛ فهي غير مرسومة في الشكل، الذي يركز على منطقة α المنخفضة التي يتم تقريبها بشكل أقل جودة.
أثبت بيرغ وبيدرسن أيضًا العديد من خصائص الوسيط، موضحين أنه دالة محدبة لـ α ، [ 12 ] وأن السلوك التقاربي بالقرب منيكون(حيث γ هو ثابت أويلر-ماسكيروني )، وذلك لجميعالوسيط محدود بـ[ 11 ]
حد أعلى خطي أدق، لـتم تقديمها فقط في عام 2021 بواسطة غونت وميركل، [ 13 ] بالاعتماد على نتيجة بيرغ وبيدرسن التي مفادها أن ميليكون في كل مكان أقل من 1: ل(مع المساواة في)) والتي يمكن توسيعها لتشمل حدًا لجميعبأخذ القيمة القصوى مع الوتر الموضح في الشكل، حيث ثبت أن الوسيط محدب. [ 12 ]
تقريب للوسيط يكون دقيقًا تقاربيًا عند قيم α العالية ومعقولًا حتىأو أقل قليلاً، ويتبع ذلك من تحويل ويلسون-هيلفرتي : والتي تصبح سالبة لـ.
في عام 2021، اقترح ليون عدة تقريبات من الشكللقد افترض قيمًا لـ A و B التي يكون عندها هذا التقريب حدًا أعلى أو أدنى محكمًا تقاربيًا لجميع[ 14 ] على وجه الخصوص ، اقترح هذه الحدود المغلقة الشكل، والتي أثبتها في عام 2023: [ 15 ]
يمثل حدًا أدنى، محكمًا تقاربًا كمايمثل حدًا أعلى، محكمًا تقاربًا كما
أظهر ليون أيضًا (بشكل غير رسمي في عام 2021، وبشكل دقيق في عام 2023) حدين أدنى آخرين ليسا تعبيرين مغلقين ، بما في ذلك هذا الحد الذي يتضمن دالة غاما ، استنادًا إلى حل التعبير التكاملي باستبدال 1 بـ: (الاقتراب من المساواة كـ) وخط المماس عندحيث تم العثور على المشتقة: (مع المساواة في)) حيث Ei هو التكامل الأسي . [ 14 ] [ 15 ]
بالإضافة إلى ذلك، أظهر أن عمليات الاستيفاء بين الحدود يمكن أن توفر تقريبات ممتازة أو حدودًا أدق للوسيط، بما في ذلك تقريب دقيق عند(أينويبلغ الحد الأقصى للخطأ النسبي أقل من 0.6%. جميع التقريبات والحدود المُستكملة هي من الشكل التالي: أينهي دالة استيفاء تعمل بشكل رتيب من 0 عند قيم α المنخفضة إلى 1 عند قيم α العالية ، وهي تقارب دالة استيفاء مثالية أو دقيقة.: بالنسبة لأبسط دالة استيفاء تم النظر فيها، دالة كسرية من الدرجة الأولى أضيق حد أدنى لديه وأضيق حد أعلى لديه تم رسم الحدود المُستكمَلة (معظمها داخل المنطقة الصفراء) في الرسم البياني اللوغاريتمي الموضح. تتوفر حدود أدق باستخدام دوال استيفاء مختلفة، ولكن ليس عادةً بمعاملات مغلقة الشكل كهذه. [ 14 ]
الخلاصة
إذا كان للمتغير العشوائي Xᵢ توزيع جاما ( αᵢ , θ ) حيث i = 1, 2, ..., N (أي أن جميع التوزيعات لها نفس معامل المقياس θ )، فإن
بشرط أن تكون جميع المتغيرات X i مستقلة .
بالنسبة للحالات التي تكون فيها المتغيرات X i مستقلة ولكن لها معلمات قياس مختلفة، انظر ماثاي [ 16 ] أو موشوبولوس [ 17 ] .
يُظهر توزيع جاما قابلية قسمة لا نهائية .
التوسع
لو
ثم، لأي قيمة لـ c > 0 ،
باستخدام دوال توليد العزوم،
أو ما يعادل ذلك، إذا
(تحديد معلمات الشكل والمعدل)
في الواقع، نعلم أنه إذا كان X متغيرًا عشوائيًا أسيًا بمعدل β ، فإن cX هو متغير عشوائي أسي بمعدل β / c ؛ وينطبق الشيء نفسه على متغيرات جاما (ويمكن التحقق من ذلك باستخدام دالة توليد العزوم ، انظر، على سبيل المثال، هذه الملاحظات ، 10.4-(ii)): الضرب في ثابت موجب c يقسم المعدل (أو، بشكل مكافئ، يضرب المقياس).
عائلة الدوال الأسية
توزيع جاما هو عائلة أسية ذات معلمتين مع معلمات طبيعية α − 1 و −1/ θ (بشكل مكافئ، α − 1 و − β )، وإحصائيات طبيعية X و ln X.
إذا تم تثبيت معلمة الشكل α ، فإن عائلة التوزيعات ذات المعلمة الواحدة الناتجة هي عائلة أسية طبيعية .
التوقع والتباين اللوغاريتميان
يمكن للمرء أن يثبت ذلك
أو ما يعادل ذلك،
حيث ψ هي دالة ديغاما . وبالمثل،
أينهي دالة تريغاما .
يمكن اشتقاق ذلك باستخدام صيغة العائلة الأسية لدالة توليد العزوم للإحصائية الكافية ، لأن إحدى الإحصائيات الكافية لتوزيع جاما هي ln x .
إنتروبيا المعلومات
إنتروبيا المعلومات هي
في تمثيل المعاملات α و θ ، تُعطى إنتروبيا المعلومات بواسطة
تباعد كولباك-لايبير

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