توزيع بيرسون

رسم تخطيطي لنظام بيرسون، يوضح توزيعات الأنواع الأول والثالث والسادس والخامس والرابع من حيث β 1 (التواء التربيع) و β 2 (التفرطح التقليدي).

يُعد توزيع بيرسون عائلة من التوزيعات الاحتمالية المستمرة . وقد نُشر لأول مرة بواسطة كارل بيرسون في عام 1895، ثم قام بتوسيعه لاحقًا في عامي 1901 و1916 في سلسلة من المقالات حول الإحصاء الحيوي .

تاريخ

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

في بحثه الأصلي،  حدد بيرسون (1895، ص 360) أربعة أنواع من التوزيعات (مرقمة من الأول إلى الرابع) بالإضافة إلى التوزيع الطبيعي (الذي كان يُعرف في الأصل بالنوع الخامس). اعتمد التصنيف على ما إذا كانت التوزيعات مدعومة على فترة محدودة، أو على نصف خط، أو على خط الأعداد الحقيقية بأكمله ؛ وما إذا كانت منحرفة أو متناظرة بالضرورة. عالج بحث ثانٍ (بيرسون 1901) نقصين: فقد أعاد تعريف توزيع النوع الخامس (الذي كان يُعرف في الأصل بالتوزيع الطبيعي فقط ، ولكنه يُعرف الآن بتوزيع غاما المعكوس ) وقدم توزيع النوع السادس. يغطي البحثان الأولان معًا الأنواع الخمسة الرئيسية لنظام بيرسون (الأول، والثالث، والرابع، والخامس، والسادس). في بحث ثالث، قدم بيرسون (1916) حالات خاصة وأنواعًا فرعية إضافية (من السابع إلى الثاني عشر).

ابتكر رايند (1909، الصفحات  430-432) طريقةً بسيطةً لتصوّر فضاء معلمات نظام بيرسون، والتي اعتمدها بيرسون لاحقًا (1916، اللوحة 1 والصفحات  430 وما بعدها، 448 وما بعدها). تتميز أنواع بيرسون بكميتين، يُشار إليهما عادةً بـ β1 و β2 . الأولى هي مربع معامل الالتواء : β1 = γ1 / حيث γ1 هو معامل الالتواء، أو العزم المعياري الثالث . أما الثاني فهو التفرطح التقليدي ، أو العزم المعياري الرابع: β 2 = γ 2 + 3. (تُعرّف المعالجات الحديثة التفرطح γ 2 بدلالة التراكمات بدلاً من العزوم، بحيث يكون لدينا في التوزيع الطبيعي γ 2 = 0 و β 2 = 3. هنا نتبع السابقة التاريخية ونستخدم β 2. ) يوضح الرسم البياني نوع بيرسون الذي ينتمي إليه توزيع محدد (يُحدد بنقطة (β 1 ، β 2 )).

العديد من التوزيعات الملتوية أو غير المتوسطة التفرطح، المألوفة لدى الإحصائيين اليوم، كانت لا تزال مجهولة في أوائل تسعينيات القرن التاسع عشر. ما يُعرف الآن بتوزيع بيتا، استخدمه توماس بايز كتوزيع احتمالي لاحق لمعلمة توزيع برنولي في عمله عام ١٧٦٣ حول الاحتمال العكسي . اكتسب توزيع بيتا شهرةً واسعةً لانضمامه إلى نظام بيرسون، وكان يُعرف حتى أربعينيات القرن العشرين بتوزيع بيرسون من النوع الأول. [ 1 ] (يُعدّ توزيع بيرسون من النوع الثاني حالة خاصة من النوع الأول، ولكنه لم يعد يُشار إليه عادةً بشكل منفصل). نشأ توزيع غاما من أعمال بيرسون (بيرسون 1893، ص  331؛ بيرسون 1895، ص  357، 360، 373-376) وكان يُعرف باسم توزيع بيرسون من النوع الثالث، قبل أن يكتسب اسمه الحديث في ثلاثينيات وأربعينيات القرن العشرين. [ 2 ] قدّم بيرسون في بحثه المنشور عام 1895 توزيع النوع الرابع، الذي يتضمن توزيع t للطالب كحالة خاصة، وذلك قبل استخدام ويليام سيلي جوسيت له لاحقًا بعدة سنوات. وقدّم في بحثه المنشور عام 1901 توزيع غاما المعكوس (النوع الخامس) وتوزيع بيتا برايم (النوع السادس).

تعريف

يتم تعريف كثافة بيرسون p على أنها أي حل صالح للمعادلة التفاضلية (انظر بيرسون 1902، ص  277).

ص(x)ص(x)+أ(x)ب(x)=0(1){\displaystyle {\frac {p'(x)}{p(x)}}+{\frac {A(x)}{B(x)}}=0\qquad (1)}

أين:

  • A ( x ) = a1 + a2x هي كثيرة حدود من الدرجة 1 على الأكثر ،
  • B ( x ) = b 0 + b 1 x + b 2 x 2 حيث B ≠ 0 هي متعددة حدود غير صفرية من الدرجة 2 على الأكثر،
  • ( a2 , b2 )(0, 0) :
    • إذا كان (a2 , b2 ) = ( 0, 0) ولكن ( a1 , b1 ) ≠ (0, 0) ، اضرب A و B في x + u مع ub0 / b1 قيمة محدودة اختيارية .
    • إذا كان (a2 , b2 ) = ( 0, 0) و ( a1 , b1 ) = (0, 0) (وبالتالي A = 0 )، اضرب A و B في ( x + u )( x + v ) مع u و v محدودين بشكل عشوائي بحيث uv .
  • المجال I لـ p هو فترة مفتوحة ( x min , x max ) محدودة باللانهايات و/أو الأصفار الحقيقية لـ B ولكنها لا تحتوي على أي منها ، أي :
    • إذا لم يكن للمتغير B أصفار حقيقية، فإنأنا=R{\displaystyle I=\mathbb {R} }،
    • إذا كان للمتغير B جذر حقيقي واحد x 0 ، فإن I = (−∞, x 0 ) أو I = ( x 0 , +∞) ،
    • إذا كان للمتغير B صفران حقيقيان x < x + ، فإن I = (−∞, x ) أو I = ( x , x + ) أو I = ( x + , +∞) .

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

يُظهر التحليل البُعدي أن المعاملات a k و b k لها بُعد [AX k ]، حيث [A] بُعد عشوائي و [X] بُعد x ، ومن ثم التحول في مؤشرات a k .

إذا كان للتوزيع لحظات حتى الرتبة 4 على الأقل، بمتوسط ​​μ وانحراف معياري σ ، فإننا نحصل على (انظر بيرسون 1916، ص  437، وكارفر 1924، ص  103-104).

ص(x)ص(x)+1σج1+ج2x-μσد0+د1x-μσ+د2(x-μσ)2=0{\displaystyle {\frac {p'(x)}{p(x)}}+{\frac {1}{\sigma }}{\frac {c_{1}+c_{2}{\frac {x-\mu }{\sigma }}}{d_{0}+d_{1}{\frac {x-\mu }{\sigma }}+d_{2}\left({\frac {x-\mu }{\sigma }}\right)^{2}}}=0}

مع

{ج2=2(5β2-6β1-9)د0=4β2-3β1ج1=د1=γ1(β2+3)د2=2β2-3β1-6{\displaystyle {\begin{cases}{\begin{aligned}c_{2}&=2(5\beta _{2}-6\beta _{1}-9)\\d_{0}&=4\beta _{2}-3\beta _{1}\\c_{1}=d_{1}&=\gamma _{1}(\beta _{2}+3)\\d_{2}&=2\beta _{2}-3\beta _{1}-6\end{aligned}}\end{cases}}}

يتم تعريفها حتى ثابت ضربي غير صفري للتناسب. علاوة على ذلك، فإن البعد الطبيعي للمعاملات a k و b k هو [X 8− k ].

دليل

البرهان أدناه هو توسيع للبرهان الموجز في كارفر 1924، الصفحات 103-104. وبالتحديد، تم افتراض أن [ q ( x ) D ( x ) xr ] J = 0 دون أي تبرير.

[...] إذا كانت دالة التردد، عند ضربها في t n ، تتلاشى عند حدود التوزيع [...]


ليكن X هو التوزيع ذو الكثافة p الذي يحقق المعادلة (1) ، مع المتوسط ​​μ ، والانحراف المعياري σ ، والالتواء γ 1 ، والتفرطح المحدود β 2 ، والمجال I = ( x min ، x max ) وفقًا للشروط المذكورة أعلاه فيما يتعلق بأصفار B : إذا كانت x min (على التوالي x max ) محدودة، فإن B ( x min ) = 0 (على التوالي B ( x max ) = 0 ). ليكن Z = Xμ / σ هو معياره ، والذي بالتالي له كثافة q ( x ) = σ p ( σx + μ ) ، ومجال J = ( x minμ / σ ، x maxμ / σ ) ، وعزومه الخام لها نفس القيم العددية (ولكن ليس نفس الأبعاد ) مثل العزوم المعيارية لـ X : J x r q ( x ) dx = α r ، مع α 0 = 1 ، α 1 = 0 ، α 2 = 1 ، α 3 = γ 1 ، و α 4 = β 2 .

كثافة Z تحقق الشرط التالي :

q(x)q(x)=σ2ص(σx+μ)σص(σx+μ)=-σأ(σx+μ)ب(σx+μ){\displaystyle {\frac {q'(x)}{q(x)}}={\frac {\sigma ^{2}p'(\sigma x+\mu )}{\sigma p(\sigma x+\mu )}}=-\sigma {\frac {A(\sigma x+\mu )}{B(\sigma x+\mu )}}}

والتي لها نفس بنية القسمة الخطية على التربيعية كما في المعادلة (1) . وبالتالي يمكننا وضع

q(x)q(x)+ج(x)د(x)=0(1){\displaystyle {\frac {q'(x)}{q(x)}}+{\frac {C(x)}{D(x)}}=0\qquad (1')}

حيث C ( x ) = c1 + c2x و D ( x ) = d0 + d1x + d2x2 . وبما أن D ( x ) تتناسب مع B ( σx + μ ) ، فإن الحدود المحدودة للمجال J هي أصفار D.

على العكس من ذلك، عندما p ( x ) = σ −1 q ( xμ / σ ) ، نحصل على

ص(x)ص(x)=σ-2q(x-μσ)σ-1q(x-μσ)=-1σج(x-μσ)د(x-μσ)،{\displaystyle {\frac {p'(x)}{p(x)}}={\frac {\sigma ^{-2}q'\left({\frac {x-\mu }{\sigma }}\right)}{\sigma ^{-1}q\left({\frac {x-\mu }{\sigma }}\right)}}=-{\frac {1}{\sigma }}{\frac {C\left({\frac {x-\mu }{\sigma }}\right)}{D\left({\frac {x-\mu }{\sigma }}\right)}},}

لذلك

ص(x)ص(x)+1σج1+ج2x-μσد0+د1x-μσ+د2(x-μσ)2=0{\displaystyle {\frac {p'(x)}{p(x)}}+{\frac {1}{\sigma }}{\frac {c_{1}+c_{2}{\frac {x-\mu }{\sigma }}}{d_{0}+d_{1}{\frac {x-\mu }{\sigma }}+d_{2}\left({\frac {x-\mu }{\sigma }}\right)^{2}}}=0}

مع تحديد المعاملات c k و d k .

بفرض عدد طبيعي r ∈ ⟦0, 3⟧ ، نضرب كل حد من حدود المعادلة (1') في q ( x ) D ( x ) xr ، فنحصل على

q(x)د(x)xر+q(x)(ج1xر+ج2xر+1)=0.{\displaystyle q'(x)D(x)x^{r}+q(x)(c_{1}x^{r}+c_{2}x^{r+1})=0.}

بإجراء التكامل على الحدين في المعادلة بالنسبة إلى J ، وإجراء التكامل بالتجزئة للحد الأول ، نحصل على

0=ج0دx=[q(x)د(x)xر]ج-جq(x)[د(x)xر]دx+جq(x)(ج1xر+ج2xر+1)دx=[q(x)د(x)xر]ج+جq(x)[ج1xر+ج2xر+1-رد0xر-1-(ر+1)د1xر-(ر+2)د2xر+1]دx=[q(x)د(x)xر]ج+ج1αر+ج2αر+1-رد0αر-1-(ر+1)د1αر-(ر+2)د2αر+1.{\displaystyle {\begin{aligned}0&=\int _{J}0\,dx\\&={\Big [}q(x)D(x)x^{r}{\Big ]}_{J}-\int _{J}q(x)[D(x)x^{r}]'\,dx+\int _{J}q(x)(c_{1}x^{r}+c_{2}x^{r+1})\,dx\\&={\Big [}q(x)D(x)x^{r}{\Big ]}_{J}+\int _{J}q(x)[c_{1}x^{r}+c_{2}x^{r+1}-rd_{0}x^{r-1}-(r+1)d_{1}x^{r}-(r+2)d_{2}x^{r+1}]\,dx\\&={\Big [}q(x)D(x)x^{r}{\Big ]}_{J}+c_{1}\alpha _{r}+c_{2}\alpha _{r+1}-rd_{0}\alpha _{r-1}-(r+1)d_{1}\alpha _{r}-(r+2)d_{2}\alpha _{r+1}.\end{aligned}}}

ليكن t حدًا، محدودًا أو غير محدود، للمجال J :

  • إذا كانت قيمة t محدودة:
    • إن تقارب J q ( x ) dx = 1 يعني أن لدينا q ( x ) = o (( xt ) −1 ) عندما xt ، وبالتالي فإن lim xt q ( x ) ( xt ) = 0 ،
    • t هو جذر لكثير الحدود التربيعي لذا فإن D ( x ) / xt هو كثير حدود خطي يأخذ قيمة محدودة عند t ،
    • t r محدود،
    • وبالتالي فإن lim xt q ( x ) D ( x ) x r = 0 ؛
  • إذا كانت قيمة t لانهائية:
    • إن وجود التفرطح لـ X ، وبالتالي لـ Z ، يعني أن لدينا x η q ( x ) = o ( x −1 ) على الأقل η ≤ 4 عندما xt ، وبالتالي lim xt q ( x ) x n = 0 على الأقل n ≤ 5 ،
    • D ( x ) x r هي متعددة حدود من الدرجة r +2 ≤ 5 على الأكثر ،
    • وبالتالي فإن lim xt q ( x ) D ( x ) x r = 0 .

وبالتالي، لدينا دائمًا [ q ( x ) D ( x ) x r ] J = 0 .

بجعل r تأخذ جميع القيم في ⟦0, 3⟧ ، نحصل بالتالي على نظام من 4 معادلات خطية :

{0=ج1α0+ج2α1-0د0α-1-1د1α0-2د2α1=ج1-د10=ج1α1+ج2α2-1د0α0-2د1α1-3د2α2=ج2-د0-3د20=ج1α2+ج2α3-2د0α1-3د1α2-4د2α3=ج1+γ1ج2-3د1-4γ1د20=ج1α3+ج2α4-3د0α2-4د1α3-5د2α4=γ1ج1+β2ج2-3د0-4γ1د1-5β2د2.{\displaystyle {\begin{cases}{\begin{alignedat}{3}0&=c_{1}\alpha _{0}+c_{2}\alpha _{1}-0d_{0}\alpha _{-1}-1d_{1}\alpha _{0}-2d_{2}\alpha _{1}&&=c_{1}-d_{1}\\0&=c_{1}\alpha _{1}+c_{2}\alpha _{2}-1d_{0}\alpha _{0}-2d_{1}\alpha _{1}-3d_{2}\alpha _{2}&&=c_{2}-d_{0}-3d_{2}\\0&=c_{1}\alpha _{2}+c_{2}\alpha _{3}-2d_{0}\alpha _{1}-3d_{1}\alpha _{2}-4d_{2}\alpha _{3}&&=c_{1}+\gamma _{1}c_{2}-3d_{1}-4\gamma _{1}d_{2}\\0&=c_{1}\alpha _{3}+c_{2}\alpha _{4}-3d_{0}\alpha _{2}-4d_{1}\alpha _{3}-5d_{2}\alpha _{4}&&=\gamma _{1}c_{1}+\beta _{2}c_{2}-3d_{0}-4\gamma _{1}d_{1}-5\beta _{2}d_{2}.\end{alignedat}}\end{cases}}}

وبالتالي ، لدينا c1 = d1 و c2 = d0 + 3d2 ، مما ينتج عنه

{0=د1+γ1(د0+3د2)-3د1-4γ1د2=γ1د0-2د1-γ1د20=γ1د1+β2(د0+3د2)-3د0-4γ1د1-5β2د2=(β2-3)د0-3γ1د1-2β2د2.{\displaystyle {\begin{cases}{\begin{alignedat}{3}0&=d_{1}+\gamma _{1}(d_{0}+3d_{2})-3d_{1}-4\gamma _{1}d_{2}&&=\gamma _{1}d_{0}-2d_{1}-\gamma _{1}d_{2}\\0&=\gamma _{1}d_{1}+\beta _{2}(d_{0}+3d_{2})-3d_{0}-4\gamma _{1}d_{1}-5\beta _{2}d_{2}&&=(\beta _{2}-3)d_{0}-3\gamma _{1}d_{1}-2\beta _{2}d_{2}.\end{alignedat}}\end{cases}}}

إذن d 1 = γ 1 ( d 0d 2 ) 2 ، وبالتالي

0=2β2د0-6د0-3γ12(د0-د2)-4β2د2=(2β2-6-3γ12)د0+(3γ12-4β2)د2=(2β2-3β1-6)د0-(4β2-3β1)د2.{\displaystyle {\begin{aligned}0&=2\beta _{2}d_{0}-6d_{0}-3\gamma _{1}^{2}(d_{0}-d_{2})-4\beta _{2}d_{2}\\&=(2\beta _{2}-6-3\gamma _{1}^{2})d_{0}+(3\gamma _{1}^{2}-4\beta _{2})d_{2}\\&=(2\beta _{2}-3\beta _{1}-6)d_{0}-(4\beta _{2}-3\beta _{1})d_{2}.\end{aligned}}}

بما أن المعاملات c k و d k يمكن ضربها جميعًا بثابت غير صفري دون تغيير المعادلة (1') ، فيمكننا بعد ذلك تحديدها بحرية

{د0=4β2-3β1د2=2β2-3β1-6{\displaystyle {\begin{cases}{\begin{aligned}d_{0}&=4\beta _{2}-3\beta _{1}\\d_{2}&=2\beta _{2}-3\beta _{1}-6\end{aligned}}\end{cases}}}

مما أدى إلى

{ج2=د0+3د2=2(5β2-6β1-9)ج1=د1=γ1(د0-د2)/2=γ1(β2+3).{\displaystyle {\begin{cases}{\begin{alignedat}{3}c_{2}&=d_{0}+3d_{2}&&=2(5\beta _{2}-6\beta _{1}-9)\\c_{1}=d_{1}&=\gamma _{1}(d_{0}-d_{2})/2&&=\gamma _{1}(\beta _{2}+3).\end{alignedat}}\end{cases}}}

بتوسيع γ 1 = α 3 = μ 3 σ −3 و β 2 = α 4 = μ 4 σ −4 وباستخدام μ 2 = σ 2 ، نحصل على

{ج2=2(5μ4σ-4-6μ32σ-6-9)=2σ-6(5μ4μ2-6μ32-9μ23)د0=4μ4σ-4-3μ32σ-6=σ-6(μ4μ2-3μ32)ج1=د1=μ3σ-3(μ4σ-4+3)=σ-7μ3(μ4+3μ22)د2=2μ4σ-4-3μ32σ-6-6=σ-6(2μ4μ2-3μ32-6μ23){\displaystyle {\begin{cases}{\begin{alignedat}{3}c_{2}&=2(5\mu _{4}\sigma ^{-4}-6\mu _{3}^{2}\sigma ^{-6}-9)&&=2\sigma ^{-6}(5\mu _{4}\mu _{2}-6\mu _{3}^{2}-9\mu _{2}^{3})\\d_{0}&=4\mu _{4}\sigma ^{-4}-3\mu _{3}^{2}\sigma ^{-6}&&=\sigma ^{-6}(\mu _{4}\mu _{2}-3\mu _{3}^{2})\\c_{1}=d_{1}&=\mu _{3}\sigma ^{-3}(\mu _{4}\sigma ^{-4}+3)&&=\sigma ^{-7}\mu _{3}(\mu _{4}+3\mu _{2}^{2})\\d_{2}&=2\mu _{4}\sigma ^{-4}-3\mu _{3}^{2}\sigma ^{-6}-6&&=\sigma ^{-6}(2\mu _{4}\mu _{2}-3\mu _{3}^{2}-6\mu _{2}^{3})\end{alignedat}}\end{cases}}}

وبالتالي

أ1+أ2xب0+ب1x+ب2x2=1σج1+ج2x-μσد0+د1x-μσ+د2(x-μσ)2=σ-7μ3(μ4+3μ22)+2σ-7(5μ4μ2-6μ32-9μ23)(x-μ)σ-5(μ4μ2-3μ32)+σ-7μ3(μ4+3μ22)(x-μ)+σ-7(2μ4μ2-3μ32-6μ23)(x-μ)2=μ3(μ4+3μ22)+2(5μ4μ2-6μ32-9μ23)(x-μ)μ2(μ4μ2-3μ32)+μ3(μ4+3μ22)(x-μ)+(2μ4μ2-3μ32-6μ23)(x-μ)2.{\displaystyle {\begin{aligned}{\frac {a_{1}+a_{2}x}{b_{0}+b_{1}x+b_{2}x^{2}}}&={\frac {1}{\sigma }}{\frac {c_{1}+c_{2}{\frac {x-\mu }{\sigma }}}{d_{0}+d_{1}{\frac {x-\mu }{\sigma }}+d_{2}\left({\frac {x-\mu }{\sigma }}\right)^{2}}}\\&={\frac {\sigma ^{-7}\mu _{3}(\mu _{4}+3\mu _{2}^{2})+2\sigma ^{-7}(5\mu _{4}\mu _{2}-6\mu _{3}^{2}-9\mu _{2}^{3})(x-\mu )}{\sigma ^{-5}(\mu _{4}\mu _{2}-3\mu _{3}^{2})+\sigma ^{-7}\mu _{3}(\mu _{4}+3\mu _{2}^{2})(x-\mu )+\sigma ^{-7}(2\mu _{4}\mu _{2}-3\mu _{3}^{2}-6\mu _{2}^{3})(x-\mu )^{2}}}\\&={\frac {\mu _{3}(\mu _{4}+3\mu _{2}^{2})+2(5\mu _{4}\mu _{2}-6\mu _{3}^{2}-9\mu _{2}^{3})(x-\mu )}{\mu _{2}(\mu _{4}\mu _{2}-3\mu _{3}^{2})+\mu _{3}(\mu _{4}+3\mu _{2}^{2})(x-\mu )+(2\mu _{4}\mu _{2}-3\mu _{3}^{2}-6\mu _{2}^{3})(x-\mu )^{2}}}.\end{aligned}}}

بما أن μ k لها بُعد [X k ] (و μ لها بُعد [X])، نرى أن إزالة σ s التي تظهر في المقامات ينتج عنها معاملات a k و b k ذات بُعد [X 8− k ].

بوضع ξ 4 = β 2 + 3 و ξ 6 = 3 β 2 − 3 β 1 − 3 ، نحصل على تعبيرات بدون ثوابت جمعية كما يلي:

{ج2=2(2ξ6-ξ4)د0=ξ6+ξ4ج1=د1=γ1ξ4د2=ξ6-ξ4.{\displaystyle {\begin{cases}{\begin{aligned}c_{2}&=2(2\xi _{6}-\xi _{4})\\d_{0}&=\xi _{6}+\xi _{4}\\c_{1}=d_{1}&=\gamma _{1}\xi _{4}\\d_{2}&=\xi _{6}-\xi _{4}.\end{aligned}}\end{cases}}}

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

في المعادلة (1) ، تحدد المعلمتان a1 و a2 نقطة ثابتة إذا كانت a2 0 ، وبالتالي تحددان، في ظل بعض الشروط، نمط التوزيع ، لأن

ص(-أ1أ2)=0{\displaystyle p'\left(-{\frac {a_{1}}{a_{2}}}\right)=0}

ويستنتج ذلك مباشرة من المعادلة التفاضلية.

بما أننا نواجه معادلة تفاضلية خطية من الدرجة الأولى ذات معاملات متغيرة ، فإن حلها مباشر:

ص(x)خبرة(-أ2x+أ1ب2x2+ب1x+ب0دx).{\displaystyle p(x)\propto \exp \left(-\int {\frac {a_{2}x+a_{1}}{b_{2}x^{2}+b_{1}x+b_{0}}}\,dx\right).}

يتبسط التكامل في هذا الحل بشكل كبير عند النظر في حالات خاصة معينة للدالة المراد تكاملها.  وقد ميز بيرسون (1895، ص 367) حالتين رئيسيتين، تحددهما إشارة المميز ( وبالتالي عدد الجذور الحقيقية ) للدالة التربيعية.

ب(x)=ب2x2+ب1x+ب0.(2){\displaystyle B(x)=b_{2}x^{2}+b_{1}x+b_{0}.\qquad (2)}

أنواع معينة من التوزيع

الحالة 1، التمييز السلبي

توزيع بيرسون من النوع الرابع

إذا كان مميز الدالة التربيعية (2) سالبًا (ب12-4ب2ب0<0{\displaystyle b_{1}^{2}-4b_{2}b_{0}<0})، ليس لها جذور حقيقية. ثم عرّف

y=x+ب12ب2،α=4ب2ب0-ب122ب2.{\displaystyle {\begin{aligned}y&=x+{\frac {b_{1}}{2b_{2}}},\\[5pt]\alpha &={\frac {\sqrt {4b_{2}b_{0}-b_{1}^{2}}}{2b_{2}}}.\end{aligned}}}

لاحظ أن α عدد حقيقي محدد جيدًا و α ≠ 0 ، لأنه بحسب الفرضية4ب2ب0-ب12>0{\displaystyle 4b_{2}b_{0}-b_{1}^{2}>0}وبالتالي فإن 0. بتطبيق هذه التعويضات، تتحول الدالة التربيعية (2) إلى

ب(x)=ب2(y2+α2).{\displaystyle B(x)=b_{2}(y^{2}+\alpha ^{2}).}

إن غياب الجذور الحقيقية واضح من هذه الصيغة، لأن α 2 موجب بالضرورة.

نعبر الآن عن حل المعادلة التفاضلية (1) كدالة لـ y :

ص(y)خبرة(-1ب2y-ب12ب2+أy2+α2دy).{\displaystyle p(y)\propto \exp \left(-{\frac {1}{b_{2}}}\int {\frac {y-{\frac {b_{1}}{2b_{2}}}+a}{y^{2}+\alpha ^{2}}}\,dy\right).}

أطلق بيرسون (1895، ص  362) على هذه الحالة اسم " الحالة المثلثية"، لأن التكامل

y-2ب2أ-ب12ب2y2+α2دy=12ln(y2+α2)-2ب2أ-ب12ب2αدالة الظل العكسي(yα)+ج0{\displaystyle \int {\frac {y-{\frac {2b_{2}a-b_{1}}{2b_{2}}}}{y^{2}+\alpha ^{2}}}\,dy={\frac {1}{2}}\ln(y^{2}+\alpha ^{2})-{\frac {2b_{2}a-b_{1}}{2b_{2}\alpha }}\arctan \left({\frac {y}{\alpha }}\right)+C_{0}}

يتضمن ذلك دالة الظل العكسي المثلثي . ثم

ص(y)خبرة[-12ب2ln(1+y2α2)-lnαب2+2ب2أ-ب12ب22αدالة الظل العكسي(yα)+ج1].{\displaystyle p(y)\propto \exp \left[-{\frac {1}{2b_{2}}}\ln \left(1+{\frac {y^{2}}{\alpha ^{2}}}\right)-{\frac {\ln \alpha }{b_{2}}}+{\frac {2b_{2}a-b_{1}}{2b_{2}^{2}\alpha }}\arctan \left({\frac {y}{\alpha }}\right)+C_{1}\right].}

وأخيراً، دع

م=12ب2،ν=-2ب2أ-ب12ب22α.{\displaystyle {\begin{aligned}m&={\frac {1}{2b_{2}}},\\[5pt]\nu &=-{\frac {2b_{2}a-b_{1}}{2b_{2}^{2}\alpha }}.\end{aligned}}}

بتطبيق هذه الاستبدالات، نحصل على الدالة البارامترية:

ص(y)[1+y2α2]-مخبرة[-νدالة الظل العكسي(yα)].{\displaystyle p(y)\propto \left[1+{\frac {y^{2}}{\alpha ^{2}}}\right]^{-m}\exp \left[-\nu \arctan \left({\frac {y}{\alpha }}\right)\right].}

تتمتع هذه الكثافة غير المعيارية بدعم على كامل خط الأعداد الحقيقية . وهي تعتمد على معامل المقياس α > 0 ومعاملي الشكل m > 1/2 و ν . فُقد أحد المعاملات عندما اخترنا إيجاد حل المعادلة التفاضلية (1) كدالة لـ y بدلاً من x . لذلك، نعيد إدخال معامل رابع، وهو معامل الموقع λ . وبذلك نكون قد استنتجنا كثافة توزيع بيرسون من النوع الرابع . 

ص(x)=|Γ(م+ν2أنا)Γ(م)|2αب(م-12،12)[1+(x-λα)2]-مخبرة[-νدالة الظل العكسي(x-λα)].{\displaystyle p(x)={\frac {\left|{\frac {\operatorname {\Gamma } \left(m+{\frac {\nu }{2}}i\right)}{\Gamma (m)}}\right|^{2}}{\alpha \operatorname {\mathrm {B} } \left(m-{\frac {1}{2}},{\frac {1}{2}}\right)}}\left[1+\left({\frac {x-\lambda }{\alpha }}\right)^{2}\right]^{-m}\exp \left[-\nu \arctan \left({\frac {x-\lambda }{\alpha }}\right)\right].}

يتضمن ثابت التطبيع دالة غاما المركبة (Γ) ودالة بيتا (B). لاحظ أن معامل الموقع λ هنا ليس هو نفسه معامل الموقع الأصلي المُقدم في الصيغة العامة، ولكنه مرتبط به عبر 

λ=λoرأنازأنانأل+αν2(م-1).{\displaystyle \lambda =\lambda _{original}+{\frac {\alpha \nu }{2(m-1)}}.}

توزيع بيرسون من النوع السابع

رسم بياني لكثافات بيرسون من النوع السابع مع λ = 0، σ = 1، و: γ 2 = ∞ (أحمر)؛ γ 2 = 4 (أزرق)؛ و γ 2 = 0 (أسود)

يتحكم معامل الشكل ν لتوزيع بيرسون من النوع الرابع في درجة التواء التوزيع . إذا ثبّتنا قيمته عند الصفر، نحصل على عائلة متناظرة ذات ثلاثة معاملات. تُعرف هذه الحالة الخاصة بتوزيع بيرسون من النوع السابع (انظر بيرسون 1916، ص  450). كثافته هي

ص(x)=1αب(م-12،12)[1+(x-λα)2]-م،{\displaystyle p(x)={\frac {1}{\alpha \operatorname {\mathrm {B} } \left(m-{\frac {1}{2}},{\frac {1}{2}}\right)}}\left[1+\left({\frac {x-\lambda }{\alpha }}\right)^{2}\right]^{-m},}

حيث B هي دالة بيتا .

يتم الحصول على معلمات بديلة (وتخصيص طفيف) لتوزيع النوع السابع عن طريق وضع

α=σ2م-3،{\displaystyle \alpha =\sigma {\sqrt {2m-3}},}

وهذا يتطلب أن تكون قيمة m أكبر من 3/2. وهذا يُؤدي إلى فقدان طفيف للعمومية، ولكنه يضمن وجود تباين التوزيع ومساواته لـ σ² . الآن، يتحكم المعامل m فقط في تفرطح التوزيع. إذا اقتربت قيمة m من اللانهاية مع ثبات قيمتي λ و σ ، فإن التوزيع الطبيعي يظهر كحالة خاصة.

ليمم1σ2م-3ب(م-12،12)[1+(x-λσ2م-3)2]-م=1σ2Γ(12)ليممΓ(م)Γ(م-12)م-32ليمم[1+(x-λσ)22م-3]-م=1σ2π1خبرة[-12(x-λσ)2].{\displaystyle {\begin{aligned}&\lim _{m\to \infty }{\frac {1}{\sigma {\sqrt {2m-3}}\,\operatorname {\mathrm {B} } \left(m-{\frac {1}{2}},{\frac {1}{2}}\right)}}\left[1+\left({\frac {x-\lambda }{\sigma {\sqrt {2m-3}}}}\right)^{2}\right]^{-m}\\[5pt]={}&{\frac {1}{\sigma {\sqrt {2}}\,\operatorname {\Gamma } \left({\frac {1}{2}}\right)}}\cdot \lim _{m\to \infty }{\frac {\Gamma (m)}{\operatorname {\Gamma } \left(m-{\frac {1}{2}}\right){\sqrt {m-{\frac {3}{2}}}}}}\cdot \lim _{m\to \infty }\left[1+{\frac {\left({\frac {x-\lambda }{\sigma }}\right)^{2}}{2m-3}}\right]^{-m}\\[5pt]={}&{\frac {1}{\sigma {\sqrt {2\pi }}}}\cdot 1\cdot \exp \left[-{\frac {1}{2}}\left({\frac {x-\lambda }{\sigma }}\right)^{2}\right].\end{aligned}}}

هذه هي كثافة التوزيع الطبيعي بمتوسط ​​λ وانحراف معياري σ .

من الملائم اشتراط أن تكون قيمة m أكبر من 5/2 والسماح

م=52+3γ2.{\displaystyle m={\frac {5}{2}}+{\frac {3}{\gamma _{2}}}.}

هذا تخصص آخر، وهو يضمن وجود اللحظات الأربع الأولى للتوزيع. وبشكل أكثر تحديدًا، فإن توزيع بيرسون من النوع السابع، المُعَلم بدلالة (λ، σ، γ² ) ، له متوسط ​​λ ، وانحراف معياري σ ، والتواء يساوي صفرًا، وتفرطح زائد موجب γ² .

توزيع t للطالب

إن توزيع بيرسون من النوع السابع يعادل توزيع t للطالب غير المعياري مع المعلمات ν > 0، μ، σ 2 عن طريق تطبيق الاستبدالات التالية على معلماته الأصلية:

λ=μ،α=νσ2،م=ν+12،{\displaystyle {\begin{aligned}\lambda &=\mu ,\\[5pt]\alpha &={\sqrt {\nu \sigma ^{2}}},\\[5pt]m&={\frac {\nu +1}{2}},\end{aligned}}}

لاحظ أن القيد m > 1/2 مُحقق.

الكثافة الناتجة هي

ص(x|μ،σ2،ν)=1νσ2ب(ν2،12)(1+1ν(x-μ)2σ2)-ν+12،{\displaystyle p(x\mid \mu ,\sigma ^{2},\nu )={\frac {1}{{\sqrt {\nu \sigma ^{2}}}\,\operatorname {\mathrm {B} } \left({\frac {\nu }{2}},{\frac {1}{2}}\right)}}\left(1+{\frac {1}{\nu }}{\frac {(x-\mu )^{2}}{\sigma ^{2}}}\right)^{-{\frac {\nu +1}{2}}},}

والتي يمكن التعرف عليها بسهولة على أنها كثافة توزيع t للطالب .

هذا يعني أن توزيع بيرسون من النوع السابع يشمل توزيع t القياسي للطالب ، وكذلك توزيع كوشي القياسي . وبالتحديد، يظهر توزيع t القياسي للطالب كحالة فرعية عندما يكون μ = 0 و σ² = 1 ، وهو ما يعادل الاستبدالات التالية:

λ=0،α=ν،م=ν+12،{\displaystyle {\begin{aligned}\lambda &=0,\\[5pt]\alpha &={\sqrt {\nu }},\\[5pt]m&={\frac {\nu +1}{2}},\end{aligned}}}

كثافة هذه العائلة المقيدة ذات المعلمة الواحدة هي توزيع t القياسي للطالب :

ص(x)=1νب(ν2،12)(1+x2ν)-ν+12،{\displaystyle p(x)={\frac {1}{{\sqrt {\nu }}\,\operatorname {\mathrm {B} } \left({\frac {\nu }{2}},{\frac {1}{2}}\right)}}\left(1+{\frac {x^{2}}{\nu }}\right)^{-{\frac {\nu +1}{2}}},}

الحالة الثانية، التمييز غير السالب

إذا كانت الدالة التربيعية (2) لها مميز غير سالب (ب12-4ب2ب00{\displaystyle b_{1}^{2}-4b_{2}b_{0}\geq 0}، ولها جذور حقيقية a 1 و a 2 (ليست بالضرورة متميزة):

أ1=-ب1-ب12-4ب2ب02ب2،أ2=-ب1+ب12-4ب2ب02ب2.{\displaystyle {\begin{aligned}a_{1}&={\frac {-b_{1}-{\sqrt {b_{1}^{2}-4b_{2}b_{0}}}}{2b_{2}}},\\[5pt]a_{2}&={\frac {-b_{1}+{\sqrt {b_{1}^{2}-4b_{2}b_{0}}}}{2b_{2}}}.\end{aligned}}}

في حالة وجود جذور حقيقية، يمكن كتابة الدالة التربيعية (2) على النحو التالي:

ب(x)=ب2(x-أ1)(x-أ2)،{\displaystyle B(x)=b_{2}(x-a_{1})(x-a_{2}),}

وبالتالي فإن حل المعادلة التفاضلية هو

ص(x)خبرة(-1ب2x-أ(x-أ1)(x-أ2)دx).{\displaystyle p(x)\propto \exp \left(-{\frac {1}{b_{2}}}\int {\frac {x-a}{(x-a_{1})(x-a_{2})}}\,dx\right).}

أطلق بيرسون (1895، ص  362) على هذه الحالة اسم "الحالة اللوغاريتمية"، لأن التكامل

x-أ(x-أ1)(x-أ2)دx=(أ1-أ)ln(x-أ1)-(أ2-أ)ln(x-أ2)أ1-أ2+ج{\displaystyle \int {\frac {x-a}{(x-a_{1})(x-a_{2})}}\,dx={\frac {(a_{1}-a)\ln(x-a_{1})-(a_{2}-a)\ln(x-a_{2})}{a_{1}-a_{2}}}+C}

لا تتضمن سوى دالة اللوغاريتم وليس دالة الظل العكسي كما في الحالة السابقة.

باستخدام الاستبدال

ν=1ب2(أ1-أ2)،{\displaystyle \nu ={\frac {1}{b_{2}(a_{1}-a_{2})}},}

نحصل على الحل التالي للمعادلة التفاضلية (1):

ص(x)(x-أ1)-ν(أ1-أ)(x-أ2)ν(أ2-أ).{\displaystyle p(x)\propto (x-a_{1})^{-\nu (a_{1}-a)}(x-a_{2})^{\nu (a_{2}-a)}.}

بما أن هذه الكثافة معروفة فقط حتى ثابت تناسب خفي، فإنه يمكن تغيير هذا الثابت وكتابة الكثافة على النحو التالي:

ص(x)(1-xأ1)-ν(أ1-أ)(1-xأ2)ν(أ2-أ).{\displaystyle p(x)\propto \left(1-{\frac {x}{a_{1}}}\right)^{-\nu (a_{1}-a)}\left(1-{\frac {x}{a_{2}}}\right)^{\nu (a_{2}-a)}.}

توزيع بيرسون من النوع الأول

ينشأ توزيع بيرسون من النوع الأول ( وهو تعميم لتوزيع بيتا إلى منطقة دعم محدودة أكثر عمومية) عندما تكون جذور المعادلة التربيعية (2) ذات إشارات متعاكسة، أيأ1<0<أ2{\displaystyle a_{1}<0<a_{2}}إذن، يكون الحل p مدعومًا على الفترة(أ1،أ2){\displaystyle (a_{1},a_{2})}قم بتطبيق الاستبدال

x=أ1+y(أ2-أ1)،{\displaystyle x=a_{1}+y(a_{2}-a_{1}),}

أين0<y<1{\displaystyle 0<y<1}، مما ينتج عنه حل بدلالة y مدعوم على الفترة (0، 1):

ص(y)(أ1-أ2أ1y)(-أ1+أ)ν(أ2-أ1أ2(1-y))(أ2-أ)ν.{\displaystyle p(y)\propto \left({\frac {a_{1}-a_{2}}{a_{1}}}y\right)^{(-a_{1}+a)\nu }\left({\frac {a_{2}-a_{1}}{a_{2}}}(1-y)\right)^{(a_{2}-a)\nu }.}

يمكن تعريف ما يلي:

م1=أ-أ1ب2(أ1-أ2)،م2=أ-أ2ب2(أ2-أ1).{\displaystyle {\begin{aligned}m_{1}&={\frac {a-a_{1}}{b_{2}(a_{1}-a_{2})}},\\[5pt]m_{2}&={\frac {a-a_{2}}{b_{2}(a_{2}-a_{1})}}.\end{aligned}}}

بإعادة تجميع الثوابت والمعاملات، يتبسط هذا إلى

ص(y)yم1(1-y)م2،{\displaystyle p(y)\propto y^{m_{1}}(1-y)^{m_{2}},}

هكذاx-λ-أ1أ2-أ1{\displaystyle {\frac {x-\lambda -a_{1}}{a_{2}-a_{1}}}}يتبع توزيع بيتاب(م1+1،م2+1){\displaystyle \mathrm {B} (m_{1}+1,m_{2}+1)}معλ=μ1-(أ2-أ1)م1+1م1+م2+2-أ1{\displaystyle \lambda =\mu _{1}-(a_{2}-a_{1}){\frac {m_{1}+1}{m_{1}+m_{2}+2}}-a_{1}}. اتضح أن m 1 ، m 2 > −1 ضروري وكافٍ لكي تكون p دالة كثافة احتمالية مناسبة.

توزيع بيرسون من النوع الثاني

يُعد توزيع بيرسون من النوع الثاني حالة خاصة من عائلة توزيعات بيرسون من النوع الأول، ويقتصر على التوزيعات المتماثلة. باستخدام الصيغ الواردة في قسم النوع الأول، معم1=م2=م{\displaystyle m_{1}=m_{2}=m}و-أ1=أ2=أ{\displaystyle -a_{1}=a_{2}=a}، على الفترة (−a, a) يمكن كتابتها على النحو التالي:

ص(x)(1-x2أ2)م.{\displaystyle p(x)\propto \left(1-{\frac {x^{2}}{a^{2}}}\right)^{m}.}

أو مع

x=-أ+2yأ،{\displaystyle x=-a+2ya,}

y{\displaystyle y}يتم توزيعها وفقًا لتوزيع بيتا على الفترة (0، 1)،

ص(y)(1-4(y-12)2)مyم(1-y)م.{\displaystyle p(y)\propto \left(1-4\left(y-{\frac {1}{2}}\right)^{2}\right)^{m}\propto y^{m}(1-y)^{m}.}

مع ثابت تناسب مناسب، تصبح دالة كثافة الاحتمال

ص(y)=yم(1-y)مΓ(2م+2)Γ(م+1)2.{\displaystyle p(y)=y^{m}(1-y)^{m}{\frac {\Gamma (2m+2)}{\Gamma (m+1)^{2}}}.}

توزيع بيرسون من النوع الثالث

تعريف

λ=μ1+ب0ب1-(م+1)ب1،{\displaystyle \lambda =\mu _{1}+{\frac {b_{0}}{b_{1}}}-(m+1)b_{1},}

ب0+ب1(x-λ){\displaystyle b_{0}+b_{1}(x-\lambda )}يكونجاما(م+1،ب12){\displaystyle \operatorname {Gamma} (m+1,b_{1}^{2})}. توزيع بيرسون من النوع الثالث هو توزيع جاما أو توزيع مربع كاي .

توزيع بيرسون من النوع الخامس

تحديد معايير جديدة:

ج1=ب12ب2،λ=μ1-أ-ج11-2ب2،{\displaystyle {\begin{aligned}C_{1}&={\frac {b_{1}}{2b_{2}}},\\\lambda &=\mu _{1}-{\frac {a-C_{1}}{1-2b_{2}}},\end{aligned}}}

x-λ{\displaystyle x-\lambda }يتبعغاما المعكوسة(1ب2-1،أ-ج1ب2){\displaystyle \operatorname {InverseGamma} ({\frac {1}{b_{2}}}-1,{\frac {a-C_{1}}{b_{2}}})}توزيع بيرسون من النوع الخامس هو توزيع جاما معكوس .

توزيع بيرسون من النوع السادس

تعريف

λ=μ1+(أ2-أ1)م2+1م2+م1+2-أ2،{\displaystyle \lambda =\mu _{1}+(a_{2}-a_{1}){\frac {m_{2}+1}{m_{2}+m_{1}+2}}-a_{2},}

x-λ-أ2أ2-أ1{\displaystyle {\frac {x-\lambda -a_{2}}{a_{2}-a_{1}}}}يتبعβ(م2+1،-م2-م1-1){\displaystyle \beta ^{\prime }(m_{2}+1,-m_{2}-m_{1}-1)}. توزيع بيرسون من النوع السادس هو توزيع بيتا برايم أو توزيع F.

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

تشمل عائلة بيرسون التوزيعات التالية، من بين توزيعات أخرى:

اعتبارًا من عام 2025، فإن النوعين الوحيدين اللذين ليس لهما اسم هما النوع الرابع (انظر أعلاه ) والنوع الثاني عشر ( توزيع بيتا معα+β=2{\displaystyle \alpha +\beta =2}).

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

أمثلة على البدائل الحديثة لمخطط بيرسون للالتواء مقابل التفرطح هي: (1) https://github.com/SchildCode/PearsonPlot و (2) "مخطط كولين وفري" في التطبيق الإحصائي R.

التطبيقات

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

في الولايات المتحدة، يُعد توزيع لوغاريتم جاما (المسمى تاريخياً بتوزيع لوغاريتم بيرسون الثالث) التوزيع الافتراضي لتحليل تكرار الفيضانات. [ 4 ]

ظهرت مؤخراً بدائل لتوزيعات بيرسون تتميز بمرونة أكبر وسهولة أكبر في ملاءمتها للبيانات. انظر إلى التوزيعات الميتالوجية .

ملحوظات

  1. ميلر، جيف؛ وآخرون  . (9 يوليو 2006). "توزيع بيتا" . أقدم الاستخدامات المعروفة لبعض مصطلحات الرياضيات . تم الاسترجاع في 9 ديسمبر 2006 .
  2. ميلر، جيف؛ وآخرون . (7 ديسمبر 2006). "توزيع جاما" . أقدم الاستخدامات المعروفة لبعض مصطلحات الرياضيات . تم الاسترجاع في 9 ديسمبر 2006 . 
  3. Ord JK (1972) ص. 2
  4. "إرشادات لتحديد تواتر تدفق الفيضانات" (ملف PDF) . هيئة المسح الجيولوجي الأمريكية للمياه . مارس 1982. تم الاطلاع عليه بتاريخ 14 يونيو 2019 .

مصادر

المصادر الأولية

مصادر ثانوية

مراجع

  • إلديرتون، السير دبليو بي، جونسون، إن إل (1969) أنظمة منحنيات التردد . مطبعة جامعة كامبريدج.
  • Ord JK (1972) عائلات التوزيعات التكرارية . Griffin، لندن.