توزيع مستقر

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

من بين المعايير الأربعة التي تحدد العائلة، تم التركيز بشكل أكبر على معيار الاستقرار.α{\displaystyle \alpha }(انظر اللوحة). التوزيعات المستقرة لها0<α2{\displaystyle 0<\alpha \leq 2}، مع حد أعلى يتوافق مع التوزيع الطبيعي ، ويقترب من دالة ديراك دلتا في النهاية عندماα0{\displaystyle \alpha \rightarrow 0}تتميز التوزيعات بتباين غير محدد لـα<2{\displaystyle \alpha <2}، ومتوسط ​​غير محدد لـα1{\displaystyle \alpha \leq 1}.

تكمن أهمية التوزيعات الاحتمالية المستقرة في كونها " جاذبات " لمجاميع معيارية مناسبة لمتغيرات عشوائية مستقلة ومتطابقة التوزيع ( iid ). يُعرّف التوزيع الطبيعي عائلة من التوزيعات المستقرة. وبحسب نظرية النهاية المركزية الكلاسيكية ، فإن المجموع المعياري المناسب لمجموعة من المتغيرات العشوائية، لكل منها تباين محدود، سيميل نحو التوزيع الطبيعي مع ازدياد عدد المتغيرات. وبدون افتراض التباين المحدود، قد تكون النهاية توزيعًا مستقرًا غير طبيعي. وقد أشار ماندلبروت إلى هذه التوزيعات باسم "توزيعات باريتو المستقرة" [ 3 ] [ 4 ] [ 5 ] نسبةً إلى فيلفريدو باريتو . وعلى وجه الخصوص، أشار إلى تلك التوزيعات ذات الانحراف الأقصى في الاتجاه الموجب.1<α<2{\displaystyle 1<\alpha <2}[ 1 ] والتي اعتبرها وصفاً أفضل لأسعار الأسهم والسلع من التوزيعات الطبيعية . [ 6 ]

تعريف

يُعتبر التوزيع غير المنحط توزيعًا مستقرًا إذا حقق الخاصية التالية:

ليكن X1 و X2 نسختين مستقلتين من المتغير العشوائي X. يُقال إن X مستقر إذا كان للمتغير العشوائي aX1 + bX2 ، لأي ثابتين a > 0 و b > 0 ، نفس توزيع cX + d لبعض الثوابت c > 0 و d . ويُقال إن التوزيع مستقر تمامًا إذا تحقق هذا الشرط مع d = 0. [ 7 ]

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

تشكل هذه التوزيعات عائلة من أربعة معلمات من التوزيعات الاحتمالية المستمرة ، والتي يتم تحديدها بواسطة معلمات الموقع والمقياس μ و c على التوالي، ومعلمتي الشكل.β{\displaystyle \beta }وα{\displaystyle \alpha }، والتي تتوافق تقريبًا مع مقاييس عدم التماثل والتركيز، على التوالي (انظر الأشكال).

الدالة المميزةφ{\displaystyle \varphi }لتوزيع احتمالي بدالة كثافةو{\displaystyle f}هو تحويل فورييه لـو.{\displaystyle f.}تكون دالة الكثافة عندئذٍ هي التحويل العكسي لفورييه للدالة المميزة: [ 8 ]φ(ت)=-و(x)هـأناxتدx.{\displaystyle \varphi (t)=\int _{-\infty }^{\infty }f(x)e^{ixt}\,dx.}

على الرغم من أنه لا يمكن كتابة دالة كثافة الاحتمال لتوزيع مستقر عام تحليليًا، إلا أنه يمكن التعبير عن الدالة المميزة العامة تحليليًا. يُطلق على المتغير العشوائي X اسم "مستقر" إذا أمكن كتابة دالته المميزة على النحو التالي [ 7 ] [ 9 ].φ(ت؛α،β،ج،μ)=خبرة(أناتμ-|جت|α(1-أناβعلامة(ت)Φ)){\displaystyle \varphi (t;\alpha ,\beta ,c,\mu )=\exp \left(it\mu -|ct|^{\alpha }\left(1-i\beta \operatorname {sgn}(t)\Phi \right)\right)} حيث sgn ( t ) هي إشارة t وΦ={لون برونزي(πα2)α1-2πسجل|ت|α=1{\displaystyle \Phi ={\begin{cases}\tan \left({\frac {\pi \alpha }{2}}\right)&\alpha \neq 1\\-{\frac {2}{\pi }}\log |t|&\alpha =1\end{cases}}}μR هو معامل إزاحة،β[-1،1]{\displaystyle \beta \in [-1,1]}يُطلق على هذا المقياس اسم معامل الالتواء ، وهو مقياس لعدم التماثل. لاحظ أنه في هذا السياق، لا يكون مفهوم الالتواء المعتاد مُعرَّفًا بدقة، كما هو الحال بالنسبة لـα<2{\displaystyle \alpha <2}لا يسمح التوزيع بوجود لحظات من الدرجة الثانية أو أعلى ، والتعريف المعتاد للانحراف هو اللحظة المركزية الثالثة .

يرجع سبب استقرار هذا التوزيع إلى أن الدالة المميزة لمجموع متغيرين عشوائيين مستقلين تساوي حاصل ضرب الدالتين المميزتين المناظرتين لهما. فجمع متغيرين عشوائيين من توزيع مستقر يعطينا توزيعًا له نفس قيم .α{\displaystyle \alpha }وβ{\displaystyle \beta }، ولكن ربما تكون هناك قيم مختلفة لـ μ و c .

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

في أبسط الحالاتβ=0{\displaystyle \beta =0}، الدالة المميزة هي مجرد دالة أسية ممتدة ؛ التوزيع متناظر حول μ ويشار إليه باسم توزيع ألفا المستقر المتناظر (ليفي) ، وغالبًا ما يتم اختصاره إلى SαS .

متىα<1{\displaystyle \alpha <1}وβ=1{\displaystyle \beta =1}، يتم دعم التوزيع على [ μ ، ∞).

المعامل c > 0 هو عامل قياس يمثل مقياسًا لعرض التوزيع بينماα{\displaystyle \alpha }يمثل الأس أو الدليل للتوزيع ويحدد السلوك التقاربي للتوزيع.

المعلمات

إنّ تمثيل التوزيعات المستقرة ليس فريدًا. وقد قام نولان [ 10 ] بتلخيص 11 تمثيلًا واردًا في الأدبيات، وقدم صيغًا للتحويل. التمثيلان الأكثر شيوعًا هما التمثيل المذكور أعلاه (تمثيل نولان "1") والتمثيل الذي يليه مباشرة (تمثيل نولان "0").

إنّ النموذج البارامتري المذكور أعلاه هو الأسهل استخدامًا في العمل النظري، لكن كثافة احتماله غير متصلة في المعاملات عندα=1{\displaystyle \alpha =1}[ 11 ] إن استخدام المعلمات المستمرة، وهو أفضل للعمل العددي، هو [ 7 ]φ(ت؛α،β،γ،دلتا)=خبرة(أناتدلتا-|γت|α(1-أناβعلامة(ت)Φ)){\displaystyle \varphi (t;\alpha ,\beta ,\gamma ,\delta )=\exp \left(it\delta -|\gamma t|^{\alpha }\left(1-i\beta \operatorname {sgn}(t)\Phi \right)\right)} أين: Φ={(1-|γت|1-α)لون برونزي(πα2)α1-2πسجل|γت|α=1//

نطاقاتα{\displaystyle \alpha }وβ{\displaystyle \beta }هي نفسها كما كانت من قبل، يجب أن تكون γ (مثل c ) موجبة، ويجب أن تكون δ (مثل μ ) حقيقية.

في أي من الصيغتين، يمكن إجراء تحويل خطي للمتغير العشوائي للحصول على متغير عشوائي تكون كثافتهو(y؛α،β،1،0){\displaystyle f(y;\alpha ,\beta ,1,0)}في عملية تحديد المعلمات الأولى ، يتم ذلك عن طريق تعريف المتغير الجديد: y={x-μγα1x-μγ-β2πlnγα=1{\displaystyle y={\begin{cases}{\frac {x-\mu }{\gamma }}&\alpha \neq 1\\[1ex]{\frac {x-\mu }{\gamma }}-\beta {\frac {2}{\pi }}\ln \gamma &\alpha =1\end{cases}}}

بالنسبة للمعامل الثاني، استخدم ببساطة y=x-دلتاγ{\displaystyle y={\frac {x-\delta }{\gamma }}} مستقل عنα{\displaystyle \alpha }في المعايرة الأولى، إذا كان المتوسط ​​موجودًا (أي،α>1{\displaystyle \alpha >1}ثم يكون مساوياً لـ μ ، بينما في المعامل الثاني عندما يكون المتوسط ​​موجوداً يكون مساوياً لـدلتا-βγلون برونزي(πα2).{\displaystyle \delta -\beta \gamma \tan \left({\tfrac {\pi \alpha }{2}}\right).}

التوزيع

وبالتالي، يُحدد التوزيع المستقر بالمعلمات الأربع المذكورة أعلاه. ويمكن إثبات أن أي توزيع مستقر غير منحل له دالة كثافة سلسة (قابلة للتفاضل بلا حدود). [ 7 ] إذاو(x؛α،β،ج،μ){\displaystyle f(x;\alpha ,\beta ,c,\mu )}يرمز إلى كثافة X و Y هو مجموع النسخ المستقلة من X : Y=أنا=1شمالكأنا(Xأنا-μ){\displaystyle Y=\sum _{i=1}^{N}k_{i}(X_{i}-\mu )} إذن، فإن Y لها الكثافة1sو(y/s؛α،β،ج،0){\displaystyle {\tfrac {1}{s}}f(y/s;\alpha ,\beta ,c,0)}مع s=(أنا=1شمال|كأنا|α)1α{\displaystyle s=\left(\sum _{i=1}^{N}|k_{i}|^{\alpha }\right)^{\frac {1}{\alpha }}}

يتم وصف السلوك التقاربي، لـα<2{\displaystyle \alpha <2}بواسطة: [ 7 ]و(x)1|x|1+α(جα(1+علامة(x)β)الخطيئة(πα2)Γ(α+1)π){\displaystyle f(x)\sim {\frac {1}{|x|^{1+\alpha }}}\left(c^{\alpha }(1+\operatorname {sgn}(x)\beta )\sin \left({\frac {\pi \alpha }{2}}\right){\frac {\Gamma (\alpha +1)}{\pi }}\right)} حيث Γ هي دالة غاما (باستثناء عندماα1{\displaystyle \alpha \geq 1}وβ=±1{\displaystyle \beta =\pm 1}لا يختفي الذيل إلى يسار أو يمين μ ، على التوالي، على الرغم من أن التعبير أعلاه يساوي صفرًا. يتسبب سلوك " الذيل الثقيل " هذا في أن يكون تباين التوزيعات المستقرة لانهائيًا لجميعα<2{\displaystyle \alpha <2}يتم توضيح هذه الخاصية في الرسوم البيانية اللوغاريتمية أدناه.

متىα=2{\displaystyle \alpha =2}، التوزيع هو توزيع غاوسي (انظر أدناه)، مع ذيول مقاربة لـ exp(− x 2 /4 c 2 )/(2 c π ).

ملكيات

تكون التوزيعات المستقرة مغلقة تحت عملية الالتفاف لقيمة ثابتة لـα{\displaystyle \alpha }بما أن الالتفاف يكافئ ضرب الدالة المحولة باستخدام تحويل فورييه، فإنه يترتب على ذلك أن حاصل ضرب دالتين مميزتين مستقرتين لهما نفسα{\displaystyle \alpha }سينتج عن ذلك دالة مميزة أخرى مماثلة. ويُعطى حاصل ضرب دالتين مميزتين مستقرتين بالصيغة التالية: خبرة[أنات(μ1+μ2)-|ج1ت|α-|ج2ت|α+أنا(β1|ج1ت|α+β2|ج2ت|α)علامة(ت)Φ]{\displaystyle \exp \left[it\left(\mu _{1}+\mu _{2}\right)-|c_{1}t|^{\alpha }-|c_{2}t|^{\alpha }+i\left(\beta _{1}|c_{1}t|^{\alpha }+\beta _{2}|c_{2}t|^{\alpha }\right)\operatorname {sgn}(t)\Phi \right]}

بما أن Φ ليست دالة لـ μ أو c أوβ{\displaystyle \beta }وبناءً على المتغيرات، فإن هذه المعاملات للدالة الملتفة تُعطى بالصيغة التالية: μ=μ1+μ2ج=(ج1α+ج2α)1αβ=β1ج1α+β2ج2αج1α+ج2α{\displaystyle {\begin{aligned}\mu &=\mu _{1}+\mu _{2}\\c&=\left(c_{1}^{\alpha }+c_{2}^{\alpha }\right)^{\frac {1}{\alpha }}\\[6pt]\beta &={\frac {\beta _{1}c_{1}^{\alpha }+\beta _{2}c_{2}^{\alpha }}{c_{1}^{\alpha }+c_{2}^{\alpha }}}\end{aligned}}}

في كل حالة، يمكن إثبات أن المعلمات الناتجة تقع ضمن الفترات المطلوبة لتوزيع مستقر.

نظرية النهاية المركزية المعممة

كانت نظرية النهاية المركزية المعممة (GCLT) ثمرة جهد مشترك بين العديد من علماء الرياضيات ( بيرنشتاين ، ليندبيرغ ، ليفي ، فيلر ، كولموغوروف ، وغيرهم) خلال الفترة من عام 1920 إلى عام 1937. [ 12 ] نُشر أول برهان كامل (باللغة الفرنسية) لنظرية النهاية المركزية المعممة عام 1937 على يد بول ليفي . [ 13 ] وتتوفر نسخة باللغة الإنجليزية من البرهان الكامل لنظرية النهاية المركزية المعممة في ترجمة كتاب غنيدينكو وكولموغوروف الصادر عام 1954. [ 14 ]

وجاء في بيان المحكمة العامة لقانون الشركات ما يلي: [ 10 ]

نظرية النهاية المركزية المعممة — يكون المتغير العشوائي غير المنحط Z مستقرًا من النوع α لبعض القيم 0 < α ≤ 2 إذا وفقط إذا وُجدت متتالية مستقلة ومتطابقة التوزيع من المتغيرات العشوائية X1 ، X2 ، X3 ، ... وثوابت an > 0 ، bn مع

a n ( X 1 + ... + X n ) − b nZ.

هنا → تعني أن سلسلة مجاميع المتغيرات العشوائية تتقارب في التوزيع؛ أي أن التوزيعات المقابلة تحقق F n ( y ) → F ( y ) عند جميع نقاط استمرارية F.

بمعنى آخر، إذا تقاربت مجاميع المتغيرات العشوائية المستقلة والمتطابقة التوزيع في التوزيع إلى قيمة Z معينة ، فإن Z يجب أن يكون توزيعًا مستقرًا.

حالات خاصة

رسم بياني لوغاريتمي لتوزيعات الاحتمالية المتماثلة المركزية المستقرة، يُظهر سلوك قانون القوة لقيم x الكبيرة . ويتضح سلوك قانون القوة من خلال ظهور خط مستقيم في دالة كثافة الاحتمالية لقيم x الكبيرة ، حيث يساوي الميل 1.-(α+1){\displaystyle -(\alpha +1)}(الاستثناء الوحيد هو لـα=2{\displaystyle \alpha =2}(باللون الأسود، وهو توزيع طبيعي.)
رسم بياني لوغاريتمي لدوال كثافة الاحتمال لتوزيعات مستقرة مركزية ملتوية، يوضح سلوك قانون القوة لقيم x الكبيرة . مرة أخرى، ميل الأجزاء الخطية يساوي-(α+1){\displaystyle -(\alpha +1)}

لا يوجد حل تحليلي عام لشكل الدالة f ( x ). ومع ذلك، توجد ثلاث حالات خاصة يمكن التعبير عنها بدلالة الدوال الأولية ، كما يتضح من فحص الدالة المميزة : [ 7 ] [ 9 ] [ 15 ]

  • لα=2{\displaystyle \alpha =2}يتحول التوزيع إلى توزيع غاوسي بتباين σ² = 2c² ومتوسط ​​μ ؛ معامل الالتواءβ{\displaystyle \beta }ليس له أي تأثير.
  • لα=1{\displaystyle \alpha =1}وβ=0{\displaystyle \beta =0}يتحول التوزيع إلى توزيع كوشي بمعامل مقياس c ومعامل إزاحة μ .
  • لα=1/2{\displaystyle \alpha =1/2}وβ=1{\displaystyle \beta =1}يتحول التوزيع إلى توزيع ليفي بمعامل مقياس c ومعامل إزاحة μ .

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

صيغة عامة مغلقة لدوال كثافة الاحتمال المستقرة ذات القيم النسبية لـα{\displaystyle \alpha }تتوفر هذه الدوال بدلالة دوال ماير G. [ 17 ] كما يمكن استخدام دوال فوكس H للتعبير عن دوال كثافة الاحتمال المستقرة. بالنسبة للأعداد النسبية البسيطة، غالبًا ما يكون التعبير المغلق بدلالة دوال خاصة أقل تعقيدًا . تتوفر العديد من التعبيرات المغلقة ذات الصيغ البسيطة نسبيًا بدلالة الدوال الخاصة. في الجدول أدناه، يُشار إلى دوال كثافة الاحتمال التي يمكن التعبير عنها باستخدام الدوال الأولية بالحرف E ، بينما يُشار إلى تلك التي يمكن التعبير عنها باستخدام الدوال الخاصة بالحرف s . [ 16 ]

α{\displaystyle \alpha }
1/3نصف2/314/33/22
β{\displaystyle \beta }0sssهـssهـ
1sهـsلs

بعض الحالات الخاصة تُعرف بأسماء محددة:

  • لα=1{\displaystyle \alpha =1}وβ=1{\displaystyle \beta =1}، التوزيع هو توزيع لاندو ( L ) الذي له استخدام محدد في الفيزياء تحت هذا الاسم.
  • لα=3/2{\displaystyle \alpha =3/2}وβ=0{\displaystyle \beta =0}يختزل التوزيع إلى توزيع هولتسمارك بمعامل مقياس c ومعامل إزاحة μ .

أيضًا، في النهاية عندما يقترب c من الصفر أو عندما يقترب α من الصفر، سيقترب التوزيع من دالة ديراك دلتا δ ( x μ )  .

تمثيل المسلسل

يمكن إعادة صياغة التوزيع المستقر كجزء حقيقي من تكامل أبسط: [ 18 ]و(x؛α،β،ج،μ)=1π[0هـأنات(x-μ)هـ-(جت)α(1-أناβΦ)دت].{\displaystyle f(x;\alpha ,\beta ,c,\mu )={\frac {1}{\pi }}\Re \left[\int _{0}^{\infty }e^{it(x-\mu )}e^{-(ct)^{\alpha }(1-i\beta \Phi )}\,dt\right].}

وبتعبير الدالة الأسية الثانية على شكل متسلسلة تايلور ، ينتج عن ذلك ما يلي: و(x؛α،β،ج،μ)=1π[0هـأنات(x-μ)ن=0(-qتα)نن!دت]{\displaystyle f(x;\alpha ,\beta ,c,\mu )={\frac {1}{\pi }}\Re \left[\int _{0}^{\infty }e^{it(x-\mu )}\sum _{n=0}^{\infty }{\frac {(-qt^{\alpha })^{n}}{n!}}\,dt\right]} أينq=جα(1-أناβΦ){\displaystyle q=c^{\alpha }(1-i\beta \Phi )}. بعكس ترتيب التكامل والجمع، وإجراء التكامل ينتج ما يلي: و(x؛α،β،ج،μ)=1π[ن=1(-q)نن!(أناx-μ)αن+1Γ(αن+1)]{\displaystyle f(x;\alpha ,\beta ,c,\mu )={\frac {1}{\pi }}\Re \left[\sum _{n=1}^{\infty }{\frac {(-q)^{n}}{n!}}\left({\frac {i}{x-\mu }}\right)^{\alpha n+1}\Gamma (\alpha n+1)\right]} وهذا سيكون صالحًا عندما يكون x μ ، وسيتقارب لقيم مناسبة للمعاملات. (لاحظ أنه تم حذف الحد n = 0 الذي ينتج عنه دالة دلتا في xμ ). إن التعبير عن الدالة الأسية الأولى كمتسلسلة سينتج عنه متسلسلة أخرى في قوى موجبة لـ xμ ، وهي أقل فائدة بشكل عام.       

بالنسبة للتوزيع المستقر أحادي الجانب، يلزم تعديل متسلسلة التوسيع المذكورة أعلاه، لأنq=خبرة(-أناαπ/2){\displaystyle q=\exp(-i\alpha \pi /2)}وqأناα=1{\displaystyle qi^{\alpha }=1}لا يوجد جزء حقيقي للجمع. بدلاً من ذلك، يجب إجراء تكامل الدالة المميزة على المحور السالب، مما ينتج عنه: [ 19 ] [ 20 ]لα(x)=1π[ن=1(-q)نن!(-أناx)αن+1Γ(αن+1)]=1πن=1-الخطيئة(ن(α+1)π)ن!(1x)αن+1Γ(αن+1){\displaystyle {\begin{aligned}L_{\alpha }(x)&={\frac {1}{\pi }}\Re \left[\sum _{n=1}^{\infty }{\frac {(-q)^{n}}{n!}}\left({\frac {-i}{x}}\right)^{\alpha n+1}\Gamma (\alpha n+1)\right]\\[1ex]&={\frac {1}{\pi }}\sum _{n=1}^{\infty }{\frac {-\sin(n(\alpha +1)\pi )}{n!}}\left({\frac {1}{x}}\right)^{\alpha n+1}\Gamma (\alpha n+1)\end{aligned}}}

تقدير المعلمات

بالإضافة إلى الاختبارات الحالية للتوزيع الطبيعي وتقدير المعلمات اللاحق ، طوّر ماكولوتش طريقة عامة تعتمد على الكميات المئوية، وهي تعمل مع كل من التوزيعات المتماثلة والملتوية المستقرة ومعامل الاستقرار.0.5<α2{\displaystyle 0.5<\alpha \leq 2}[ 21 ]

محاكاة المتغيرات المستقرة

لا توجد تعابير تحليلية للمعكوسF-1(x){\displaystyle F^{-1}(x)}ولا صندوق الدفاع الكنديF(x){\displaystyle F(x)}لذا، لا يمكن استخدام طريقة الانعكاس لتوليد متغيرات ذات توزيع مستقر. [ 11 ] تتطلب الطرق القياسية الأخرى، مثل طريقة الرفض، حسابات مطولة. وقد اقترح تشامبرز ومالوز وستوك (CMS) حلاً أنيقًا وفعالًا، [ 22 ] حيث لاحظوا أن صيغة تكاملية معينة [ 23 ] تُنتج الخوارزمية التالية: [ 24 ]

  • توليد متغير عشوائييو{\displaystyle U}موزعة بشكل متساوٍ على(-π2،π2){\displaystyle \left(-{\tfrac {\pi }{2}},{\tfrac {\pi }{2}}\right)}ومتغير عشوائي أسي مستقلدبليو{\displaystyle W}بمتوسط ​​1؛
  • لα1{\displaystyle \alpha \neq 1}حساب:X=(1+ζ2)12αالخطيئة(α(يو+ξ))(كوس(يو))1α(كوس(يو-α(يو+ξ))دبليو)1-αα،{\displaystyle X=\left(1+\zeta ^{2}\right)^{\frac {1}{2\alpha }}{\frac {\sin(\alpha (U+\xi ))}{(\cos(U))^{\frac {1}{\alpha }}}}\left({\frac {\cos(U-\alpha (U+\xi ))}{W}}\right)^{\frac {1-\alpha }{\alpha }},}
  • لα=1{\displaystyle \alpha =1}حساب:X=1ξ{(π2+βيو)لون برونزييو-βسجل(π2دبليوكوسيوπ2+βيو)}،{\displaystyle X={\frac {1}{\xi }}\left\{\left({\frac {\pi }{2}}+\beta U\right)\tan U-\beta \log \left({\frac {{\frac {\pi }{2}}W\cos U}{{\frac {\pi }{2}}+\beta U}}\right)\right\},}أينζ=-βلون برونزيπα2،ξ={1αدالة الظل العكسي(-ζ)α1π2α=1{\displaystyle \zeta =-\beta \tan {\frac {\pi \alpha }{2}},\qquad \xi ={\begin{cases}{\frac {1}{\alpha }}\arctan(-\zeta )&\alpha \neq 1\\{\frac {\pi }{2}}&\alpha =1\end{cases}}}

تُنتج هذه الخوارزمية متغيرًا عشوائيًاXSα(β،1،0){\displaystyle X\sim S_{\alpha }(\beta ,1,0)}للاطلاع على برهان مفصل، انظر [ 25 ] .

لمحاكاة متغير عشوائي مستقر لجميع القيم المسموح بها للمعاملاتα{\displaystyle \alpha }،ج{\displaystyle c}،β{\displaystyle \beta }وμ{\displaystyle \mu }استخدم الخاصية التالية: إذاXSα(β،1،0){\displaystyle X\sim S_{\alpha }(\beta ,1,0)}ثم Y={جX+μα1جX+μ+2πβجسجلجα=1{\displaystyle Y={\begin{cases}cX+\mu &\alpha \neq 1\\cX+\mu +{\frac {2}{\pi }}\beta c\log c&\alpha =1\end{cases}}} يكونSα(β،ج،μ){\displaystyle S_{\alpha }(\beta ,c,\mu )}. لα=2{\displaystyle \alpha =2}β=0{\displaystyle \beta =0}تُختزل طريقة CMS إلى تحويل بوكس-مولر المعروف لتوليد متغيرات عشوائية غاوسية . [ 26 ] في حين تم اقتراح مناهج أخرى في الأدبيات، بما في ذلك تطبيق متسلسلات بيرجستروم [ 27 ] وليباج [ 28 ] ، تُعتبر طريقة CMS الأسرع والأكثر دقة.

التطبيقات

تستمد التوزيعات المستقرة أهميتها، نظريًا وعمليًا، من تعميم نظرية النهاية المركزية على المتغيرات العشوائية التي لا تمتلك عزمًا من الرتبة الثانية (وربما الأولى)، وما يصاحب ذلك من تشابه ذاتي للعائلة المستقرة. وقد دفع الانحراف الظاهر عن التوزيع الطبيعي، إلى جانب الحاجة إلى نموذج متشابه ذاتيًا للبيانات المالية (أي أن شكل توزيع تغيرات أسعار الأصول السنوية يجب أن يشبه شكل توزيع تغيرات الأسعار اليومية أو الشهرية المكونة له)، بينوا ماندلبروت إلى اقتراح أن أسعار القطن تتبع توزيعًا مستقرًا من نوع ألفا.α{\displaystyle \alpha }يساوي 1.7. [ 6 ] تُستخدم توزيعات ليفي بشكل متكرر في تحليل السلوك الحرج والبيانات المالية. [ 9 ] [ 29 ]

كما أنها موجودة في التحليل الطيفي كتعبير عام عن خط طيفي متسع بفعل الضغط شبه الساكن . [ 18 ]

تم توضيح توزيع ليفي لفترات انتظار التوهجات الشمسية (الفترة الزمنية بين أحداث التوهج) لتوهجات الأشعة السينية الصلبة التي رصدها مرصد CGRO BATSE في ديسمبر 2001. وكشف تحليل البصمة الإحصائية لليفي عن وجود بصمتين مختلفتين للذاكرة؛ إحداهما مرتبطة بالدورة الشمسية، والأخرى يبدو أن أصلها مرتبط بتأثيرات موضعية أو مجموعة من التأثيرات الموضعية للمناطق النشطة الشمسية. [ 30 ]

حالات تحليلية أخرى

هناك عدد من حالات التوزيعات المستقرة القابلة للتعبير التحليلي. لنفترض أن التوزيع المستقر يُعبَّر عنه بالصيغة التالية:و(x؛α،β،ج،μ){\displaystyle f(x;\alpha ,\beta ,c,\mu )}، ثم:

  • يُعطى توزيع كوشي بالصيغة التاليةو(x؛1،0،1،0).{\displaystyle f(x;1,0,1,0).}
  • يُعطى توزيع ليفي بالصيغة التاليةو(x؛12،1،1،0).{\displaystyle f(x;{\tfrac {1}{2}},1,1,0).}
  • يُعطى التوزيع الطبيعي بالصيغة التالية :و(x؛2،0،1،0).{\displaystyle f(x;2,0,1,0).}
  • يتركSμ،ν(z){\displaystyle S_{\mu ,\nu }(z)}إذا كانت دالة لوميل ، فإن: [ 31 ]و(x؛13،0،1،0)=(2هـ-أناπ433π1x3S0،13(2هـأناπ4331x)){\displaystyle f{\left(x;{\tfrac {1}{3}},0,1,0\right)}=\Re \left({\frac {2e^{-{\frac {i\pi }{4}}}}{3{\sqrt {3}}\pi }}{\frac {1}{\sqrt {x^{3}}}}S_{0,{\frac {1}{3}}}{\left({\frac {2e^{\frac {i\pi }{4}}}{3{\sqrt {3}}}}{\frac {1}{\sqrt {x}}}\right)}\right)}
  • يتركS(x){\displaystyle S(x)}وج(x){\displaystyle C(x)}إذا رمزنا إلى تكاملات فرينل ، فإن: [ 32 ]و(x؛12،0،1،0)=(12π|x|3)1/2(الخطيئة(14|x|)[12-S(12π|x|)]+كوس(14|x|)[12-ج(12π|x|)]){\displaystyle f{\left(x;{\tfrac {1}{2}},0,1,0\right)}=\left({\tfrac {1}{2\pi \left|x\right|^{3}}}\right)^{1/2}\left(\sin \left({\tfrac {1}{4|x|}}\right)\left[{\tfrac {1}{2}}-S{\left({\tfrac {1}{\sqrt {2\pi |x|}}}\right)}\right]+\cos \left({\tfrac {1}{4|x|}}\right)\left[{\tfrac {1}{2}}-C{\left({\tfrac {1}{\sqrt {2\pi |x|}}}\right)}\right]\right)}
  • يترككv(x){\displaystyle K_{v}(x)}إذا كانت دالة بيسل المعدلة من النوع الثاني، فإن: [ 32 ]و(x؛13،1،1،0)=252374π1x3ك13(2523941x){\displaystyle f{\left(x;{\tfrac {1}{3}},1,1,0\right)}={\frac {2^{\frac {5}{2}}}{3^{\frac {7}{4}}\pi }}{\frac {1}{\sqrt {x^{3}}}}K_{\frac {1}{3}}{\left({\frac {2^{\frac {5}{2}}}{3^{\frac {9}{4}}}}{\frac {1}{\sqrt {x}}}\right)}}
  • يتركمFن{\displaystyle {}_{m}F_{n}}إذا رمزنا للدوال الهندسية الفائقة ، فإن: [ 31 ]و(x؛43،0،1،0)=354252π12Γ(712)Γ(1112)Γ(612)Γ(812)2F2(712،1112؛612،812؛33x444)-3114x32132π12Γ(1312)Γ(1712)Γ(1812)Γ(1512)2F2(1312،1712؛1812،1512؛33x444){\displaystyle {\begin{aligned}f{\left(x;{\tfrac {4}{3}},0,1,0\right)}&={\frac {3^{\frac {5}{4}}}{2^{\frac {5}{2}}\pi ^{\frac {1}{2}}}}{\frac {\Gamma {\left({\tfrac {7}{12}}\right)}\,\Gamma {\left({\tfrac {11}{12}}\right)}}{\Gamma {\left({\tfrac {6}{12}}\right)}\,\Gamma {\left({\tfrac {8}{12}}\right)}}}\;{}_{2}F_{2}{\left({\tfrac {7}{12}},{\tfrac {11}{12}};{\tfrac {6}{12}},{\tfrac {8}{12}};{\tfrac {3^{3}x^{4}}{4^{4}}}\right)}\\[2pt]&\quad -{\frac {3^{\frac {11}{4}}x^{3}}{2^{\frac {13}{2}}\pi ^{\frac {1}{2}}}}{\frac {\Gamma {\left({\tfrac {13}{12}}\right)}\,\Gamma {\left({\tfrac {17}{12}}\right)}}{\Gamma {\left({\tfrac {18}{12}}\right)}\,\Gamma {\left({\tfrac {15}{12}}\right)}}}\;{}_{2}F_{2}{\left({\tfrac {13}{12}},{\tfrac {17}{12}};{\tfrac {18}{12}},{\tfrac {15}{12}};{\tfrac {3^{3}x^{4}}{4^{4}}}\right)}\end{aligned}}}و(x؛32،0،1،0)=Γ(53)π2F3(512،1112؛13،12،56؛-22x636)-x23π3F4(34،1،54؛23،56،76،43؛-22x636)+7x4Γ(43)34π22F3(1312،1912؛76،32،53؛-22x636){\displaystyle {\begin{aligned}f{\left(x;{\tfrac {3}{2}},0,1,0\right)}&={\frac {\Gamma {\left({\tfrac {5}{3}}\right)}}{\pi }}{}_{2}F_{3}{\left({\tfrac {5}{12}},{\tfrac {11}{12}};{\tfrac {1}{3}},{\tfrac {1}{2}},{\tfrac {5}{6}};-{\tfrac {2^{2}x^{6}}{3^{6}}}\right)}\\[2pt]&\quad -{\frac {x^{2}}{3\pi }}\,{}_{3}F_{4}{\left({\tfrac {3}{4}},1,{\tfrac {5}{4}};{\tfrac {2}{3}},{\tfrac {5}{6}},{\tfrac {7}{6}},{\tfrac {4}{3}};-{\tfrac {2^{2}x^{6}}{3^{6}}}\right)}\\[2pt]&\quad +{\frac {7x^{4}\Gamma {\left({\tfrac {4}{3}}\right)}}{3^{4}\pi ^{2}}}{}_{2}F_{3}{\left({\tfrac {13}{12}},{\tfrac {19}{12}};{\tfrac {7}{6}},{\tfrac {3}{2}},{\tfrac {5}{3}};-{\tfrac {2^{2}x^{6}}{3^{6}}}\right)}\end{aligned}}}أما الأخيرة فهي شركة التوزيع هولتسمارك .
  • يتركدبليوك،μ(z){\displaystyle W_{k,\mu }(z)}إذا كانت دالة ويتاكر ، فإن: [ 33 ] [ 34 ] [ 35 ]و(x؛23،0،1،0)=36π|x|خبرة(227x-2)دبليو-12،16(427x-2)و(x؛23،1،1،0)=3π|x|خبرة(-1627x-2)دبليو12،16(3227x-2)و(x؛32،1،1،0)={3π|x|خبرة(127x3)دبليو12،16(-227x3)x<036π|x|خبرة(127x3)دبليو-12،16(227x3)x0{\displaystyle {\begin{aligned}f\left(x;{\tfrac {2}{3}},0,1,0\right)&={\frac {\sqrt {3}}{6{\sqrt {\pi }}|x|}}\exp \left({\tfrac {2}{27}}x^{-2}\right)W_{-{\frac {1}{2}},{\frac {1}{6}}}\left({\tfrac {4}{27}}x^{-2}\right)\\[8pt]f\left(x;{\tfrac {2}{3}},1,1,0\right)&={\frac {\sqrt {3}}{{\sqrt {\pi }}|x|}}\exp \left(-{\tfrac {16}{27}}x^{-2}\right)W_{{\frac {1}{2}},{\frac {1}{6}}}\left({\tfrac {32}{27}}x^{-2}\right)\\[8pt]f\left(x;{\tfrac {3}{2}},1,1,0\right)&={\begin{cases}{\frac {\sqrt {3}}{{\sqrt {\pi }}|x|}}\exp \left({\frac {1}{27}}x^{3}\right)W_{{\frac {1}{2}},{\frac {1}{6}}}\left(-{\frac {2}{27}}x^{3}\right)&x<0\\{}\\{\frac {\sqrt {3}}{6{\sqrt {\pi }}|x|}}\exp \left({\frac {1}{27}}x^{3}\right)W_{-{\frac {1}{2}},{\frac {1}{6}}}\left({\frac {2}{27}}x^{3}\right)&x\geq 0\end{cases}}\end{aligned}}}

انظر أيضاً

تطبيقات البرمجيات

  • يتوفر برنامج STABLE لنظام التشغيل Windows على صفحة الويب الخاصة بجون نولان: http://www.robustanalysis.com/public/stable.html . يحسب البرنامج دالة الكثافة (pdf) ودالة التوزيع التراكمي (cdf) والكميات لتوزيع مستقر عام، ويُجري تقدير الاحتمال الأقصى للمعلمات المستقرة، بالإضافة إلى بعض تقنيات تحليل البيانات الاستكشافية لتقييم مدى ملاءمة مجموعة البيانات.
  • تحتوي مكتبة GNU العلمية المكتوبة بلغة C على حزمة randist ، والتي تتضمن من بين توزيعات Gaussian و Cauchy أيضًا تطبيقًا لتوزيع Levy alpha-stable، سواء مع أو بدون معلمة الانحراف.
  • libstable هو تطبيق C لوظائف توزيع الاحتمال الثابت، ودالة التوزيع التراكمي، والرقم العشوائي، والكمية، والتوفيق (إلى جانب حزمة تكرار قياسية وحزمة R).
  • حزمة R المسماة 'stabledist' من تطوير ديتلم وورتز ومارتن ماخلر وأعضاء الفريق الأساسي لـ Rmetrics. تحسب هذه الحزمة الكثافة المستقرة والاحتمالية والكميات والأعداد العشوائية.
  • يوجد تطبيق بايثون في scipy.stats.levy_stable في حزمة SciPy .
  • توفر لغة جوليا حزمة StableDistributions.jl التي تتضمن طرقًا لتوليد التوزيعات المستقرة، وملاءمتها، وحساب كثافة الاحتمال، ودالة التوزيع التراكمي، ودوال توليد الخصائص والعزوم، ودوال الكميات وما يتصل بها، والتحويلات الالتفافية والخطية. وتستخدم هذه الحزمة خوارزميات مُحسّنة من قِبل جون ب. نولان. [ 10 ]

مراجع

  1. 1 2 ماندلبروت، ب. (1960). "قانون باريتو-ليفي وتوزيع الدخل". المجلة الاقتصادية الدولية . 1 (2): 79-106 . doi : 10.2307/2525289 . JSTOR 2525289 . 
  2. ^ ليفي ، بول (1925). حساب الاحتمالات . باريس: غوتييه فيلار. او سي ال سي 1417531 . 
  3. ماندلبروت، ب. (1961). "دوال باريتو العشوائية المستقرة والتغير المضاعف للدخل". إيكونومتريكا . 29 (4): 517-543 . doi : 10.2307/1911802 . JSTOR 1911802 . 
  4. ماندلبروت، ب. (1963). "تغير بعض أسعار المضاربة". مجلة الأعمال . 36 (4): 394-419 . doi : 10.1086/294632 . JSTOR 2350970 . 
  5. فاما، يوجين ف. (1963). "ماندلبروت وفرضية باريتو المستقرة". مجلة الأعمال . 36 (4): 420-429 . doi : 10.1086/294633 . JSTOR 2350971 . 
  6. 1 2 ماندلبروت، ب. (1963). "أساليب جديدة في الاقتصاد الإحصائي". مجلة الاقتصاد السياسي . 71 (5): 421-440 . doi : 10.1086/258792 . S2CID 53004476 . 
  7. 1 2 3 4 5 6 نولان، جون ب. "التوزيعات المستقرة - نماذج للبيانات ذات الذيول السميكة" (ملف PDF) . مؤرشف من الأصل (ملف PDF) بتاريخ 17 يوليو 2011. تم الاطلاع عليه بتاريخ 21 فبراير 2009 .
  8. سيغريست، كايل. "التوزيعات المستقرة" . www.randomservices.org . تم الاطلاع عليه بتاريخ 18-10-2018 .
  9. 1 2 3 فويت، يوهانس (2005). باليان، ر؛ بيغلبوك، و؛ غروس، هـ؛ ثيرينغ، و (محررون). الميكانيكا الإحصائية للأسواق المالية - سبرينغر . نصوص ودراسات في الفيزياء. سبرينغر. doi : 10.1007/b137351 . ISBN 978-3-540-26285-5.
  10. 1 2 3 نولان، جون ب. (2020). التوزيعات المستقرة أحادية المتغير، نماذج للبيانات ذات الذيول الثقيلة . سلسلة سبرينغر في بحوث العمليات والهندسة المالية. سويسرا: سبرينغر. doi : 10.1007/978-3-030-52915-4 . ISBN 978-3-030-52914-7. S2CID 226648987 . 
  11. 1 2 نولان، جون ب. (1997). "الحساب العددي للكثافات المستقرة ودوال التوزيع". الاتصالات في الإحصاء. النماذج العشوائية . 13 (4): 759-774 . doi : 10.1080/15326349708807450 . ISSN 0882-0287 . 
  12. لو كام، ل. (فبراير 1986). "نظرية النهاية المركزية حوالي عام 1935". العلوم الإحصائية . 1 (1): 78-91 . JSTOR 2245503 . 
  13. ^ ليفي ، بول (1937). Theorie de l'addition des Variables aleatoires [نظرية الجمع بين المتغيرات غير المتوقعة] . باريس: غوتييه فيلار.
  14. غنيدينكو، بوريس فلاديميروفيتش؛ كولوغوروف، أندريه نيكولايفيتش؛ دوب، جوزيف ل.؛ هسو، باو-لو (1968). التوزيعات الحدية لمجاميع المتغيرات العشوائية المستقلة . ريدينغ، ماساتشوستس: أديسون-ويسلي.
  15. سامورودنيتسكي، ج.؛ تاكو، م.س. (1994). العمليات العشوائية غير الغاوسية المستقرة: نماذج عشوائية ذات تباين لانهائي . مطبعة سي آر سي. رقم ISBN 978-0-412-05171-5.
  16. 1 2 لي، واي ها (2010). الخصائص المستمرة والمتقطعة للعمليات العشوائية . أطروحة دكتوراه، جامعة نوتنغهام.
  17. زولوتاريف، ف. (1995). "حول تمثيل كثافات القوانين المستقرة بواسطة دوال خاصة". نظرية الاحتمالات وتطبيقاتها . 39 (2): 354-362 . doi : 10.1137/1139025 . ISSN 0040-585X . 
  18. 1 2 بيتش، ج. (1981). "نظرية اتساع وانزياح خطوط الطيف بفعل الضغط". التقدم في الفيزياء . 30 (3): 367-474 . Bibcode : 1981AdPhy..30..367P . doi : 10.1080/00018738100101467 . ISSN 0001-8732 . 
  19. بولارد، هوارد (1946). "تمثيل e^{-x^\lambda} كتكامل لابلاس" . نشرة الجمعية الأمريكية للرياضيات 52 : 908. doi : 10.1090/S0002-9904-1946-08672-3 .
  20. بينسون، ك.أ.؛ غورسكا، ك. (17-11-2010). "كثافات الاحتمال الدقيقة والصريحة لتوزيعات ليفي المستقرة أحادية الجانب". رسائل المراجعة الفيزيائية . 105 (21) 210604. arXiv : 1007.0193 . Bibcode : 2010PhRvL.105u0604P . doi : 10.1103/PhysRevLett.105.210604 . PMID 21231282. S2CID 27497684 .  
  21. ماكولوتش، ج. هوستون (1986). "مُقدِّرات بسيطة ومتسقة لمعلمات التوزيع المستقر" (ملف PDF) . الاتصالات في الإحصاء. المحاكاة والحساب . 15 (4): 1109-1136 . doi : 10.1080/03610918608812563 .
  22. تشامبرز، جيه إم؛ مالوز، سي إل؛ ستوك، بي دبليو (1976). "طريقة لمحاكاة المتغيرات العشوائية المستقرة". مجلة الجمعية الإحصائية الأمريكية . 71 (354): 340-344 . doi : 10.1080/01621459.1976.10480344 . ISSN 0162-1459 . 
  23. زولوتاريف، ف.م. (1986). التوزيعات المستقرة أحادية البعد . الجمعية الأمريكية للرياضيات. ISBN 978-0-8218-4519-6.
  24. ميسيوريك، آدم؛ ويرون، رافال (2012). جنتل، جيمس إي؛ هاردل، وولفغانغ كارل؛ موري، يويتشي (محررون). التوزيعات ذات الذيول الثقيلة في حسابات القيمة المعرضة للخطر (ملف PDF) . كتيبات سبرينغر للإحصاءات الحاسوبية. سبرينغر برلين هايدلبرغ. الصفحات 1025-1059 . doi : 10.1007/978-3-642-21551-3_34 . ISBN  978-3-642-21550-6.
  25. فيرون، رافال (1996). "حول طريقة تشامبرز-مالوز-ستوك لمحاكاة المتغيرات العشوائية المستقرة المنحرفة". رسائل الإحصاء والاحتمالات . 28 (2): 165-171 . CiteSeerX 10.1.1.46.3280 . doi : 10.1016/0167-7152(95)00113-1 . S2CID 9500064 .  
  26. يانيكي، ألكسندر؛ ويرون، ألكسندر (1994). محاكاة وسلوك فوضوي للعمليات العشوائية المستقرة ألفا . مطبعة سي آر سي. رقم ISBN 978-0-8247-8882-7.
  27. مانتينيا، روزاريو نونزيو (1994). "خوارزمية سريعة ودقيقة للمحاكاة العددية للعمليات العشوائية المستقرة وفقًا لقانون ليفي". مجلة Physical Review E. 49 ( 5): 4677–4683 . Bibcode : 1994PhRvE..49.4677M . doi : 10.1103/PhysRevE.49.4677 . PMID 9961762 . 
  28. يانيكي، ألكسندر؛ كوكوشكا، بيوتر (1992). "دراسة حاسوبية لمعدل تقارب متسلسلات ليباج إلى متغيرات عشوائية مستقرة من النوع ألفا". الإحصاء . 23 (4): 365-373 . doi : 10.1080/02331889208802383 . ISSN 0233-1888 . 
  29. راتشيف، سفيتلوزار ت.؛ ميتنيك، ستيفان (2000). نماذج باريتو المستقرة في التمويل . وايلي. ISBN 978-0-471-95314-2.
  30. ليدون، د.، دراسة إحصائية للتوهجات الشمسية ذات الأشعة السينية الصلبة
  31. 1 2 غاروني، تي إم؛ فرانكل، إن إي (2002). "رحلات ليفي: نتائج دقيقة وسلوك تقاربي يتجاوز جميع الرتب". مجلة الفيزياء الرياضية . 43 (5): 2670-2689 . Bibcode : 2002JMP....43.2670G . doi : 10.1063/1.1467095 .
  32. 1 2 هوبكرافت، كي آي؛ جاكمان، إي؛ تانر، آر إم جيه (1999). "المسارات العشوائية لليفي ذات عدد الخطوات المتذبذب والسلوك متعدد المقاييس". مجلة Physical Review E. 60 ( 5): 5327-5343 . Bibcode : 1999PhRvE..60.5327H . doi : 10.1103/physreve.60.5327 . PMID 11970402 . 
  33. أوتشايكين، في في؛ زولوتاريف، في إم (1999). "الصدفة والاستقرار - التوزيعات المستقرة وتطبيقاتها". VSP .
  34. زلوتاريف، ف.م. (1961). "التعبير عن كثافة التوزيع المستقر ذي الأس ألفا الأكبر من واحد باستخدام التردد ذي الأس 1/ألفا". مختارات من الترجمات في الإحصاء الرياضي والاحتمالات (مترجم من المقالة الروسية: دوكل. أكاد. ناوك إس إس إس آر. 98، 735-738 (1954)) . 1 : 163-167 .
  35. زاليابين، الرابع؛ كاغان، واي واي؛ شونبيرغ، إف بي (2005). "تقريب توزيع مجاميع باريتو" . الجيوفيزياء البحتة والتطبيقية . 162 (6): 1187-1228 . Bibcode : 2005PApGe.162.1187Z . doi : 10.1007/s00024-004-2666-3 . S2CID 18754585 .