توزيع لابلاس

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

التعريفات

دالة كثافة الاحتمال

للمتغير العشوائيلابلاس(μ،ب){\displaystyle \operatorname {Laplace} (\mu ,b)}التوزيع إذا كانت دالة كثافة احتماله هي

و(x|μ،ب)=12بهـ-|x-μ|ب،{\displaystyle f(x\mid \mu ,b)={\frac {1}{2b}}e^{-{\frac {|x-\mu |}{b}}},}

أينμ{\displaystyle \mu }هو مُعامل الموقع ، وب>0{\displaystyle b>0}، والذي يُشار إليه أحيانًا باسم "التنوع"، هو مُعامل قياس . إذاμ=0{\displaystyle \mu =0}وب=1{\displaystyle b=1}، يكون نصف الخط الموجب بالضبط توزيعًا أسيًا مضروبًا في 1/2. [ 2 ]

تُشبه دالة كثافة الاحتمال لتوزيع لابلاس دالة التوزيع الطبيعي ؛ إلا أنه بينما يُعبَّر عن التوزيع الطبيعي بدلالة مربع الفرق عن المتوسط،μ{\displaystyle \mu }تُعبّر دالة كثافة لابلاس عن الفرق المطلق عن المتوسط. ونتيجةً لذلك، يتميز توزيع لابلاس بذيل أكثر سمكًا من التوزيع الطبيعي. وهو حالة خاصة من التوزيع الطبيعي المعمم والتوزيع الزائدي . ومن التوزيعات المتناظرة المستمرة التي لها ذيول أسية، مثل توزيع لابلاس، ولكن دوال كثافة احتمالها قابلة للتفاضل عند المنوال، التوزيع اللوجستي ، وتوزيع القاطع الزائدي ، وتوزيع شامبرنون .

دالة التوزيع التراكمي

يسهل حساب تكامل توزيع لابلاس (إذا تم التمييز بين حالتين متماثلتين) نظرًا لاستخدام دالة القيمة المطلقة . ودالة التوزيع التراكمي له هي كما يلي: F(x)=-xو(u)دu={12خبرة(x-μب)لو x<μ1-12خبرة(-x-μب)لو xμ=12+12علامة(x-μ)(1-خبرة(-|x-μ|ب)).{\displaystyle {\begin{aligned}F(x)&=\int _{-\infty }^{x}\!\!f(u)\,\mathrm {d} u={\begin{cases}{\frac {1}{2}}\exp \left({\frac {x-\mu }{b}}\right)&{\mbox{if }}x<\mu \\1-{\frac {1}{2}}\exp \left(-{\frac {x-\mu }{b}}\right)&{\mbox{if }}x\geq \mu \end{cases}}\\&={\tfrac {1}{2}}+{\tfrac {1}{2}}\operatorname {sgn}(x-\mu )\left(1-\exp \left(-{\frac {|x-\mu |}{b}}\right)\right).\end{aligned}}}

دالة التوزيع التراكمي العكسي معطاة بالصيغة التالية

F-1(ص)=μ-بعلامة(ص-0.5)ln(1-2|ص-0.5|).{\displaystyle F^{-1}(p)=\mu -b\,\operatorname {sgn}(p-0.5)\,\ln(1-2|p-0.5|).}

ملكيات

لحظات

μر=12ك=0ر[ر!(ر-ك)!بكμ(ر-ك){1+(-1)ك}].{\displaystyle \mu _{r}'={\frac {1}{2}}\sum _{k=0}^{r}\left[{\frac {r!}{(rk)!}}b^{k}\mu ^{(rk)}\left\{1+(-1)^{k}\right\}\right].}

  • لوXلابلاس(μ،ب){\displaystyle X\sim {\textrm {Laplace}}(\mu ,b)}ثمكX+جلابلاس(كμ+ج،|ك|ب){\displaystyle kX+c\sim {\textrm {Laplace}}(k\mu +c,|k|b)}.
  • لوXلابلاس(0،1){\displaystyle X\sim {\textrm {Laplace}}(0,1)}ثمبXلابلاس(0،ب){\displaystyle bX\sim {\textrm {Laplace}}(0,b)}.
  • لوXلابلاس(0،ب){\displaystyle X\sim {\textrm {Laplace}}(0,b)}ثم|X|النمو الأسي(ب-1){\displaystyle \left|X\right|\sim {\textrm {Exponential}}\left(b^{-1}\right)}( التوزيع الأسي ).
  • لوX،Yالنمو الأسي(λ){\displaystyle X,Y\sim {\textrm {Exponential}}(\lambda )}ثمX-Yلابلاس(0،λ-1){\displaystyle XY\sim {\textrm {Laplace}}\left(0,\lambda ^{-1}\right)}
  • لوXلابلاس(μ،ب){\displaystyle X\sim {\textrm {Laplace}}(\mu ,b)}ثم|X-μ|النمو الأسي(ب-1){\displaystyle \left|X-\mu \right|\sim {\textrm {Exponential}}(b^{-1})}.
  • لوXلابلاس(μ،ب){\displaystyle X\sim {\textrm {Laplace}}(\mu ,b)}ثمXEPD(μ،ب،1){\displaystyle X\sim {\textrm {EPD}}(\mu ,b,1)}( توزيع القوة الأسي ).
  • لوX1،...،X4شمال(0،1){\displaystyle X_{1},...,X_{4}\sim {\textrm {N}}(0,1)}( التوزيع الطبيعي ) إذنX1X2-X3X4لابلاس(0،1){\displaystyle X_{1}X_{2}-X_{3}X_{4}\sim {\textrm {Laplace}}(0,1)}و(X12-X22+X32-X42)/2لابلاس(0،1){\displaystyle (X_{1}^{2}-X_{2}^{2}+X_{3}^{2}-X_{4}^{2})/2\sim {\textrm {Laplace}}(0,1)}.
  • لوXأنالابلاس(μ،ب){\displaystyle X_{i}\sim {\textrm {Laplace}}(\mu ,b)}ثم2بأنا=1ن|Xأنا-μ|χ2(2ن){\displaystyle {\frac {\displaystyle 2}{b}}\sum _{i=1}^{n}|X_{i}-\mu |\sim \chi ^{2}(2n)}( توزيع مربع كاي ).
  • لوX،Yلابلاس(μ،ب){\displaystyle X,Y\sim {\textrm {Laplace}}(\mu ,b)}ثم|X-μ||Y-μ|F(2،2){\displaystyle {\tfrac {|X-\mu |}{|Y-\mu |}}\sim \operatorname {F} (2,2)}( توزيع F )
  • لوX،Yيو(0،1){\displaystyle X,Y\sim {\textrm {U}}(0,1)}( التوزيع المنتظم ) إذنسجل(X/Y)لابلاس(0،1){\displaystyle \log(X/Y)\sim {\textrm {Laplace}}(0,1)}.
  • لوXالنمو الأسي(λ){\displaystyle X\sim {\textrm {Exponential}}(\lambda )}وYبرنولي(0.5){\displaystyle Y\sim {\textrm {Bernoulli}}(0.5)}( توزيع برنولي ) مستقل عنX{\displaystyle X}، ثمX(2Y-1)لابلاس(0،λ-1){\displaystyle X(2Y-1)\sim {\textrm {Laplace}}\left(0,\lambda ^{-1}\right)}.
  • لوXالنمو الأسي(λ){\displaystyle X\sim {\textrm {Exponential}}(\lambda )}وYالنمو الأسي(ν){\displaystyle Y\sim {\textrm {الأسي}}(\nu )}مستقل عنX{\displaystyle X}، ثمλX-νYلابلاس(0،1){\displaystyle \lambda X-\nu Y\sim {\textrm {Laplace}}(0,1)}
  • لوX{\displaystyle X}لديها توزيع رادماخر وYالنمو الأسي(λ){\displaystyle Y\sim {\textrm {Exponential}}(\lambda )}ثمXYلابلاس(0،1/λ){\displaystyle XY\sim {\textrm {Laplace}}(0,1/\lambda )}.
  • لوVالنمو الأسي(1){\displaystyle V\sim {\textrm {Exponential}}(1)}وZشمال(0،1){\displaystyle Z\sim N(0,1)}مستقل عنV{\displaystyle V}، ثمX=μ+ب2VZلأصلأجهـ(μ،ب){\displaystyle X=\mu +b{\sqrt {2V}}Z\sim \mathrm {Laplace} (\mu ,b)}.
  • لوXمستقر هندسيًا(2،0،λ،0){\displaystyle X\sim {\textrm {GeometricStable}}(2,0,\lambda ,0)}( التوزيع الهندسي المستقر ) ثمXلابلاس(0،λ){\displaystyle X\sim {\textrm {Laplace}}(0,\lambda )}.
  • يُعد توزيع لابلاس حالة حدية للتوزيع الزائدي .
  • لوX|Yشمال(μ،Y2){\displaystyle X|Y\sim {\textrm {N}}(\mu ,Y^{2})}معYرايلي(ب){\displaystyle Y\sim {\textrm {Rayleigh}}(b)}( توزيع رايلي )Xلابلاس(μ،ب){\displaystyle X\sim {\textrm {Laplace}}(\mu ,b)}لاحظ أنه إذاYرايلي(ب){\displaystyle Y\sim {\textrm {Rayleigh}}(b)}، ثمY2جاما(1،2ب2){\displaystyle Y^{2}\sim {\textrm {Gamma}}(1,2b^{2})}معهـ(Y2)=2ب2{\displaystyle {\textrm {E}}(Y^{2})=2b^{2}}وهذا بدوره يساوي التوزيع الأسيخبرة(1/(2ب2)){\displaystyle {\textrm {Exp}}(1/(2b^{2}))}.
  • بفرض عدد صحيحن1{\displaystyle n\geq 1}، لوXأنا،YأناΓ(1ن،ب){\displaystyle X_{i},Y_{i}\sim \Gamma \left({\frac {1}{n}},b\right)}( توزيع جاما ، باستخدامك،θ{\displaystyle k,\theta }(التوصيف)، ثمأنا=1ن(μن+Xأنا-Yأنا)لابلاس(μ،ب){\displaystyle \sum _{i=1}^{n}\left({\frac {\mu }{n}}+X_{i}-Y_{i}\right)\sim {\textrm {Laplace}}(\mu ,b)}( قابلية القسمة اللانهائية ) [ 3 ]
  • إذا كان X له توزيع لابلاس، فإن Y = e X له توزيع لوغاريتمي لابلاس ؛ وعلى العكس من ذلك، إذا كان X له توزيع لوغاريتمي لابلاس، فإن لوغاريتمه له توزيع لابلاس.

احتمال أن يكون عدد لابلاس أكبر من عدد آخر

يترك X،Y{\displaystyle X,Y}أن تكون متغيرات عشوائية مستقلة من نوع لابلاس:Xلابلاس(μX،بX){\displaystyle X\sim {\textrm {Laplace}}(\mu _{X},b_{X})}وYلابلاس(μY،بY){\displaystyle Y\sim {\textrm {Laplace}}(\mu _{Y},b_{Y})}ونريد أن نحسبP(X>Y){\displaystyle P(X>Y)}.

احتماليةP(X>Y){\displaystyle P(X>Y)}يمكن اختزالها (باستخدام الخصائص أدناه) إلىP(μ+بZ1>Z2){\displaystyle P(\mu +bZ_{1}>Z_{2})}، أينZ1،Z2لابلاس(0،1){\displaystyle Z_{1},Z_{2}\sim {\textrm {Laplace}}(0,1)}هذا الاحتمال يساوي

P(μ+بZ1>Z2)={ب2هـμ/ب-هـμ2(ب2-1)،متى μ<01-ب2هـ-μ/ب-هـ-μ2(ب2-1)،متى μ>0{\displaystyle P(\mu +bZ_{1}>Z_{2})={\begin{cases}{\frac {b^{2}e^{\mu /b}-e^{\mu }}{2(b^{2}-1)}},&{\text{when }}\mu <0\\1-{\frac {b^{2}e^{-\mu /b}-e^{-\mu }}{2(b^{2}-1)}},&{\text{when }}\mu >0\\\end{cases}}}

متىب=1{\displaystyle b=1}يتم استبدال كلا التعبيرين بنهايتهما كما يلي:ب1{\displaystyle b\to 1}:

P(μ+Z1>Z2)={هـμ(2-μ)4،متى μ<01-هـ-μ(2+μ)4،متى μ>0{\displaystyle P(\mu +Z_{1}>Z_{2})={\begin{cases}e^{\mu }{\frac {(2-\mu )}{4}},&{\text{when }}\mu <0\\1-e^{-\mu }{\frac {(2+\mu )}{4}},&{\text{when }}\mu >0\\\end{cases}}}

لحساب الحالة لـμ>0{\displaystyle \mu >0}، لاحظ أنP(μ+Z1>Z2)=1-P(μ+Z1<Z2)=1-P(-μ-Z1>-Z2)=1-P(-μ+Z1>Z2){\displaystyle P(\mu +Z_{1}>Z_{2})=1-P(\mu +Z_{1}<Z_{2})=1-P(-\mu -Z_{1}>-Z_{2})=1-P(-\mu +Z_{1}>Z_{2})}

منذZ-Z{\displaystyle Z\sim -Z}متىZلابلاس(0،1){\displaystyle Z\sim {\textrm {Laplace}}(0,1)}.

العلاقة بالتوزيع الأسي

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

لنفترض وجود متغيرين عشوائيين مستقلين ومتطابقين في التوزيع.X،Yالنمو الأسي(λ){\displaystyle X,Y\sim {\textrm {Exponential}}(\lambda )}الدوال المميزة لـX،-Y{\displaystyle X,-Y}نكون

λ-أنات+λ،λأنات+λ{\displaystyle {\frac {\lambda }{-it+\lambda }},\quad {\frac {\lambda }{it+\lambda }}}

على التوالي. عند ضرب هذه الدوال المميزة (المكافئة للدالة المميزة لمجموع المتغيرات العشوائية)X+(-Y){\displaystyle X+(-Y)}والنتيجة هي

λ2(-أنات+λ)(أنات+λ)=λ2ت2+λ2.{\displaystyle {\frac {\lambda ^{2}}{(-it+\lambda )(it+\lambda )}}={\frac {\lambda ^{2}}{t^{2}+\lambda ^{2}}}.}

وهذا هو نفسه الدالة المميزة لـZلابلاس(0،1/λ){\displaystyle Z\sim {\textrm {Laplace}}(0,1/\lambda )}، وهو

11+ت2λ2.{\displaystyle {\frac {1}{1+{\frac {t^{2}}{\lambda ^{2}}}}}.}

توزيعات سارجان

تُعدّ توزيعات سارجان نظامًا من التوزيعات، ويُعتبر توزيع لابلاس أحد أعضائه الأساسيين.ص{\displaystyle p}توزيع سارجان من الرتبة th له كثافة [ 4 ] [ 5 ]وص(x)=12خبرة(-α|x|)1+ج=1صβجαج|x|ج1+ج=1صج!βج،{\displaystyle f_{p}(x)={\tfrac {1}{2}}\exp(-\alpha |x|){\frac {\displaystyle 1+\sum _{j=1}^{p}\beta _{j}\alpha ^{j}|x|^{j}}{\displaystyle 1+\sum _{j=1}^{p}j!\beta _{j}}},} بالنسبة للمعاملاتα0،βج0{\displaystyle \alpha \geq 0,\beta _{j}\geq 0}نتائج توزيع لابلاس لـص=0{\displaystyle p=0}.

الاستدلال الإحصائي

منحن{\displaystyle n}عينات مستقلة وموزعة بشكل متطابقx1،x2،...،xن{\displaystyle x_{1},x_{2},...,x_{n}}، مُقدِّر الاحتمال الأقصى (MLE) لـμ{\displaystyle \mu }هو الوسيط العيني ، [ 6 ]μ^=مهـد(x).{\displaystyle {\hat {\mu }}=\mathrm {med} (x).}

مقدِّر الاحتمال الأقصى لـب{\displaystyle b}هو متوسط ​​الانحراف المطلق عن الوسيط،ب^=1نأنا=1ن|xأنا-μ^|.{\displaystyle {\hat {b}}={\frac {1}{n}}\sum _{i=1}^{n}|x_{i}-{\hat {\mu }}|.} يكشف ذلك عن وجود صلة بين توزيع لابلاس وأقل الانحرافات المطلقة . ويمكن تطبيق تصحيح للعينات الصغيرة على النحو التالي: ب^*=ب^ن/(ن-2){\displaystyle {\hat {b}}^{*}={\hat {b}}\cdot n/(n-2)} (انظر: التوزيع الأسي#تقدير المعلمات ).

التواجد والتطبيقات

تم استخدام توزيع لابلاس في التعرف على الكلام لنمذجة المعلومات المسبقة حول معاملات DFT [ 7 ] وفي ضغط صور JPEG لنمذجة معاملات AC [ 8 ] التي تم إنشاؤها بواسطة DCT .

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

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

بافتراض وجود متغير عشوائييو{\displaystyle U}تم سحبها من التوزيع المنتظم في الفترة(-1/2،1/2){\displaystyle \left(-1/2,1/2\right)}المتغير العشوائي

X=μ-بعلامة(يو)ln(1-2|يو|){\displaystyle X=\mu -b\,\operatorname {sgn}(U)\,\ln(1-2|U|)}

له توزيع لابلاس بمعاملاتμ{\displaystyle \mu }وب{\displaystyle b}وهذا يتبع من دالة التوزيع التراكمي العكسي المذكورة أعلاه.

ألابلاس(0،ب){\displaystyle {\textrm {Laplace}}(0,b)}يمكن أيضًا إنشاء المتغير كفرق بين متغيرين مستقلين ومتطابقين في التوزيع.النمو الأسي(1/ب){\displaystyle {\textrm {Exponential}}(1/b)}المتغيرات العشوائية. أو بصورة مكافئة،لابلاس(0،1){\displaystyle {\textrm {Laplace}}(0,1)}ويمكن أيضًا توليدها كلوغاريتم نسبة متغيرين عشوائيين مستقلين ومتطابقين وموزعين توزيعًا منتظمًا.

تاريخ

يُشار إلى هذا التوزيع غالبًا باسم "قانون لابلاس الأول للأخطاء". وقد نشره عام 1774، حيث صوّر تكرار الخطأ كدالة أسية لقيمته بمجرد تجاهل إشارته. لاحقًا، استبدل لابلاس هذا النموذج بـ"قانونه الثاني للأخطاء"، المستند إلى التوزيع الطبيعي، بعد اكتشاف نظرية النهاية المركزية . [ 13 ] [ 14 ]

نشر كينز ورقة بحثية في عام 1911 استنادًا إلى أطروحته السابقة حيث أظهر أن توزيع لابلاس يقلل من الانحراف المطلق عن الوسيط. [ 15 ]

انظر أيضاً

مراجع

  1. تشاتامفيلي، راجان؛ شانموغام، رامالينغام (2021)، "توزيع لابلاس" ، التوزيعات المستمرة في الهندسة والعلوم التطبيقية - الجزء الثاني ، تشام: دار نشر سبرينغر الدولية، ص 189-199 ، doi : 10.1007/978-3-031-02435-1_4 ، ISBN  978-3-031-01307-2تم الاطلاع عليه بتاريخ 2025-04-04
  2. هوانغ، يونفي؛ وآخرون (2022). "الاستدلال المتفرق والتعلم النشط للمعادلات التفاضلية العشوائية من البيانات" . التقارير العلمية . 12 (1): 21691. arXiv : 2203.11010 . Bibcode : 2022NatSR..1221691H . doi : 10.1038/ s41598-022-25638-9 . PMC 9755218. PMID 36522347 .   
  3. 1 2 كوتز، صموئيل؛ كوزوبوفسكي، توماش ج.؛ بودغورسكي، كريستوف (2001). توزيع لابلاس وتعميماته: إعادة نظر مع تطبيقات في الاتصالات والاقتصاد والهندسة والتمويل . بيركهاوزر. ص 23 (الفرضية 2.2.2، المعادلة 2.2.8). ISBN  9780817641665.
  4. إيفريت، بي إس (2002) قاموس كامبريدج للإحصاء ، مطبعة جامعة كامبريدج. رقم ISBN 0-521-81099-X
  5. جونسون، إن إل، كوتز، إس، بالاكريشنان، إن (1994) التوزيعات أحادية المتغير المستمرة ، وايلي. ISBN 0-471-58495-9ص 60
  6. روبرت م. نورتون (مايو 1984). "التوزيع الأسي المزدوج: استخدام حساب التفاضل والتكامل لإيجاد مُقدِّر الاحتمال الأقصى". الإحصائي الأمريكي . 38 (2). الجمعية الإحصائية الأمريكية: 135-136 . doi : 10.2307/2683252 . JSTOR 2683252 . 
  7. إلتوفت، ت.؛ تايسو كيم؛ تي-وون لي (2006). "حول توزيع لابلاس متعدد المتغيرات" (ملف PDF) . رسائل معالجة الإشارات IEEE . 13 (5): 300-303 . رمز Bibcode : 2006ISPL...13..300E . doi : 10.1109/LSP.2006.870353 . S2CID 1011487. مؤرشف من الأصل (ملف PDF) بتاريخ 2013-06-06 . تم الاسترجاع بتاريخ 2012-07-04 . 
  8. مينغيلون، ج.؛ بوجول، ج. (2001). "نمذجة خطأ التكميم الموحد لمعيار JPEG مع تطبيقات على أوضاع التشغيل التسلسلي والتدريجي" (ملف PDF) . مجلة التصوير الإلكتروني . 10 (2): 475-485 . Bibcode : 2001JEI....10..475M . doi : 10.1117/1.1344592 . hdl : 10609/6263 .
  9. باردو، سكوت (2020). التحليل الإحصائي لأساليب البيانات التجريبية في العلوم التطبيقية . سبرينغر. ص 58. ISBN  978-3-030-43327-7.
  10. كو، إس جي (8 أغسطس 2002). "نموذج الانتشار القفزي لتسعير الخيارات" . مجلة علوم الإدارة . 48 (8): 1086-1101 . doi : 10.1287/mnsc.48.8.1086.166 . JSTOR 822677. تاريخ الاسترجاع: 1 مارس 2022 . 
  11. تشين، جيان (2018). طريقة تسعير خيارات التوازن العام: دراسة نظرية وتجريبية . سبرينغر. ص 70. ISBN  9789811074288.
  12. مجموعة من التوزيعات المركبة
  13. ^ لابلاس، ملاحظة. (1774). احتفظ باحتمالية الأسباب للأحداث. مذكرات الأكاديمية الملكية للعلوم المقدمة للغواصين سافان، 6، 621-656
  14. ويلسون، إدوين بيدويل (1923). "القانونان الأول والثاني للخطأ". مجلة الجمعية الإحصائية الأمريكية . 18 (143). إنفورما المملكة المتحدة المحدودة: 841-851 . doi : 10.1080/01621459.1923.10502116 . ISSN 0162-1459 . المجال العامتتضمن هذه المقالة نصًا من هذا المصدر، وهو متاح للعموم .
  15. كينز، ج. م. (1911). "المتوسطات الرئيسية وقوانين الخطأ المؤدية إليها" . مجلة الجمعية الإحصائية الملكية . 74 (3). JSTOR: 322–331 . doi : 10.2307/2340444 . ISSN 0952-8385 . JSTOR 2340444 .