توزيع فيشر الهندسي الفائق غير المركزي

دالة الكتلة الاحتمالية لتوزيع فيشر الهندسي الفائق غير المركزي لقيم مختلفة لنسبة الاحتمالات ω. m1 = 80، m2 = 60، n = 100، ω = 0.01، ...، 1000
عالم الأحياء والإحصائي رونالد فيشر

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

يمكن توضيح التوزيع باستخدام نموذج الجرة التالي. لنفترض ، على سبيل المثال، أن جرة تحتوي على m⁻¹ كرة حمراء و m⁻² كرة بيضاء، ليصبح المجموع N = m⁻¹ + m⁻² كرة . وزن كل كرة حمراء هو ω⁻¹ ، ووزن كل كرة بيضاء هو ω⁻² . نسبة الاحتمالات هي ω = ω⁻¹ / ω⁻² . الآن ، نسحب الكرات عشوائيًا بحيث يكون احتمال سحب كرة معينة متناسبًا مع وزنها، ولكنه مستقل عما يحدث للكرات الأخرى. يتبع عدد الكرات المسحوبة من لون معين التوزيع ذي الحدين . إذا كان العدد الإجمالي n للكرات المسحوبة معروفًا، فإن التوزيع الشرطي لعدد الكرات الحمراء المسحوبة لقيمة n معينة هو توزيع فيشر الهندسي الفائق غير المركزي. لتوليد هذا التوزيع تجريبيًا، علينا تكرار التجربة حتى نحصل على n كرة.

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

كلا التوزيعين متساويان مع التوزيع الهندسي الفائق (المركزي) عندما تكون نسبة الاحتمالات 1.

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

أُطلق على توزيع فيشر الهندسي الفائق غير المركزي اسم التوزيع الهندسي الفائق الممتد (هاركنس، 1965)، ولا يزال بعض المؤلفين يستخدمون هذا الاسم حتى اليوم.

التوزيع أحادي المتغير

يتم عرض دالة الاحتمال والمتوسط ​​والتباين في الجدول المجاور.

هناك تعبير بديل للتوزيع يتضمن كلاً من عدد الكرات المأخوذة من كل لون وعدد الكرات التي لم يتم أخذها كمتغيرات عشوائية، وبذلك يصبح التعبير عن الاحتمالية متناظرًا.

قد يطول وقت حساب دالة الاحتمال عندما يحتوي مجموع P₀ على العديد من الحدود. ويمكن تقليل وقت الحساب بحساب حدود المجموع بشكل تكراري نسبةً إلى الحد الخاص بـ y = مع تجاهل الحدود المهملة في الأطراف (Liao and Rosen, 2001).

يمكن تقريب المتوسط ​​من خلال:

μب-ب2-4أج2أ{\displaystyle \mu \approx {\frac {b-{\sqrt {b^{2}-4ac}}}{2a}}\,}،

أينأ=ω-1{\displaystyle a=\omega -1}،ب=م1+ن-شمال-(م1+ن)ω{\displaystyle b=m_{1}+nN-(m_{1}+n)\omega }،ج=م1نω{\displaystyle c=m_{1}n\omega }.

يمكن تقريب التباين بواسطة:

σ2شمالشمال-1/(1μ+1م1-μ+1ن-μ+1μ+م2-ن){\displaystyle \sigma ^{2}\approx {\frac {N}{N-1}}{\bigg /}\left({\frac {1}{\mu }}+{\frac {1}{m_{1}-\mu }}+{\frac {1}{n-\mu }}+{\frac {1}{\mu +m_{2}-n}}\right)}.

قدّم كلٌّ من ليفين (1984، 1990)، وماكولاغ ونيلدر (1989)، ولياو (1992)، وإيسينغا وبيلزر (2011) تقريبات أفضل للمتوسط ​​والتباين. وتُقدّم طرق نقطة السرج لتقريب المتوسط ​​والتباين، التي اقترحها إيسينغا وبيلزر (2011)، نتائج دقيقة للغاية.

ملكيات

تنطبق علاقات التناظر التالية:

fnchypg(x;ن،م1،شمال،ω)=fnchypg(ن-x;ن،م2،شمال،1/ω).{\displaystyle \operatorname {fnchypg} (x;n,m_{1},N,\omega )=\operatorname {fnchypg} (nx;n,m_{2},N,1/\omega )\,.}
fnchypg(x;ن،م1،شمال،ω)=fnchypg(x;م1،ن،شمال،ω).{\displaystyle \operatorname {fnchypg} (x;n,m_{1},N,\omega )=\operatorname {fnchypg} (x;m_{1},n,N,\omega )\,.}
fnchypg(x;ن،م1،شمال،ω)=fnchypg(م1-x;شمال-ن،م1،شمال،1/ω).{\displaystyle \operatorname {fnchypg} (x;n,m_{1},N,\omega )=\operatorname {fnchypg} (m_{1}-x;Nn,m_{1},N,1/\omega )\,.}

علاقة التكرار:

fnchypg(x;ن،م1،شمال،ω)=fnchypg(x-1;ن،م1،شمال،ω)(م1-x+1)(ن-x+1)x(م2-ن+x)ω.{\displaystyle \operatorname {fnchypg} (x;n,m_{1},N,\omega )=\operatorname {fnchypg} (x-1;n,m_{1},N,\omega ){\frac {(m_{1}-x+1)(n-x+1)}{x(m_{2}-n+x)}}\omega \,.}

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

الاشتقاق

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

لنفترض أن X و Y متغيران عشوائيان موزعان توزيعًا ثنائيًا، يمثلان عدد المستجيبين في مجموعتين متناظرتين بحجم m X و m Y على التوالي.

Xسلة المهملات(مX،πX)،Yسلة المهملات(مY،πY){\displaystyle X\sim \operatorname {Bin} (m_{X},\pi _{X}),\quad Y\sim \operatorname {Bin} (m_{Y},\pi _{Y})\,}.

نسبة احتمالاتهم هي

ω=ωXωY=πX/(1-πX)πY/(1-πY){\displaystyle \omega ={\frac {\omega _{X}}{\omega _{Y}}}={\frac {\pi _{X}/(1-\pi _{X})}{\pi _{Y}/(1-\pi _{Y})}}}.

انتشار المستجيبينπأنا{\displaystyle \pi _{i}}يتم تحديدها بالكامل من حيث الاحتمالاتωأنا{\displaystyle \omega _{i}}،أنا{X،Y}{\displaystyle i\in \{X,Y\}}والتي تتوافق مع تحيز أخذ العينات في مخطط الجرة المذكور أعلاه، أي

πأنا=ωأنا1+ωأنا{\displaystyle \pi _{i}={\frac {\omega _{i}}{1+\omega _{i}}}}.

يمكن تلخيص التجربة وتحليلها من خلال جدول الطوارئ التالي.

      علاج
  مجموعة
المستجيبغير مستجيبالمجموع
Xx.م س
Yy.لي
المجموعن.شمال

في الجدول،ن=x+y{\displaystyle n=x+y}يمثل هذا العدد الإجمالي للمستجيبين في جميع المجموعات، بينما يمثل N العدد الإجمالي للمرضى الذين تم تجنيدهم في التجربة. تشير النقاط إلى عدد التكرارات المقابلة التي لا علاقة لها بالموضوع.

توزيع عينة المستجيبين في المجموعة X بناءً على نتيجة التجربة ومعدلات الانتشار، Pر(X=x|X+Y=ن،مX،مY،ωX،ωY){\displaystyle Pr(X=x\;|\;X+Y=n,m_{X},m_{Y},\omega _{X},\omega _{Y})}، هي دالة هندسية فائقة غير مركزية:

F(X،ω):=Pر(X=x|X+Y=ن،مX،مY،ωX،ωY)=Pر(X=x،X+Y=ن|مX،مY،ωX،ωY)Pر(X+Y=ن|مX،مY،ωX،ωY)=Pر(X=x|مX،ωX)Pر(Y=ن-x|مY،ωY،X=x)Pر(X+Y=ن|مX،مY،ωX،ωY)=(مXx)πXx(1-πX)مX-x(مYن-x)πYن-x(1-πY)مY-(ن-x)Pر(X+Y=ن|مX،مY،ωX،ωY)=(مXx)ωXx(1-πX)مX(مYن-x)ωYن-x(1-πY)مYPر(X+Y=ن|مX،مY،ωX،ωY)=(مXx)(مYن-x)ωx(1-πX)مXωYن(1-πY)مY(1-πX)مXωYن(1-πY)مYu=الأعلى(0،ن-مY)مين(مX،ن)(مXu)(مYن-u)ωu=(مXx)(مYن-x)ωxu=الأعلى(0،ن-مY)مين(مX،ن)(مXu)(مYن-u)ωu{\displaystyle {\begin{aligned}F(X,\omega ):&=Pr(X=x\;|\;X+Y=n,m_{X},m_{Y},\omega _{X},\omega _{Y})\\&={\frac {Pr(X=x,X+Y=n\;|\;m_{X},m_{Y},\omega _{X},\omega _{Y})}{Pr(X+Y=n\;|\;m_{X},m_{Y},\omega _{X},\omega _{Y})}}\\&={\frac {Pr(X=x\;|\;m_{X},\omega _{X})Pr(Y=n-x\;|\;m_{Y},\omega _{Y},X=x)}{Pr(X+Y=n\;|\;m_{X},m_{Y},\omega _{X},\omega _{Y})}}\\&={\frac {{\binom {m_{X}}{x}}\pi _{X}^{x}(1-\pi _{X})^{m_{X}-x}{\binom {m_{Y}}{n-x}}\pi _{Y}^{n-x}(1-\pi _{Y})^{m_{Y}-(n-x)}}{Pr(X+Y=n\;|\;m_{X},m_{Y},\omega _{X},\omega _{Y})}}\\&={\frac {{\binom {m_{X}}{x}}\omega _{X}^{x}(1-\pi _{X})^{m_{X}}{\binom {m_{Y}}{n-x}}\omega _{Y}^{n-x}(1-\pi _{Y})^{m_{Y}}}{Pr(X+Y=n\;|\;m_{X},m_{Y},\omega _{X},\omega _{Y})}}\\&={\frac {{\binom {m_{X}}{x}}{\binom {m_{Y}}{n-x}}\omega ^{x}(1-\pi _{X})^{m_{X}}\omega _{Y}^{n}(1-\pi _{Y})^{m_{Y}}}{(1-\pi _{X})^{m_{X}}\omega _{Y}^{n}(1-\pi _{Y})^{m_{Y}}\sum _{u=\max(0,n-m_{Y})}^{\min(m_{X},n)}{\binom {m_{X}}{u}}{\binom {m_{Y}}{n-u}}\omega ^{u}}}\\&={\frac {{\binom {m_{X}}{x}}{\binom {m_{Y}}{n-x}}\omega ^{x}}{\sum _{u=\max(0,n-m_{Y})}^{\min(m_{X},n)}{\binom {m_{X}}{u}}{\binom {m_{Y}}{n-u}}\omega ^{u}}}\end{aligned}}}

لاحظ أن المقام هو في الأساس نفس البسط، بعد جمعه على جميع أحداث فضاء العينة المشترك.(X،Y){\displaystyle (X,Y)}ولذلك يرى أنX+Y=ن{\displaystyle X+Y=n}يمكن إخراج الحدود المستقلة عن X من المجموع وحذفها من البسط.

التوزيع متعدد المتغيرات

يمكن توسيع التوزيع ليشمل أي عدد من ألوان الكرات (c) في الجرة. ويُستخدم التوزيع متعدد المتغيرات عندما يكون هناك أكثر من لونين.

تُعرض دالة الاحتمال وتقريب بسيط للمتوسط ​​على اليمين. وقد قدم ماكولاغ ونيلدر (1989) تقريبات أفضل للمتوسط ​​والتباين.

ملكيات

ترتيب الألوان عشوائي بحيث يمكن تبديل أي لون.

يمكن تغيير حجم الأوزان بشكل تعسفي:

mfnchypg(x;ن،م،ω)=mfnchypg(x;ن،م،رω){\displaystyle \operatorname {mfnchypg} (\mathbf {x} ;n,\mathbf {m} ,{\boldsymbol {\omega }})=\operatorname {mfnchypg} (\mathbf {x} ;n,\mathbf {m} ,r{\boldsymbol {\omega }})\,\,}للجميعرR+.{\displaystyle r\in \mathbb {R} _{+}.}

يمكن حذف الألوان ذات الرقم الصفري ( m i = 0) أو الوزن الصفري (ω i = 0) من المعادلات.

يمكن دمج الألوان ذات الوزن نفسه:

mfnchypg(x;ن،م،(ω1،...،ωج-1،ωج-1))=mfnchypg((x1،...،xج-1+xج);ن،(م1،...،مج-1+مج)،(ω1،...،ωج-1))هيبج(xج;xج-1+xج،مج،مج-1+مج){\displaystyle {\begin{aligned}&{}\operatorname {mfnchypg} \left(\mathbf {x} ;n,\mathbf {m} ,(\omega _{1},\ldots ,\omega _{c-1},\omega _{c-1})\right)\\&{}=\operatorname {mfnchypg} \left((x_{1},\ldots ,x_{c-1}+x_{c});n,(m_{1},\ldots ,m_{c-1}+m_{c}),(\omega _{1},\ldots ,\omega _{c-1})\right)\,\cdot \\&\qquad \operatorname {hypg} (x_{c};x_{c-1}+x_{c},m_{c},m_{c-1}+m_{c})\end{aligned}}}

أينهيبج(x;ن،م،شمال){\displaystyle \operatorname {hypg} (x;n,m,N)}هو احتمال التوزيع الهندسي الفائق (أحادي المتغير، المركزي).

التطبيقات

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

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

البرامج المتاحة

انظر أيضاً

مراجع

بريسلو، نبراسكا؛ داي، نبراسكا (1980)، الأساليب الإحصائية في أبحاث السرطان ، ليون: الوكالة الدولية لأبحاث السرطان.

إيسينغا، ر.؛ بيلزر، ب. (2011)، "تقريبات نقطة السرج لمتوسط ​​وتباين التوزيع الهندسي الفائق الممتد" (ملف PDF) ، مجلة الإحصاء الهولندية ، المجلد  65، العدد  1، الصفحات 22-31 ، doi : 10.1111/j.1467-9574.2010.00468.x .

فوج، أ. (2007)، نظرية الأعداد العشوائية.

فوج، أ. (2008)، "طرق أخذ العينات لتوزيعات والينيوس وفيشر غير المركزية فوق الهندسية"، الاتصالات في الإحصاء والمحاكاة والحساب ، المجلد  37، العدد  2، الصفحات 241-257 ، doi : 10.1080/03610910701790236 ، S2CID 14904723  .

جونسون، إن إل؛ كيمب، إيه دبليو ؛ كوتز، إس. (2005)، التوزيعات المنفصلة أحادية المتغير ، هوبوكين، نيو جيرسي: وايلي وأولاده.

ليفين، ب. (1984)، "تحسينات بسيطة على تقريب كورنفيلد لمتوسط ​​متغير عشوائي غير مركزي ذي توزيع هندسي فائق"، مجلة Biometrika ، المجلد  71، العدد  3، الصفحات 630-632 ، doi : 10.1093/biomet/71.3.630 .

ليفين، ب. (1990)، "تصحيح نقطة السرج في تحليل الاحتمالية اللوجستية الشرطية"، مجلة Biometrika ، المجلد  77، العدد  2، [مطبعة جامعة أكسفورد، مؤسسة Biometrika]، الصفحات 275-285 ، doi : 10.1093/biomet/77.2.275 ، JSTOR 2336805  .

لياو، ج. (1992)، "خوارزمية لحساب المتوسط ​​والتباين للتوزيع الهندسي الفائق غير المركزي"، القياسات الحيوية ، المجلد  48، العدد  3، [وايلي، الجمعية الدولية للقياسات الحيوية]، الصفحات 889-892 ، doi : 10.2307/2532354 ، JSTOR 2532354  .

لياو، جي جي؛ روزن، أو. (2001)، "خوارزميات سريعة ومستقرة للحساب وأخذ العينات من التوزيع الهندسي الفائق غير المركزي"، الإحصائي الأمريكي ، المجلد  55، العدد  4، الصفحات 366-369 ، doi : 10.1198/000313001753272547 ، S2CID 121279235  .

مكولاج، ب.؛ نيلدر، ج. أ. (1989)، النماذج الخطية المعممة، الطبعة الثانية، لندن: تشابمان وهول.