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

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

يُعرف هذا التوزيع أيضًا باسم توزيع ماركوف-بوليا العكسي وتوزيع وارينغ المعمم [ 1 ] أو يُختصر ببساطة إلى توزيع BNB . ويُطلق على شكل مُزاح من هذا التوزيع اسم توزيع بيتا-باسكال . [ 1 ]

إذا كانت معلمات توزيع بيتاα{\displaystyle \alpha }وβ{\displaystyle \beta }وإذا

X|صشمالب(ر،ص)،{\displaystyle X\mid p\sim \mathrm {NB} (r,p),}

أين

صب(α،β)،{\displaystyle p\sim {\textrm {B}}(\alpha ,\beta ),}

ثم التوزيع الهامشي لـX{\displaystyle X}(أي التوزيع التنبؤي اللاحق ) هو توزيع بيتا ذو الحدين السالب:

Xبشمالب(ر،α،β).{\displaystyle X\sim \mathrm {BNB} (r,\alpha ,\beta ).}

في ما سبق،شمالب(ر،ص){\displaystyle \mathrm {NB} (r,p)}هو التوزيع ذو الحدين السالب وب(α،β){\displaystyle {\textrm {B}}(\alpha ,\beta )}هو توزيع بيتا .

التعريف والاشتقاق

يدل علىوX|ص(ك|q)،وص(q|α،β){\displaystyle f_{X|p}(k|q),f_{p}(q|\alpha ,\beta )}وباستخدام كثافات التوزيعات ذات الحدين السالبة وتوزيعات بيتا على التوالي، نحصل على دالة الكتلة الاحتماليةو(ك|α،β،ر){\displaystyle f(k|\alpha ,\beta ,r)} توزيع BNB عن طريق التهميش:

و(ك|α،β،ر)=01وX|ص(ك|ر،q)وص(q|α،β)دq=01(ك+ر-1ك)(1-q)كqرqα-1(1-q)β-1ب(α،β)دq=1ب(α،β)(ك+ر-1ك)01qα+ر-1(1-q)β+ك-1دq{\displaystyle {\begin{aligned}f(k|\alpha ,\beta ,r)\;=&\;\int _{0}^{1}f_{X|p}(k|r,q)\cdot f_{p}(q|\alpha ,\beta )\mathrm {d} q\\=&\;\int _{0}^{1}{\binom {k+r-1}{k}}(1-q)^{k}q^{r}\cdot {\frac {q^{\alpha -1}(1-q)^{\beta -1}}{\mathrm {B} (\alpha ,\beta )}}\mathrm {d} q\\=&\;{\frac {1}{\mathrm {B} (\alpha ,\beta )}}{\binom {k+r-1}{k}}\int _{0}^{1}q^{\alpha +r-1}(1-q)^{\beta +k-1}\mathrm {d} q\end{aligned}}}

مع ملاحظة أن قيمة التكامل تساوي:

01qα+ر-1(1-q)β+ك-1دq=Γ(α+ر)Γ(β+ك)Γ(α+β+ك+ر){\displaystyle \int _{0}^{1}q^{\alpha +r-1}(1-q)^{\beta +k-1}\mathrm {d} q={\frac {\Gamma (\alpha +r)\Gamma (\beta +k)}{\Gamma (\alpha +\beta +k+r)}}}

يمكننا التوصل إلى الصيغ التالية من خلال عمليات حسابية بسيطة نسبياً.

لور{\displaystyle r}إذا كان عددًا صحيحًا، فيمكن كتابة دالة الكتلة الاحتمالية بدلالة دالة بيتا .

و(ك|α،β،ر)=(ر+ك-1ك)ب(α+ر،β+ك)ب(α،β){\displaystyle f(k|\alpha ,\beta ,r)={\binom {r+k-1}{k}}{\frac {\mathrm {B} (\alpha +r,\beta +k)}{\mathrm {B} (\alpha ,\beta )}}}.

وبشكل أعم، يمكن كتابة دالة الكتلة الاحتمالية

و(ك|α،β،ر)=Γ(ر+ك)ك!Γ(ر)ب(α+ر،β+ك)ب(α،β){\displaystyle f(k|\alpha ,\beta ,r)={\frac {\Gamma (r+k)}{k!\;\Gamma (r)}}{\frac {\mathrm {B} (\alpha +r,\beta +k)}{\mathrm {B} (\alpha ,\beta )}}}

أو

و(ك|α،β،ر)=ب(ر+ك،α+β)ب(ر،α)Γ(ك+β)ك!Γ(β){\displaystyle f(k|\alpha ,\beta ,r)={\frac {\mathrm {B} (r+k,\alpha +\beta )}{\mathrm {B} (r,\alpha )}}{\frac {\Gamma (k+\beta )}{k!\;\Gamma (\beta )}}}.

تم التعبير عن PMF باستخدام جاما

باستخدام خصائص دالة بيتا ، دالة الكتلة الاحتمالية ذات العدد الصحيحر{\displaystyle r}يمكن إعادة كتابتها على النحو التالي:

و(ك|α،β،ر)=(ر+ك-1ك)Γ(α+ر)Γ(β+ك)Γ(α+β)Γ(α+ر+β+ك)Γ(α)Γ(β){\displaystyle f(k|\alpha ,\beta ,r)={\binom {r+k-1}{k}}{\frac {\Gamma (\alpha +r)\Gamma (\beta +k)\Gamma (\alpha +\beta )}{\Gamma (\alpha +r+\beta +k)\Gamma (\alpha )\Gamma (\beta )}}}.

وبشكل أعم، يمكن كتابة دالة الكتلة الاحتمالية على النحو التالي:

و(ك|α،β،ر)=Γ(ر+ك)ك!Γ(ر)Γ(α+ر)Γ(β+ك)Γ(α+β)Γ(α+ر+β+ك)Γ(α)Γ(β){\displaystyle f(k|\alpha ,\beta ,r)={\frac {\Gamma (r+k)}{k!\;\Gamma (r)}}{\frac {\Gamma (\alpha +r)\Gamma (\beta +k)\Gamma (\alpha +\beta )}{\Gamma (\alpha +r+\beta +k)\Gamma (\alpha )\Gamma (\beta) )}}}.

يُعبَّر عن PMF برمز بوشامر الصاعد

غالبًا ما يتم عرض دالة الكتلة الاحتمالية أيضًا بدلالة رمز بوشامر للأعداد الصحيحةر{\displaystyle r}

و(ك|α،β،ر)=ر(ك)α(ر)β(ك)ك!(α+β)(ر+ك){\displaystyle f(k|\alpha ,\beta ,r)={\frac {r^{(k)}\alpha ^{(r)}\beta ^{(k)}}{k!(\alpha +\beta )^{(r+k)}}}}

ملكيات

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

يُعرَّف العزم العاملي من الرتبة k لمتغير عشوائي ذي توزيع ثنائي سالب من نوع بيتا X على النحو التالي:ك<α{\displaystyle k<\alpha }وفي هذه الحالة يكون مساوياً لـ

هـ[(X)ك]=Γ(ر+ك)Γ(ر)Γ(β+ك)Γ(β)Γ(α-ك)Γ(α).{\displaystyle \operatorname {E} {\bigl [}(X)_{k}{\bigr ]}={\frac {\Gamma (r+k)}{\Gamma (r)}}{\frac {\Gamma (\beta +k)}{\Gamma (\beta )}}{\frac {\Gamma (\alpha -k)}{\Gamma (\alpha )}}.}

غير قابل للتحديد

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

العلاقة بالتوزيعات الأخرى

يحتوي توزيع بيتا ذو الحدين السالب على توزيع بيتا الهندسي كحالة خاصة عندما يكون إمار=1{\displaystyle r=1}أوβ=1{\displaystyle \beta =1}وبالتالي، يمكنه تقريب التوزيع الهندسي بدقة عالية. كما أنه يُقارب توزيع ذي الحدين السالب بدقة عالية أيضًا بالنسبة للقيم الكبيرة.α{\displaystyle \alpha }وبالتالي، يمكنها تقريب توزيع بواسون بشكل جيد للغاية بالنسبة للقيم الكبيرة.α{\displaystyle \alpha }،β{\displaystyle \beta }ور{\displaystyle r}.

ذي الذيل الثقيل

باستخدام تقريب ستيرلينغ لدالة بيتا، يمكن إثبات ذلك بسهولة بالنسبة للقيم الكبيرةك{\displaystyle k}

و(ك|α،β،ر)Γ(α+ر)Γ(ر)ب(α،β)كر-1(β+ك)ر+α{\displaystyle f(k|\alpha ,\beta ,r)\sim {\frac {\Gamma (\alpha +r)}{\Gamma (r)\mathrm {B} (\alpha ,\beta )}}{\frac {k^{r-1}}{(\beta +k)^{r+\alpha }}}}

وهذا يعني أن توزيع بيتا ذي الحدين السالب ذو ذيل سميك وأن العزوم أقل من أو تساويα{\displaystyle \alpha }غير موجود.

التوزيع الهندسي بيتا

يُعد التوزيع الهندسي بيتا حالة خاصة مهمة من توزيع بيتا ذي الحدين السالب، ويحدث لـر=1{\displaystyle r=1}في هذه الحالة، تتبسط دالة الكتلة الاحتمالية إلى

و(ك|α،β)=ب(α+1،β+ك)ب(α،β){\displaystyle f(k|\alpha ,\beta )={\frac {\mathrm {B} (\alpha +1,\beta +k)}{\mathrm {B} (\alpha ,\beta )}}}.

يتم استخدام هذا التوزيع في بعض نماذج "اشتر حتى تموت" (BTYD).

علاوة على ذلك، عندماβ=1{\displaystyle \beta =1}يختزل التوزيع الهندسي بيتا إلى توزيع يول-سيمون . ومع ذلك، من الشائع تعريف توزيع يول-سيمون بدلالة نسخة مُزاحة من التوزيع الهندسي بيتا. على وجه الخصوص، إذاXبجي(α،1){\displaystyle X\sim BG(\alpha ,1)}ثمX+1YS(α){\displaystyle X+1\sim YS(\alpha )}.

بيتا ذات الحدين السلبيين كنموذج لجرة بوليا

في حالة وجود المعايير الثلاثةر،α{\displaystyle r,\alpha }وβ{\displaystyle \beta }إذا كانت أعدادًا صحيحة موجبة، فيمكن أيضًا تفسير توزيع بيتا ذي الحدين السالب بنموذج الجرة - أو بشكل أكثر تحديدًا نموذج جرة بوليا الأساسي . لنفترض جرة تحتوي في البداية علىα{\displaystyle \alpha }الكرات الحمراء (لون التوقف) وβ{\displaystyle \beta }كرات زرقاء. في كل خطوة من خطوات النموذج، تُسحب كرة عشوائيًا من الجرة وتُعاد إلى مكانها، بالإضافة إلى كرة أخرى من نفس اللون. تُكرر العملية مرارًا وتكرارًا حتىر{\displaystyle r}يتم سحب كرات حمراء اللون. المتغير العشوائيX{\displaystyle X}يتم توزيع عمليات سحب الكرات الزرقاء المرصودة وفقًا لـبشمالب(ر،α،β){\displaystyle \mathrm {BNB} (r,\alpha ,\beta )}لاحظ أنه في نهاية التجربة، تحتوي الجرة دائمًا على العدد الثابتر+α{\displaystyle r+\alpha }من الكرات الحمراء التي تحتوي على الرقم العشوائيX+β{\displaystyle X+\beta }كرات زرقاء.

بموجب خاصية عدم قابلية التحديد،X{\displaystyle X}يمكن توليدها بشكل مكافئ باستخدام الجرة التي تحتوي في البداية علىα{\displaystyle \alpha }الكرات الحمراء (لون التوقف) ور{\displaystyle r}الخصيتان الزرقاوان والتوقف عندماβ{\displaystyle \beta }لوحظ وجود كرات حمراء.

انظر أيضاً

ملحوظات

  1. 1 2 جونسون وآخرون (1993)

مراجع

  • جونسون، إن إل؛ كوتز، إس؛ كيمب، إيه دبليو (1993) التوزيعات المنفصلة أحادية المتغير ، الطبعة الثانية، وايلي، رقم ISBN 0-471-54897-9 (القسم 6.2.3)
  • كيمب ، سي دي؛ كيمب، إيه دبليو (1956) "التوزيعات الهندسية الفائقة المعممة"، مجلة الجمعية الإحصائية الملكية ، السلسلة ب، 18، 202-211
  • وانغ، تشاوليانغ (2011) "توزيع ذي الحدين السالب المختلط مع تطبيق"، مجلة التخطيط والاستدلال الإحصائي ، 141 (3)، 1153-1160 doi : 10.1016/j.jspi.2010.09.020