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

התפלגות פארטו, הנקראת על שם המהנדס האזרחי, הכלכלן והסוציולוג האיטלקי וילפרדו פארטו, היא התפלגות הסתברות המשמשת לתיאור של תופעות מתחומים מגוונים: חברתיות, בקרת איכות, מדעיות, גיאופיזיות, אקטואריות ועוד. העיקרון, שיושם במקור לתיאור של חלוקת העושר בחברה, מתאים לתופעה לפיה חלק גדול מהעושר מוחזק על ידי חלק קטן מהאוכלוסייה. עקרון פארטו או "כלל 80-20", הקובע ש-80% מהתוצאות נובעות מ-20% מהסיבות, נקרא על שם פארטו, אך רק התפלגויות פארטו עם פרמטר צורה α שערכו $\log _{4}5\approx 1.16$ ישקפו במדויק את עיקרון פארטו. תצפיות אמפיריות הראו שהתפלגות 80-20 זו מתאימה למגוון רחב של תופעות טבע והתנהגות אנושיות.

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

פונקציית צפיפות ההסתברות
px
פונקציית ההסתברות המצטברת
px
מאפיינים
פרמטרים	$x_{\mathrm {m} }>0$ (ממשי) $\alpha >0$ (ממשי)
תומך	$x\in [x_{\mathrm {m} },\infty )$
פונקציית צפיפות הסתברות (pdf)	${\frac {\alpha x_{\mathrm {m} }^{\alpha }}{x^{\alpha +1}}}$
פונקציית ההסתברות המצטברת (cdf)	$1-\left({\frac {x_{\mathrm {m} }}{x}}\right)^{\alpha }$
תוחלת	${\begin{cases}\infty &{\text{for }}\alpha \leq 1\\{\dfrac {\alpha x_{\mathrm {m} }}{\alpha -1}}&{\text{for }}\alpha >1\end{cases}}$
חציון	$x_{\mathrm {m} }{\sqrt[{\alpha }]{2}}$
ערך שכיח	$x_{\mathrm {m} }$
שונות	${\begin{cases}\infty &{\text{for }}\alpha \leq 2\\{\dfrac {x_{\mathrm {m} }^{2}\alpha }{(\alpha -1)^{2}(\alpha -2)}}&{\text{for }}\alpha >2\end{cases}}$
אנטרופיה	$\log \left(\left({\frac {x_{\mathrm {m} }}{\alpha }}\right)\,e^{1+{\tfrac {1}{\alpha }}}\right)$
צידוד	${\frac {2(1+\alpha )}{\alpha -3}}{\sqrt {\frac {\alpha -2}{\alpha }}}{\text{ for }}\alpha >3$
גבנוניות	${\frac {6(\alpha ^{3}+\alpha ^{2}-6\alpha -2)}{\alpha (\alpha -3)(\alpha -4)}}{\text{ for }}\alpha >4$

הגדרות

אם X הוא משתנה מקרי עם התפלגות פארטו (סוג I), אזי ההסתברות ש- X גדול ממספר כלשהו x, כלומר, פונקציית ההישרדות (נקראת גם פונקציית זנב), נתונה על ידי

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

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

פונקציית התפלגות מצטברת

מההגדרה, פונקציית ההתפלגות המצטברת של משתנה מקרי שמתפלג פארטו עם הפרמטרים α ו- x_m היא

$F_{X}(x)={\begin{cases}1-\left({\frac {x_{\mathrm {m} }}{x}}\right)^{\alpha }&x\geq x_{\mathrm {m} },\\0&x<x_{\mathrm {m} }.\end{cases}}$

פונקציית צפיפות הסתברות

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

$f_{X}(x)={\begin{cases}{\frac {\alpha x_{\mathrm {m} }^{\alpha }}{x^{\alpha +1}}}&x\geq x_{\mathrm {m} },\\0&x<x_{\mathrm {m} }.\end{cases}}$

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

מאפיינים

מומנטים ופונקציה אופיינית

• התוחלת של משתנה מקרי שמתפלג פארטו היא

$\operatorname {E} (X)={\begin{cases}\infty &\alpha \leq 1,\\{\frac {\alpha x_{\mathrm {m} }}{\alpha -1}}&\alpha >1.\end{cases}}$

• השונות של משתנה מקרי שמתפלג פארטו היא

$\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}}$

(עבור α ≤ 1, השונות אינה מוגדרת)

• המומנטים הגולמיים הם

$\mu _{n}'={\begin{cases}\infty &\alpha \leq n,\\{\frac {\alpha x_{\mathrm {m} }^{n}}{\alpha -n}}&\alpha >n.\end{cases}}$

• הפונקציה היוצרת מוגדרת רק עבור ערכים שליליים או אפס $t\leq 0$

$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)$

$M\left(0,\alpha ,x_{\mathrm {m} }\right)=1.$

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

• הפונקציה האופיינית נתונה על ידי

$\varphi (t;\alpha ,x_{\mathrm {m} })=\alpha (-ix_{\mathrm {m} }t)^{\alpha }\Gamma (-\alpha ,-ix_{\mathrm {m} }t),$

כאשר$\Gamma (a,x)$ היא פונקציית גמא הלא שלמה.

ניתן לחלץ את הפרמטרים באמצעות שיטת המומנטים.

התפלגויות מותנות

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

${\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}}$

הביטוי מרמז כל כך שהערך התוחלת המותנית (אם היא סופית, כלומר $\alpha >1$ ) פרופורציוני ל-$x_{1}$ :

${\text{E}}(X|X\geq x_{1})\propto x_{1}$

במקרה של משתנים אקראיים המתארים את משך חייו של אובייקט, המשמעות היא שתוחלת החיים היא פרופורציונלית לגיל, והיא נקראת אפקט לינדי או חוק לינדי.

משפט אפיון

נניח כי $X_{1},X_{2},X_{3},\dotsc$ הם משתנים מקריים בלתי תלויים אשר התפלגות ההסתברות שלהם נתמכת בקטע $[x_{\text{m}},\infty )$ עבור $x_{\text{m}}>0$ כלשהו. נניח שלכל $n$, שני המשתנים האקראיים $\min\{X_{1},\dotsc ,X_{n}\}$ ו $(X_{1}+\dotsb +X_{n})/\min\{X_{1},\dotsc ,X_{n}\}$ הם בלתי תלויים. אז ההתפלגות המשותפת היא התפלגות פארטו.

ממוצע גאומטרי

הממוצע הגאומטרי (G) הוא

$G=x_{\text{m}}\exp \left({\frac {1}{\alpha }}\right).$

ממוצע הרמוני

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

$H=x_{\text{m}}\left(1+{\frac {1}{\alpha }}\right).$

ייצוג גרפי

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

$\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.$

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

הסקה סטטיסטית

הערכת פרמטרים

פונקציית הנראות של פרמטרי התפלגות פארטו α ו- x_m, בהינתן מדגם בלתי תלוי x = (x₁ , x₂ , ..., x_n), הוא

$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}}}.$

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

$\ell (\alpha ,x_{\mathrm {m} })=n\ln \alpha +n\alpha \ln x_{\mathrm {m} }-(\alpha +1)\sum _{i=1}^{n}\ln x_{i}.$

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

${\widehat {x}}_{\mathrm {m} }=\min _{i}{x_{i}}.$

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

${\frac {\partial \ell }{\partial \alpha }}={\frac {n}{\alpha }}+n\ln x_{\mathrm {m} }-\sum _{i=1}^{n}\ln x_{i}=0.$

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

${\widehat {\alpha }}={\frac {n}{\sum _{i}\ln(x_{i}/{\widehat {x}}_{\mathrm {m} })}}.$

השגיאה הסטטיסטית הצפויה היא:

$\sigma ={\frac {\widehat {\alpha }}{\sqrt {n}}}.$

דוגמאות ומקרים פרטיים

כללי

וילפרדו פארטו השתמש במקור בהתפלגות זו כדי לתאר את הקצאת העושר בין יחידים, כיוון שנראה היה שהיא מתארת היטב את האופן שבו בחברות שונות חלק גדול של העושר נמצא בבעלות אחוז קטן של האנשים. הוא גם השתמש בהתפלגות זו כדי לתאר את חלוקת ההכנסה. רעיון זה מתבטא לפעמים בצורה פשוטה יותר כעקרון פארטו או "כלל 80-20", שאומר ש-20% מהאוכלוסייה שולטת ב-80% מהעושר.

התפלגות זו אינה מוגבלת לתיאור עושר או הכנסה. הדוגמאות הבאות נתפסות לפעמים כבעלות התפלגות פארטו בקירוב:

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

התאמת התפלגות פארטו מצטברת לכמות גשמים מקסימלית של יום אחד

• עתודות הנפט בשדות נפט
• תשואות של מניות בודדות
• גדלים של חלקיקי חול
• גדלים של מטאוריטים
• בהידרולוגיה, התפלגות הפארטו מיושמת על אירועי קיצון כמו כמות גשמים מקסימלית מדי שנה של יום אחד וספיקה של נהרות. התמונה הכחולה ממחישה דוגמה להתאמת התפלגות פארטו לכמות גשמים מקסימלית שנתית מדורגת של יום אחד, המציגה גם את רווח בר סמך של 90% בהתבסס על ההתפלגות הבינומית. נתוני הגשמים מיוצגים על ידי התווית מיקומים כחלק מניתוח התדירות המצטבר.

הקשר לחוק זיף

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

$f(x)={\frac {\alpha x_{\mathrm {m} }^{\alpha }}{x^{\alpha +1}}}={\frac {1}{x^{s}H(N,s)}}$

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

קשר ל"עקרון פארטו"

חוק "80–20", לפיו 20% מכלל האנשים מקבלים 80% מכלל ההכנסה, ו-20% מ-20% האמידים ביותר מקבלים 80% מאותם 80%, וכן הלאה, מתקיים כאשר מדד פארטו הוא $\alpha =\log _{4}5={\cfrac {\log _{10}5}{\log _{10}4}}\approx 1.161$ . תוצאה זו יכולה להיגזר מנוסחת עקומת לורנץ. יתר על כן, הוכח כי הטענות הבאות שקולות מתמטית:

• ההכנסה מחולקת לפי התפלגות פארטו עם מדד α > 1.
• קיים מספר ${1 \over 2}\geq p\geq 0$ כך שחלק יחסי p מכלל האנשים מקבלים חלק יחסי $1-p$ מההכנסה., ובאופן דומה עבור כל מספר ממשי 0<n (לא בהכרח מספר שלם) חלק יחסי $p^{n}$ של כלל האנשים יקבל חלק יחסי $(1-p)^{n}$ מכלל ההכנסה. α ו- p קשורים ב-

$1-{\frac {1}{\alpha }}={\frac {\ln(1-p)}{\ln(p)}}={\frac {\ln((1-p)^{n})}{\ln(p^{n})}}$ .

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

ראו גם

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

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

מדיה וקבצים בנושא התפלגות פארטו בוויקישיתוף

• התפלגות פארטו, באתר MathWorld (באנגלית)