התפלגות פארטו

</math>
| מצטברת = <math>1-\left(\frac{x_\mathrm{m}}{x}\right)^\alpha</math>
| תוחלת = <math>\begin{cases}
 \infty & \text{for }\alpha\le 1 \\
 \dfrac{\alpha x_\mathrm{m}}{\alpha-1} & \text{for }\alpha>1
 \end{cases}</math>
| סטיית תקן =
| חציון = <math>x_\mathrm{m} \sqrt[\alpha]{2}</math>
| שכיח = <math>x_\mathrm{m}</math>
| שונות = <math>\begin{cases}
 \infty & \text{for }\alpha\le 2 \\
 \dfrac{x_\mathrm{m}^2\alpha}{(\alpha- 1)^2(\alpha-2)} & \text{for }\alpha>2
 \end{cases}</math>
| אנטרופיה = <math>\log\left(\left(\frac{x_\mathrm{m}}{\alpha}\right)\,e^{1+\tfrac{1}{\alpha}}\right) </math>
| מומנטים =
| אופיינית =
| צידוד = <math>\frac{2(1+\alpha)}{\alpha-3}\sqrt{\frac{\alpha-2}{\alpha}}\text{ for }\alpha>3</math>
| גבנוניות = <math>\frac{6(\alpha^3+\alpha^2-6\alpha-2)}{\alpha(\alpha-3)(\alpha-4)}\text{ for }\alpha>4</math>
}}
התפלגות פארטו, הנקראת על שם המהנדס האזרחי, הכלכלן והסוציולוג האיטלקי וילפרדו פארטו,<ref></ref> היא התפלגות הסתברות המשמשת לתיאור של תופעות מתחומים מגוונים: חברתיות, בקרת איכות, מדעיות, גיאופיזיות, אקטואריות ועוד. העיקרון, שיושם במקור לתיאור של חלוקת העושר בחברה, מתאים לתופעה לפיה חלק גדול מהעושר מוחזק על ידי חלק קטן מהאוכלוסייה.<ref></ref><ref name=":1">Pareto, Vilfredo, Cours d'Économie Politique: Nouvelle édition par G.-H. </ref> עקרון פארטו או "כלל 80-20", הקובע ש-80% מהתוצאות נובעות מ-20% מהסיבות, נקרא על שם פארטו, אך רק התפלגויות פארטו עם פרמטר צורה שערכו <math>\log_4 5 \approx 1.16</math> ישקפו במדויק את עיקרון פארטו. תצפיות אמפיריות הראו שהתפלגות 80-20 זו מתאימה למגוון רחב של תופעות טבע<ref></ref> והתנהגות אנושיות.<ref></ref><ref></ref>

== הגדרות ==
אם X הוא משתנה מקרי עם התפלגות פארטו (סוג I),<ref name="arnold"></ref> אזי ההסתברות ש- X גדול ממספר כלשהו x, כלומר, פונקציית ההישרדות (נקראת גם פונקציית זנב), נתונה על ידי

: <math>\overline{F}(x) = \Pr(X>x) = \begin{cases}
\left(\frac{x_\mathrm{m}}{x}\right)^\alpha & x\ge x_\mathrm{m}, \\
1 & x < x_\mathrm{m},
\end{cases}
</math>

כאשר x<sub>m</sub> הוא הערך המינימלי האפשרי (חיובי בהכרח) של X, ו- α הוא פרמטר חיובי. התפלגות פארטו מסוג I מאופיינת בפרמטר קנה-המידה x<sub>m</sub> ובפרמטר צורה α, הידוע כמדד הזנב. אם התפלגות זו משמשת למודל של התפלגות העושר, אז הפרמטר α נקרא אינדקס פארטו.

=== פונקציית התפלגות מצטברת ===
מההגדרה, פונקציית ההתפלגות המצטברת של משתנה מקרי שמתפלג פארטו עם הפרמטרים α ו- x<sub>m</sub> היא

: <math>F_X(x) = \begin{cases}
1-\left(\frac{x_\mathrm{m}}{x}\right)^\alpha & x \ge x_\mathrm{m}, \\
0 & x < x_\mathrm{m}.
\end{cases}</math>

=== פונקציית צפיפות הסתברות ===
על ידי גזירה מקבלים שפונקציית צפיפות ההסתברות היא

: <math>f_X(x)= \begin{cases} \frac{\alpha x_\mathrm{m}^\alpha}{x^{\alpha+1}} & x \ge x_\mathrm{m}, \\ 0 & x < x_\mathrm{m}. \end{cases} </math>

כאשר מציירים את ההתפלגות על צירים ליניאריים, ההתפלגות מתוארת על ידי עקומה בצורת האות J, ששואפת אסימפטוטית לשני הצירים. כל מקטעי העקומה דומים לעצמם (עם פרמטרי קנה מידה מתאימים). כאשר משרטטים בדיאגרמה שבה שני הצירים מכויילים לוגריתמית, ההתפלגות מיוצגת על ידי קו ישר.

== מאפיינים ==
=== מומנטים ופונקציה אופיינית ===
* התוחלת של משתנה מקרי שמתפלג פארטו היא

:: <math>\operatorname{E}(X)= \begin{cases} \infty & \alpha\le 1, \\
\frac{\alpha x_{\mathrm{m}}}{\alpha-1} & \alpha>1.
\end{cases}</math>
* השונות של משתנה מקרי שמתפלג פארטו היא

:: <math>\operatorname{Var}(X)= \begin{cases}
\infty & \alpha\in(1,2], \\
\left(\frac{x_\mathrm{m}}{\alpha-1}\right)^2 \frac{\alpha}{\alpha-2} & \alpha>2.
\end{cases}</math>

: (עבור α ≤ 1, השונות אינה מוגדרת)
* המומנטים הגולמיים הם

:: <math>\mu_n'= \begin{cases} \infty & \alpha\le n, \\ \frac{\alpha x_\mathrm{m}^n}{\alpha-n} & \alpha>n. \end{cases}</math>
* הפונקציה היוצרת מוגדרת רק עבור ערכים שליליים או אפס <math>t \leq 0</math>

:: <math>M\left(t;\alpha,x_\mathrm{m}\right) = \operatorname{E} \left [e^{tX} \right ] = \alpha(-x_\mathrm{m} t)^\alpha\Gamma(-\alpha,-x_\mathrm{m} t)</math>
:: <math>M\left(0,\alpha,x_\mathrm{m}\right)=1.</math>

היות שהתוחלת אינה מתכנסת בקטע הפתוח המכיל את <math>t=0</math> אנו אומרים שהפונקציה היוצרת אינה מוגדרת.
* הפונקציה האופיינית נתונה על ידי

:: <math>\varphi(t;\alpha,x_\mathrm{m})=\alpha(-ix_\mathrm{m} t)^\alpha\Gamma(-\alpha,-ix_\mathrm{m} t),</math>

: כאשר<math>\Gamma(a, x)</math> היא פונקציית גמא הלא שלמה.

ניתן לחלץ את הפרמטרים באמצעות שיטת המומנטים.<ref>S. Hussain, S.H. Bhatti (2018). </ref>

=== התפלגויות מותנות ===
התפלגות ההסתברות המותנית של משתנה מקרי המתפלג פארטו, בהינתן המאורע שערכו של המשתנה המקרי גדול או שווה למספר מסוים <math>x_1</math>, שגדול מ <math>x_\text{m}</math>, היא התפלגות פארטו עם אותו אינדקס פארטו <math>\alpha</math> אך עם מינימום <math>x_1</math> בִּמקוֹם <math>x_\text{m}</math>:

: <math>
\text{Pr}(X \geq x | X \geq x_1) =
\begin{cases}
\left(\frac{x_1}{x}\right)^\alpha & x \geq x_1, \\
1 & x < x_1.
\end{cases}
</math>

הביטוי מרמז כל כך שהערך התוחלת המותנית (אם היא סופית, כלומר <math>\alpha>1</math> ) פרופורציוני ל-<math>x_1</math> :

: <math>\text{E}(X | X \geq x_1) \propto x_1 </math>

במקרה של משתנים אקראיים המתארים את משך חייו של אובייקט, המשמעות היא שתוחלת החיים היא פרופורציונלית לגיל, והיא נקראת אפקט לינדי או חוק לינדי.<ref name=":02"></ref>

=== משפט אפיון ===
נניח כי <math>X_1, X_2, X_3, \dotsc</math> הם משתנים מקריים בלתי תלויים אשר התפלגות ההסתברות שלהם נתמכת בקטע <math>[x_\text{m},\infty)</math> עבור <math>x_\text{m}>0</math> כלשהו. נניח שלכל <math>n</math>, שני המשתנים האקראיים <math>\min\{X_1,\dotsc,X_n\}</math> ו <math>(X_1+\dotsb+X_n)/\min\{X_1,\dotsc,X_n\}</math> הם בלתי תלויים. אז ההתפלגות המשותפת היא התפלגות פארטו.

=== ממוצע גאומטרי ===
הממוצע הגאומטרי (G) הוא<ref name="Johnson1994">Johnson NL, Kotz S, Balakrishnan N (1994) Continuous univariate distributions Vol 1. </ref>

: <math> G = x_\text{m} \exp \left( \frac{1}{\alpha} \right).</math>

=== ממוצע הרמוני ===
הממוצע ההרמוני ( H ) הוא

: <math> H = x_\text{m} \left( 1 + \frac{ 1 }{ \alpha } \right).</math>

=== ייצוג גרפי ===
התפלגות ה'זנב הארוך' המעוקל האופיינית, כשהיא משורטטת בסקאלה ליניארית, אינה מגלה את הפשטות הבסיסית של הפונקציה. כשהיא משורטטת על מערכת צירים בסקאלה לוגריתמית-לוגריתמית, היא מתוארת כקו ישר עם שיפוע שלילי, כפי שניתן לראות מהנוסחה עבור פונקציית צפיפות ההסתברות עבור x ≥ x <sub>m</sub>

: <math>\log f_X(x)= \log \left(\alpha\frac{x_\mathrm{m}^\alpha}{x^{\alpha+1}}\right) = \log (\alpha x_\mathrm{m}^\alpha) - (\alpha+1) \log x.</math>

מכיוון ש- α חיובי, השיפוע <math>-(\alpha + 1)</math> הוא שלילי.

== הסקה סטטיסטית ==
=== הערכת פרמטרים ===
פונקציית הנראות של פרמטרי התפלגות פארטו α ו- x<sub>m</sub>, בהינתן מדגם בלתי תלוי x =&nbsp;(x<sub>1</sub> ,&nbsp;x<sub>2</sub> ,&nbsp;...,&nbsp;x<sub>n</sub>), הוא

: <math>L(\alpha, x_\mathrm{m}) = \prod_{i=1}^n \alpha \frac {x_\mathrm{m}^\alpha} {x_i^{\alpha+1}} = \alpha^n x_\mathrm{m}^{n\alpha} \prod_{i=1}^n \frac {1}{x_i^{\alpha+1}}.</math>

לכן, הלוגריתם הטבעי של פונקציית הנראות הוא

: <math>\ell(\alpha, x_\mathrm{m}) = n \ln \alpha + n\alpha \ln x_\mathrm{m} - (\alpha + 1) \sum_{i=1} ^n \ln x_i.</math>

אפשר לראות כי <math>\ell(\alpha, x_\mathrm{m})</math> עולה מונוטונית עם x<sub>m</sub>, כלומר ככל שהערך של x<sub>m</sub> גדל, כך גדל הערך של פונקציית הנראות. לפיכך, מכיוון ש- x ≥ x<sub>m</sub>, אנו מסיקים כי

: <math>\widehat x_\mathrm{m} = \min_i {x_i}.</math>

כדי למצוא את המשערך של α, מחשבים את הנגזרת החלקית המתאימה ומוצאים את הערך שבו היא מתאפסת:

: <math>\frac{\partial \ell}{\partial \alpha} = \frac{n}{\alpha} + n \ln x_\mathrm{m} - \sum _{i=1}^n \ln x_i = 0.</math>

לפיכך משערך הנראות המרבית ל α הוא:

: <math>\widehat \alpha = \frac{n}{\sum _i \ln (x_i/\widehat x_\mathrm{m}) }.</math>

השגיאה הסטטיסטית הצפויה היא:<ref></ref>

: <math>\sigma = \frac {\widehat \alpha} {\sqrt n}. </math>

== דוגמאות ומקרים פרטיים ==
=== כללי ===
וילפרדו פארטו השתמש במקור בהתפלגות זו כדי לתאר את הקצאת העושר בין יחידים, כיוון שנראה היה שהיא מתארת היטב את האופן שבו בחברות שונות חלק גדול של העושר נמצא בבעלות אחוז קטן של האנשים. הוא גם השתמש בהתפלגות זו כדי לתאר את חלוקת ההכנסה. רעיון זה מתבטא לפעמים בצורה פשוטה יותר כעקרון פארטו או "כלל 80-20", שאומר ש-20% מהאוכלוסייה שולטת ב-80% מהעושר.<ref>For a two-quantile population, where approximately 18% of the population owns 82% of the wealth, the Theil index takes the value 1.</ref>

התפלגות זו אינה מוגבלת לתיאור עושר או הכנסה. הדוגמאות הבאות נתפסות לפעמים כבעלות התפלגות פארטו בקירוב:
* ארבעת המשתנים של מגבלת התקציב של משק הבית: צריכה, הכנסת עבודה, הכנסת הון ועושר.<ref></ref>
* הגדלים של יישובים אנושיים (מעט ערים, כפרים/כפרים רבים)<ref name="Reed"></ref><ref name="Reed2002"></ref>
* התפלגות גודלי קבצים בתעבורת אינטרנט המשתמשת בפרוטוקול TCP ( קבצים קטנים רבים, קבצים גדולים מעטים)<ref name="Reed" />

התאמת התפלגות פארטו מצטברת לכמות גשמים מקסימלית של יום אחד
* עתודות הנפט בשדות נפט
* תשואות של מניות בודדות<ref name="Reed" />
* גדלים של חלקיקי חול<ref name="Reed" />
* גדלים של מטאוריטים<ref></ref>
* בהידרולוגיה, התפלגות הפארטו מיושמת על אירועי קיצון כמו כמות גשמים מקסימלית מדי שנה של יום אחד וספיקה של נהרות.<ref>CumFreq, software for cumulative frequency analysis and probability distribution fitting </ref> התמונה הכחולה ממחישה דוגמה להתאמת התפלגות פארטו לכמות גשמים מקסימלית שנתית מדורגת של יום אחד, המציגה גם את רווח בר סמך של 90% בהתבסס על ההתפלגות הבינומית. נתוני הגשמים מיוצגים על ידי התווית מיקומים כחלק מניתוח התדירות המצטבר.

=== הקשר לחוק זיף ===
התפלגות פארטו היא התפלגות הסתברות רציפה. חוק זיף הוא התפלגות בדידה, המפרידה בין הערכים לדירוג פשוט. שניהם הם חוקי חזקות פשוטים עם מעריך שלילי, מנורמלים כך שההתפלגות המצטברת שלהם שווה ל-1. ניתן לגזור את חוק זיף מהתפלגות פארטו אם ערכי <math>x</math> ערכים (ההכנסות) מחולקים ל <math>N</math> דליים, כך שמספר האנשים בכל דלי מתנהג כהופכי של הדירוג שלו. ההתפלגות מנורמלת על ידי ההגדרה של <math>x_m</math> כך שיתקיים <math>\alpha x_\mathrm{m}^\alpha = \frac{1}{H(N,\alpha-1)}</math>, כאשר <math>H(N,\alpha-1)</math> הוא המספר ההרמוני המוכלל. כך, פונקציית צפיפות ההסתברות של זיף מהווה מקרה פרטי של התפלגות פארטו:

: <math>f(x) = \frac{\alpha x_\mathrm{m}^\alpha}{x^{\alpha+1}} = \frac{1}{x^s H(N,s)}</math>

כאשר <math>s = \alpha-1</math>, N הוא מדרגת ההכנסה הגבוהה ביותר ו-<math>x</math> הוא מספר שלם המייצג דירוג מ-1 עד N. אז לאדם שנבחר באקראי (או מילה, קישור לאתר או עיר) מאוכלוסייה (או שפה, אינטרנט או מדינה) יש הסתברות <math>f(x)</math> לדירוג <math>x</math> .

=== קשר ל"עקרון פארטו" ===
חוק "80–20", לפיו 20% מכלל האנשים מקבלים 80% מכלל ההכנסה, ו-20% מ-20% האמידים ביותר מקבלים 80% מאותם 80%, וכן הלאה, מתקיים כאשר מדד פארטו הוא <math>\alpha = \log_4 5 = \cfrac{\log_{10} 5}{\log_{10} 4} \approx 1.161</math> . תוצאה זו יכולה להיגזר מנוסחת עקומת לורנץ. יתר על כן, הוכח כי הטענות הבאות שקולות מתמטית<ref></ref>:
* ההכנסה מחולקת לפי התפלגות פארטו עם מדד α&nbsp;>&nbsp;1.
* קיים מספר <math>{1 \over 2} \geq p \geq 0</math> כך שחלק יחסי p מכלל האנשים מקבלים חלק יחסי <math>1-p</math> מההכנסה., ובאופן דומה עבור כל מספר ממשי 0<n (לא בהכרח מספר שלם) חלק יחסי <math>p^n</math> של כלל האנשים יקבל חלק יחסי <math>(1-p)^n</math> מכלל ההכנסה. α ו- p קשורים ב-

: <math>1-\frac{1}{\alpha}=\frac{\ln(1-p)}{\ln(p)}=\frac{\ln((1-p)^n)}{\ln(p^n)}</math> .

התפלגויות פארטו, שעבורן אשר, כפי שצוין לעיל, יש להן עתוחלת אינסופית, אינן יכולות לשמש כמודל סביר להתפלגות של הכנסה.

== ראו גם ==
* אפקט מתי
* עקרון פארטו
* יעילות פארטו
* חוק זיף
* התפלגות זנב עבה

==קישורים חיצוניים==

==הערות שוליים==

קטגוריה:התפלגויות רציפות