توزيع بيتا-ذو الحدين

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

يُعدّ التوزيع بيتا-ذو الحدين نسخة أحادية البُعد من التوزيع ديريشلي-متعدد الحدود، كما أن التوزيعين ذي الحدين وبيتا هما نسختان أحاديتان من التوزيعين متعدد الحدود وديريشلي على التوالي. وتُعرف الحالة الخاصة التي يكون فيها α و β عددين صحيحين بالتوزيع الهندسي الفائق السالب .

الدافع والاشتقاق

كتوزيع مركب

يُعد توزيع بيتا توزيعًا مترافقًا لتوزيع ذي الحدين . وتؤدي هذه الحقيقة إلى توزيع مركب قابل للمعالجة التحليلية، حيث يمكن للمرء أن يفكر فيص{\displaystyle p}لنفترض أن المعلمة في التوزيع ذي الحدين يتم سحبها عشوائيًا من توزيع بيتا. لنفترض أننا مهتمون بتوقع احتمال ظهور عدد محدد من الصور،x{\displaystyle x}فين{\displaystyle n}التجارب المستقبلية. هذا ما يقدمه

و(x|ن،α،β)=01بأنان(x|ن،ص)بهـتأ(ص|α،β)دص=(نx)1ب(α،β)01صx+α-1(1-ص)ن-x+β-1دص=(نx)ب(x+α،ن-x+β)ب(α،β).{\displaystyle {\begin{aligned}f(x\mid n,\alpha ,\beta )&=\int _{0}^{1}\mathrm {Bin} (x|n,p)\mathrm {Beta} (p\mid \alpha ,\beta )\,dp\\[6pt]&={n \choose x}{\frac {1}{\mathrm {B} (\alpha ,\beta )}}\int _{0}^{1}p^{x+\alpha -1}(1-p)^{n-x+\beta -1}\,dp\\[6pt]&={n \choose x}{\frac {\mathrm {B} (x+\alpha ,n-x+\beta )}{\mathrm {B} (\alpha ,\beta )}}.\end{aligned}}}

باستخدام خصائص دالة بيتا ، يمكن كتابة ذلك بطريقة بديلة

و(x|ن،α،β)=Γ(ن+1)Γ(x+α)Γ(ن-x+β)Γ(ن+α+β)Γ(x+1)Γ(ن-x+1)Γ(α+β)Γ(α)Γ(β){\displaystyle f(x\mid n,\alpha ,\beta )={\frac {\Gamma (n+1)\Gamma (x+\alpha )\Gamma (n-x+\beta )}{\Gamma (n+\alpha +\beta )\Gamma (x+1)\Gamma (n-x+1)}}{\frac {\Gamma (\alpha +\beta )}{\Gamma (\alpha) )\جاما (\بيتا )}}}

كنموذج جرة

يمكن أيضًا تفسير توزيع بيتا-ذو الحدين باستخدام نموذج الجرة للأعداد الصحيحة الموجبة α و β ، والمعروف بنموذج جرة بوليا . تخيل جرة تحتوي على α كرة حمراء و β كرة سوداء، حيث تُجرى عمليات سحب عشوائية. إذا سُحبت كرة حمراء، تُعاد كرتان حمراوان إلى الجرة. وبالمثل، إذا سُحبت كرة سوداء، تُعاد كرتان سوداوان إلى الجرة. إذا تكرر هذا n مرة، فإن احتمال رؤية x كرة حمراء يتبع توزيع بيتا-ذو الحدين بمعاملات n و α و β . 

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

اللحظات والخصائص

اللحظات الثلاث الأولى الخام هي

μ1=نαα+βμ2=نα[ن(1+α)+β](α+β)(1+α+β)μ3=نα[ن2(1+α)(2+α)+3ن(1+α)β+β(β-α)](α+β)(1+α+β)(2+α+β)\begin{aligned}\mu_{1}=\frac{n\alpha}{\alpha +\beta}\\\mu_{2}=\frac{n\alpha[n(1+\alpha)+\beta]}{(\alpha +\beta)(1+\alpha +\beta)}\\\mu_{3}=\frac{n\alpha[n^2(1+\alpha)(2+\alpha)+3n(1+\alpha)\beta +\beta(\beta -\alpha)]}{(\alpha +\beta)(1+\alpha +\beta)(2+\alpha +\beta)}}\end{aligned}}}

والتفرطح هو

β2=(α+β)2(1+α+β)نαβ(α+β+2)(α+β+3)(α+β+ن)[(α+β)(α+β-1+6ن)+3αβ(ن-2)+6ن2-3αβن(6-ن)α+β-18αβن2(α+β)2].\beta_2 = \frac{(\alpha +\beta)^2(1+\alpha +\beta)}{n\alpha\beta(\alpha +\beta +2)(\alpha +\beta +3)(\alpha +\beta +n)} \left[(\alpha +\beta)(\alpha +\beta -1+6n)+3\alpha\beta(n-2)+6n^2 - \frac{3\alpha\betan(6-n)}{\alpha +\beta}} - \frac{18\alpha\betan^2}{(\alpha +\beta)^2} \right].

تأجيرص=αα+β{\displaystyle p={\frac {\alpha }{\alpha +\beta }}\!}نلاحظ، على سبيل الإيحاء، أن المتوسط ​​الحسابي يمكن كتابته على النحو التالي

μ=نαα+β=نص{\displaystyle \mu ={\frac {n\alpha }{\alpha +\beta }}=np\!}

والتباين كـ

σ2=نαβ(α+β+ن)(α+β)2(α+β+1)=نص(1-ص)α+β+نα+β+1=نص(1-ص)[1+(ن-1)ρ]{\displaystyle \sigma ^{2}={\frac {n\alpha \beta (\alpha +\beta +n)}{(\alpha +\beta )^{2}(\alpha +\beta +1)}}=np(1-p){\frac {\alpha +\beta +n}{\alpha +\beta +1}}=np(1-p)[1+(n-1)\rho ]\!}

أينρ=1α+β+1{\displaystyle \rho ={\tfrac {1}{\alpha +\beta +1}}\!}المعلمةρ{\displaystyle \rho \;\!} يُعرف هذا بالارتباط "داخل الفئة" أو "داخل المجموعة". هذا الارتباط الإيجابي هو الذي يُؤدي إلى التشتت الزائد. لاحظ أنه عندمان=1{\displaystyle n=1}، لا تتوفر معلومات للتمييز بين التباين بيتا والتباين ذي الحدين، والنموذجان لهما تباينات متساوية.

العزوم العاملية

العزم العاملي من الرتبة r لمتغير عشوائي بيتا-ذو حدين X هو

هـ[(X)ر]=ن!(ن-ر)!ب(α+ر،β)ب(α،β)=(ن)رب(α+ر،β)ب(α،β){\displaystyle \operatorname {E} {\bigl [}(X)_{r}{\bigr ]}={\frac {n!}{(nr)!}}{\frac {B(\alpha +r,\beta )}{B(\alpha ,\beta )}}=(n)_{r}{\frac {B(\alpha +r,\beta )}{B(\alpha ,\beta )}}}.

تقديرات نقطية

طريقة العزوم

يمكن الحصول على تقديرات طريقة العزوم من خلال ملاحظة العزم الأول والثاني لتوزيع بيتا ذي الحدين ومساواتهما بعزوم العينة.م1{\displaystyle m_{1}}وم2{\displaystyle m_{2}}نجد

α^=نم1-م2ن(م2م1-م1-1)+م1β^=(ن-م1)(ن-م2م1)ن(م2م1-م1-1)+م1.\begin{aligned}{\widehat {\alpha }}&={\frac {nm_{1}-m_{2}}{n({\frac {m_{2}}{m_{1}}}-m_{1}-1)+m_{1}}}\\[5pt]{\widehat {\beta }}&={\frac {(n-m_{1})(n-{\frac {m_{2}}{m_{1}}})}{n({\frac {m_{2}}{m_{1}}}-m_{1}-1)+m_{1}}}.\end{aligned}}}

قد تكون هذه التقديرات سلبية بشكل غير منطقي، مما يدل على أن البيانات إما غير متشتتة أو متشتتة بشكل أقل من التوزيع ذي الحدين. في هذه الحالة، يُعد كل من التوزيع ذي الحدين والتوزيع الهندسي الفائق مرشحين بديلين.

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

على الرغم من أن تقديرات الاحتمال الأقصى المغلقة غير عملية، نظرًا لأن دالة كثافة الاحتمال تتكون من دوال مشتركة ( دالة جاما و/أو دوال بيتا)، إلا أنه يمكن إيجادها بسهولة عبر التحسين العددي المباشر. يمكن حساب تقديرات الاحتمال الأقصى من البيانات التجريبية باستخدام طرق عامة لتركيب توزيعات بوليا متعددة الحدود، والتي وُصفت في (مينكا 2003). تُسهّل حزمة VGAM في لغة البرمجة R ، من خلال الدالة vglm، عبر الاحتمال الأقصى، تركيب نماذج من نوع glm مع استجابات موزعة وفقًا لتوزيع بيتا-ذو الحدين. لا يُشترط أن يكون n ثابتًا في جميع المشاهدات.

مثال: تباين نسبة الجنس

تُظهر البيانات التالية عدد الأطفال الذكور بين أول 12 طفلاً في عائلة مكونة من 13 فرداً، وذلك في 6115 عائلة مأخوذة من سجلات المستشفيات في ساكسونيا خلال القرن التاسع عشر (سوكال ورولف، ص  59 من ليندسي). وقد تم تجاهل الطفل الثالث عشر للحد من تأثير توقف العائلات عن الإنجاب بشكل غير عشوائي عند بلوغ جنس معين.

الذكور0123456789101112
العائلات324104286670103313431112829478181457

أول لحظتين من لحظات العينة هما

م1=6.23م2=42.31ن=12{\displaystyle {\begin{aligned}m_{1}&=6.23\\m_{2}&=42.31\\n&=12\end{aligned}}}

وبالتالي فإن تقديرات طريقة العزوم هي

α^=34.1350β^=31.6085.{\displaystyle {\begin{aligned}{\widehat {\alpha }}&=34.1350\\{\widehat {\beta }}&=31.6085.\end{aligned}}}

يمكن إيجاد تقديرات الاحتمال الأقصى عدديًا

α^ملهـ=34.09558β^ملهـ=31.5715{\displaystyle {\begin{aligned}{\widehat {\alpha }}_{\mathrm {mle} }&=34.09558\\{\widehat {\beta }}_{\mathrm {mle} }&=31.5715\end{aligned}}}

وأقصى احتمال لوغاريتمي هو

سجلل=-12492.9{\displaystyle \log {\mathcal {L}}=-12492.9}

ومنها نستنتج معيار معلومات أكايكي (AIC).

أأناج=24989.74.{\displaystyle {\mathit {AIC}}=24989.74.}

قيمة معيار معلومات أكايكي (AIC) للنموذج الثنائي المنافس هي  25070.34،  مما يدل على أن نموذج بيتا-ثنائي يُقدم توافقًا أفضل مع البيانات، أي أن هناك دليلًا على التشتت الزائد. ويُقدم تريفرز وويلارد تبريرًا نظريًا لعدم التجانس في الميل الجنسي بين صغار الثدييات.

تتجلى الملاءمة الفائقة بشكل خاص بين الذيل

الذكور0123456789101112
العائلات المرصودة324104286670103313431112829478181457
التوزيع المتوقع (بيتا-ذو الحدين)2.322.6104.8310.9655.71036.21257.91182.1853.6461.9177.943.85.2
القيمة المتوقعة (التوزيع الثنائي p = 0.519215)0.912.171.8258.5628.11085.21367.31265.6854.2410.0132.826.12.3

دورها في الإحصاءات البايزية

يلعب توزيع بيتا-ذو الحدين دورًا بارزًا في التقدير البايزي لاحتمالية نجاح برنوليص{\displaystyle p}والتي نرغب في تقديرها بناءً على البيانات. لنفترضX={X1،X2،Xن1}{\displaystyle \mathbf {X} =\{X_{1},X_{2},\cdots X_{n_{1}}\}}لتكن عينة من متغيرات برنولي العشوائية المستقلة والمتطابقة التوزيعXأنابرنولي(ص){\displaystyle X_{i}\sim {\text{Bernoulli}}(p)}لنفترض أن معرفتنا بـص{\displaystyle p}- على الطريقة البايزية - غير مؤكد ويتم نمذجته بواسطة التوزيع المسبقصبيتا(α،β){\displaystyle p\sim {\text{Beta}}(\alpha ,\beta )}. لوY1=أنا=1ن1Xأنا{\displaystyle Y_{1}=\sum _{i=1}^{n_{1}}X_{i}}ثم من خلال التراكم ، التوزيع التنبؤي المسبق لـ

Y1بيتا بين(ن1،α،β){\displaystyle Y_{1}\sim {\text{BetaBin}}(n_{1},\alpha ,\beta )}.

بعد الملاحظةY1{\displaystyle Y_{1}}نلاحظ أن التوزيع الاحتمالي اللاحق لـص{\displaystyle p}

و(ص|X،α،β)(أنا=1ن1صxأنا(1-ص)1-xأنا)صα-1(1-ص)β-1=جصxأنا+α-1(1-ص)ن1-xأنا+β-1=جصy1+α-1(1-ص)ن1-y1+β-1{\displaystyle {\begin{aligned}f(p|\mathbf {X} ,\alpha ,\beta )&\propto \left(\prod _{i=1}^{n_{1}}p^{x_{i}}(1-p)^{1-x_{i}}\right)p^{\alpha -1}(1-p)^{\beta -1}\\&=Cp^{\sum x_{i}+\alpha -1}(1-p)^{n_{1}-\sum x_{i}+\beta -1}\\&=Cp^{y_{1}+\alpha -1}(1-p)^{n_{1}-y_{1}+\beta -1}\end{aligned}}}

أينج{\displaystyle C}هو ثابت تطبيع . ندرك التوزيع الاحتمالي اللاحق لـص{\displaystyle p}كـبهـتأ(y1+α،ن1-y1+β){\displaystyle \mathrm {Beta} (y_{1}+\alpha ,n_{1}-y_{1}+\beta )}.

وهكذا، من خلال التراكم مرة أخرى، نجد أن التوزيع التنبؤي اللاحق للمجموعY2{\displaystyle Y_{2}}من عينة مستقبلية بحجمن2{\displaystyle n_{2}}لبهـرنouللأنا(ص){\displaystyle \mathrm {Bernoulli} (p)}المتغيرات العشوائية هي

Y2بهـتأبأنان(ن2،y1+α،ن1-y1+β){\displaystyle Y_{2}\sim \mathrm {BetaBin} (n_{2},y_{1}+\alpha ,n_{1}-y_{1}+\beta )}.

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

لرسم متغير عشوائي بيتا-ذو الحدينXبهـتأبأنان(ن،α،β){\displaystyle X\sim \mathrm {BetaBin} (n,\alpha ,\beta )}ارسم ببساطةصبهـتأ(α،β){\displaystyle p\sim \mathrm {Beta} (\alpha ,\beta )}ثم ارسمXب(ن،ص){\displaystyle X\sim \mathrm {B} (n,p)}.

  • بهـتأبأنان(1،α،β)بهـرنouللأنا(ص){\displaystyle \mathrm {BetaBin} (1,\alpha ,\beta )\sim \mathrm {Bernoulli} (p)\,}أينص=αα+β{\displaystyle p={\frac {\alpha }{\alpha +\beta }}\,}.
  • بهـتأبأنان(ن،1،1)يو(0،ن){\displaystyle \mathrm {BetaBin} (n,1,1)\sim U(0,n)\,}أينيو(أ،ب){\displaystyle U(a,b)\,}هو التوزيع المنتظم المنفصل .
  • لوXبهـتأبأنان(ن،α،β){\displaystyle X\sim \mathrm {BetaBin} (n,\alpha ,\beta )\,}ثم(ن-X)بهـتأبأنان(ن،β،α){\displaystyle (n-X)\sim \mathrm {BetaBin} (n,\beta ,\alpha )\,}
  • ليمsبهـتأبأنان(ن،صs،(1-ص)s)ب(ن،ص){\displaystyle \lim _{s\rightarrow \infty }\mathrm {BetaBin} (n,ps,(1-p)s)\sim \mathrm {B} (n,p)\,}أينص=αα+β{\displaystyle p={\frac {\alpha }{\alpha +\beta }}\,}وs=α+β{\displaystyle s=\alpha +\beta \,}وب(ن،ص){\displaystyle \mathrm {B} (n,p)\,}هو التوزيع ذو الحدين .
  • ليمنبهـتأبأنان(ن،نλ،ن2)Poأناs(λ){\displaystyle \lim _{n\rightarrow \infty }\mathrm {BetaBin} (n,n\lambda ,n^{2})\sim \mathrm {Pois} (\lambda )\,}أينPoأناs(λ){\displaystyle \mathrm {Pois} (\lambda )\,}هو توزيع بواسون .
  • ليمنبهـتأبأنان(ن،1،نص(1-ص))جيهـoم(ص){\displaystyle \lim _{n\rightarrow \infty }\mathrm {BetaBin} (n,1,{\frac {np}{(1-p)}})\sim \mathrm {Geom} (p)\,}أينجيهـoم(ص){\displaystyle \mathrm {Geom} (p)\,}هو التوزيع الهندسي .
  • ليمنبهـتأبأنان(ن،ر،نص(1-ص))شمالب(ر،ص){\displaystyle \lim _{n\rightarrow \infty }\mathrm {BetaBin} (n,r,{\frac {np}{(1-p)}})\sim \mathrm {NB} (r,p)\,}أينشمالب(ر،ص){\displaystyle \mathrm {NB} (r,p)\,}هو التوزيع ذو الحدين السالب .

انظر أيضاً

مراجع