סטטיסטי מספיק

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

== הגדרה מתמטית ==
סטטיסטי <math>T(X)</math> הוא מספיק לפרמטר <math>\theta</math>, אם ההתפלגות המותנית של המדגם בהינתן הערך של <math>T(X)</math> אינה תלויה ב-<math>\theta</math>. במילים אחרות, הפונקציה <math>Pr_{\theta}(X=x | T(X)=t)</math> תלויה אולי ב-x ו-t, אבל לא ב-<math>\theta</math>.

== משפט הפירוק של פישר וניימן ==
משפט הפירוק של פישר וניימן (באנגלית: Fisher–Neyman factorization theorem) מספק אפיון נוח לכך שסטטיסטי הוא מספיק: אם פונקציית הצפיפות היא <math>f_\theta (x)</math>, אזי לפי המשפט <math>T</math> הוא מספיק עבור <math>\theta</math> אם ורק אם קיימות פונקציות אי-שליליות g ו-h כך ש:

<big><math display="block">f_\theta (x) = g_\theta(T(x))h(x)</math></big>

כלומר, ניתן לפרק את פונקציית הצפיפות למכפלה של שני גורמים, כך שגורם אחד, h, אינו תלוי ב-<math>\theta</math> והגורם השני, g, אשר כן תלוי ב-<math>\theta</math>, תלוי ב-x רק דרך <math>T(x)</math>.

קל לראות שאם <math>f(t)</math> היא פונקציה חד-חד-ערכית ו-<math>T</math> הוא סטטיסטי מספיק, אזי <math>f(T)</math> הוא סטטיסטי מספיק גם כן. בפרט, ניתן לכפול סטטיסטי מספיק בקבוע שאינו אפס, ולקבל סטטיסטי מספיק.

=== משמעות לגבי הסקה סטטיסטית ===
אחת ההשלכות של משפט הפירוק היא שכשמשתמשים בהסקה סטטיסטית על סמך נראות, שני מאגרי נתונים בגודל זהה אשר להם בדיוק אותו ערך עבור הסטטיסטי המספיק <math>T(x)</math>, תמיד יספקו את אותה מסקנה לגבי <math>\theta</math> (למשל: רווח סמך ל-<math>\theta</math>, דחיית/אי-דחיית השערה לגבי <math>\theta</math>). לפי קריטריון הפירוק התלות של הנראות ב-<math>\theta</math> היא רק ביחד עם <math>T(x)</math>. מכיוון שזה נכון לגבי שני מאגרי הנתונים שלהם אותו ערך <math>T(x)</math>, התלות של הנראות ב-<math>\theta</math> תהיה זהה גם כן, מה שיוביל להסקה סטטיסטית זהה.

=== ניסוח מדויק ===
יהיו <math>X_1, X_2, \ldots, X_n</math> תצפיות של מדגם מקרי מהתפלגות עם פונקציית צפיפות <math>f(x, \theta)</math> עבור <math>\theta \in \Theta</math>, ויהי <math>Y_1 = u_1(X_1, X_2, \ldots, X_n)</math> סטטיסטי שפונקציית הצפיפות שלו היא <math>g_1(y_1 ; \theta)</math>. אז <math>Y_1 = u_1(X_1, X_2, \ldots, X_n)</math> הוא סטטיסטי מספיק עבור <math>\theta</math> אם ורק אם קיימת פונקציה <math>H</math> כך ש:

<math display="block"> \prod_{i=1}^n f(x_i; \theta) = g_1 \left[u_1 (x_1, x_2, \dots, x_n); \theta \right] H(x_1, x_2, \dots, x_n)\,</math>

=== הוכחה ===
ההוכחה להלן ניתנה על ידי רוברט הוג ואלן קרייג.

==== כיוון אחד ====
נניח ש:
:<math> \prod_{i=1}^n f(x_i; \theta) = g_1 \left[u_1 (x_1, x_2, \dots, x_n); \theta \right] \cdot H(x_1, x_2, \dots, x_n)</math>

ונוכיח ש-<math>Y_1</math> הוא סטטיסטי מספיק.

יהיו <math>Y_i = u_i(X_1, X_2, \ldots, X_n)</math>, <math>i = 2, 3, \ldots, n</math> סטטיסטים כלשהם שפונקציות הצפיפות שלהם הן <math>g_i(y_i ; \theta)</math> בהתאמה.

נבצע את הטרנספורמציה: <math>y_i = u_i(x_1, x_2, \ldots, x_n)</math> לכל <math>i = 1, 2, \ldots, n</math>, כשהפונקציות ההפוכות הן: <math>x_j = w_j(y_1, y_2, \ldots, y_n)</math> לכל <math>j = 1, 2, \ldots, n</math>, ויעקוביאן <math>J = [w_j / y_i]</math>. אז:
:<math>
 \prod_{j=1}^n f \left[ w_j(y_1, y_2, \dots, y_n); \theta \right] =
 |J| \cdot g_1 (y_1; \theta) H \left[ w_1(y_1, y_2, \dots, y_n), \dots, w_n(y_1, y_2, \dots, y_n) \right]
</math>

הביטוי בצד שמאל של המשוואה הוא הצפיפות המשותפת: <math>g(y_1, y_2, \dots, y_n; \theta)</math> של המשתנים (הסטטיסטים): <math>Y_1, Y_2, \ldots, Y_n</math>.

בצד ימין של המשוואה, הביטוי <math>g_1(y_1;\theta)</math> הוא פונקציית הצפיפות של <math>Y_1</math>, וכך יוצא שהביטוי: <math>H[ w_1, \dots , w_n] \cdot |J|</math> שווה ל- <math>g(y_1,\dots,y_n;\theta)</math> חלקי <math>g_1(y_1;\theta)</math>, כלומר, הוא שווה לפונקציית הצפיפות המותנה <math>h(y_2, \dots, y_n \mid y_1; \theta)</math> של <math>Y_2,\dots,Y_n</math> בהינתן <math>Y_1 = y_1</math>.

אבל לפי ההנחה, <math>H(x_1,x_2,\dots,x_n)</math>, וממילא גם <math>H\left[w_1(y_1,\dots,y_n), \dots, w_n(y_1, \dots, y_n))\right]</math> אינו תלוי ב-<math>\theta</math>. מכיוון ש-<math>\theta</math> לא הוכנס בתוך הטרנספורמציה שביצענו, וכן לא בא לידי ביטוי ביעקוביאן <math>J</math>, יוצא ש-<math>h(y_2, \dots, y_n \mid y_1; \theta)</math> אינו תלוי ב-<math>\theta</math> (אלא רק ב-<math>y_1</math>, וזאת לכל סטטיסטים <math>Y_2, Y_3, \ldots, Y_n</math> אפשריים), ולכן <math>Y_1</math> הוא סטטיסטי מספיק עבור <math>\theta</math>.

==== כיוון שני ====
נניח ש-<math>Y_1</math> הוא סטטיסטי מספיק, ונוכיח ש:

:<math> \prod_{i=1}^n f(x_i; \theta) = g_1 \left[u_1 (x_1, x_2, \dots, x_n); \theta \right] \cdot H(x_1, x_2, \dots, x_n)\,</math>

כיוון ש-ש-<math>Y_1</math> הוא סטטיסטי מספיק, ניתן לרשום:

:<math>g(y_1,\dots,y_n;\theta)=g_1(y_1; \theta) h(y_2, \dots, y_n \mid y_1)\,</math>

כש-<math>h(y_2, \dots, y_n \mid y_1)</math> אינו תלוי ב-<math>\theta</math> כיוון ש-<math>Y_2 ... Y_n</math> תלוי רק ב-<math>X_1 ... X_n</math>, אשר הם בלתי תלויים ב-<math>\theta</math> בהינתן שידוע <math>Y_1</math> שהוא סטטיסטי מספיק.

באמצעות חלוקת שני צידי המשוואה בערך המוחלט של היעקוביאן <math>J</math>, והחלפת <math>y_1, \dots, y_n</math> בפונקציות: <math>u_1(x_1, \dots, x_n), \dots, u_n(x_1,\dots, x_n)</math>, מתקבלת המשוואה:

:<math>\frac{g\left[ u_1(x_1, \dots, x_n), \dots, u_n(x_1, \dots, x_n); \theta \right]}{|J*|}=g_1\left[u_1(x_1,\dots,x_n); \theta\right] \frac{h(u_2, \dots, u_n \mid u_1)}{|J*|}</math>

כש-<math>J*</math> הוא היעקוביאן עם <math>y_1,\dots,y_n</math> שהוחלפו על ידי: <math>x_1, \dots, x_n</math>.

הצד הימני הוא בהכרח הצפיפות המשותפת <math> \prod_{i=1}^n f(x_i; \theta)</math> של <math>X_1,\dots,X_n</math>. נסמן: <math>G_1\left[u_1(x_1,\dots,x_n); \theta\right] = \frac{g_1\left[u_1(x_1,\dots,x_n); \theta\right]}{|J*|}</math>, ונקבל:

:<math> \prod_{i=1}^n f(x_i; \theta) = G_1 \left[u_1 (x_1, x_2, \dots, x_n); \theta \right] \cdot h(u_2,\dots,u_n\mid u_1)</math>

כיוון ש-<math>h(y_2,\dots,y_n\mid y_1)</math> (ולכן גם <math>h(u_2,\dots,u_n\mid u_1)</math>) אינו תלוי ב-<math>\theta</math> (לפי ההנחה), מתקבל ש:

:<math>H(x_1,\dots,x_2)=\frac{h(u_2,\dots,u_n\mid u_1)}{|J*|}</math>

היא פונקציה שאינה תלויה ב-<math>\theta</math>, כלומר:

:<math> \prod_{i=1}^n f(x_i; \theta) = G_1 \left[u_1 (x_1, x_2, \dots, x_n); \theta \right] \cdot H(x_1, x_2, \dots, x_n)</math>

כנדרש.

== סטטיסטי מספיק מינימלי ==
סטטיסטי מספיק הוא מינימלי אם הוא יכול להיות מוצג כפונקציה של כל סטטיסטי מספיק אחר. במילים אחרות, <math>S(X)</math> הוא סטטיסטי מספיק מינימלי אם ורק אם:
# <math>S(X)</math> הוא סטטיסטי מספיק
# לכל סטטיסטי מספיק <math>T(X)</math>, קיימת פונקציה <math>f</math> כך ש- <math>S(X) = f(T(X))</math>
באופן אינטואיטיבי, סטטיסטי מספיק מינימלי תופס באופן היעיל ביותר את כל המידע בנוגע לפרמטר <math>\theta</math>.

מאפיין שימושי של סטטיסטי מספיק מינימלי הוא שכאשר פונקציית הצפיפות <math>f_\theta</math> קיימת, <math>S(X)</math> הוא סטטיסטי מספיק מינימלי אם ורק אם:

<math>\frac{f_\theta(x)}{f_\theta(y)}</math> בלתי תלוי ב-<math>\theta</math> <math>\Longleftrightarrow</math> <math>S(x) = S(y)</math>

זה נובע באופן ישיר ממשפט הפירוק של פישר-נימן.

מצב בו לא קיים סטטיסטי מספיק מינימלי הוצג על ידי הסטטיסטיקאי ההודי Raghu Raj Bahadur, בשנת 1954. אולם, תחת תנאים מתונים למדי, סטטיסטי מספיק מינימלי תמיד קיים.

בפרט, במרחב אוקלידי, התנאים הללו תמיד מתקיימים אם המשתנים המקריים (הקשורים ל- <math>P_\theta</math>) הם כולם בדידים או כולם רציפים.

אם קיים סטטיסטי מספיק מינימלי, אזי סטטיסטי מספיק שלם הוא בהכרח מינימלי. בעוד שקשה למצוא מקרים בהם סטטיסטי מספיק מינימלי אינו קיים, לא קשה למצוא מצבים בהם לא קיים סטטיסטי מספיק שלם.

אוסף היחסים של פונקציות נראות <math>\left\{\frac{L(\theta_1|X)}{L(\theta_2|X)}\right\}</math> הוא סטטיסטי מספיק מינימלי אם <math>P(X|\theta)</math> הוא בדיד או בעל פונקציית צפיפות.

== דוגמאות ==
=== התפלגות ברנולי ===
אם <math>X_1,...,X_n</math> משתנים מקריים בלתי תלויים בעלי התפלגות ברנולי עם פרמטר <math>p</math>, אזי הסכום <math>T(X) = X_1+...+X_n</math> הוא סטטיסטי מספיק עבור <math>p</math>.

מאחר שהמשתנים המקריים הם בלתי תלויים, פונקציית הצפיפות המשותפת מקיימת:

:<math>f_p(X) = p^{x_1}(1-p)^{1-x_1} p^{x_2}(1-p)^{1-x_2}\cdots p^{x_n}(1-p)^{1-x_n} \,\!</math>

ועל ידי קיבוץ חזקות של <math>p</math> ושל <math>1-p</math>, מתקבל:
:<math>f_p(X) = p^{\sum x_i}(1-p)^{n-\sum x_i} = p^{T(x)}(1-p)^{n-T(x)}\,\!</math>

הצגה זו מקיימת את תנאי משפט הפירוק כאשר:

<math>h(X) = 1 ; \quad g(T(x),p) = p^{\sum x_i}(1-p)^{n-\sum x_i}) </math>

הפרמטר הבלתי ידוע, <math>\ p</math>, תלוי בנתונים (<math>X_1,...,X_n</math>) רק דרך הסטטיסטי שלהם <math>T(x) = {\sum_i x_i}</math>

=== התפלגות אחידה ===
אם <math>X_1,...,X_n</math> משתנים מקריים בלתי תלויים המתפלגים באופן אחיד בטווח <math>[0,\theta]</math> אזי <math>T(X) = \max(X_1,...,X_n)</math> הוא סטטיסטי מספיק ל-<math>\theta</math>.

מאחר שהמשתנים המקריים הם בלתי תלויים, פונקציית הצפיפות המשותפת מקיימת:

:<math>\begin{align}
f_\theta(x_1,\ldots,x_n)
 &= \frac{1}{\theta}\mathbf{1}_{\{0\leq x_1\leq\theta\}} \cdots
 \frac{1}{\theta}\mathbf{1}_{\{0\leq x_n\leq\theta\}} \\
 &= \frac{1}{\theta^n}\mathbf{1}_{\{0\leq\min\{x_i\}\}}\mathbf{1}_{\{\max\{x_i\}\leq\theta\}}
\end{align}</math>

כאשר <math>\mathbf{1}_{\{\cdots\}}</math> היא פונקציית האינדיקטור. לכן, אם נסמן: <math>g(T(x),\theta) = \frac{1}{\theta^n} \cdot \mathbf{1}_{\{\max\{x_i\}\leq \theta\}} = \frac{1}{\theta^n} \cdot \mathbf{1}_{\{T(x)\leq \theta\}}</math>

וכן <math>h(X) = \mathbf{1}_{\{0\leq \min(X_i)\}}</math>,

נקבל את תנאי משפט הפירוק, ושאכן <math>T(X) = \max(X_i)</math> הוא סטטיסטי מספיק.

=== התפלגות פואסון ===
אם <math>X_1,...,X_n</math> משתנים מקריים בלתי תלויים בעלי התפלגות פואסון עם פרמטר <math>\lambda</math>, אזי הסכום <math>T(X) = X_1+...+X_n</math> הוא סטטיסטי מספיק עבור <math>\lambda</math>.

על מנת להיווכח בכך, נתבונן בפונקציית ההסתברות המשותפת:
:<math>
\Pr(X=x)=P(X_1=x_1,X_2=x_2\ldots,X_n=x_n)\,
</math>

ומכיוון שהתצפיות בלתי תלויות:
:<math>
P(X_1=x_1,X_2=x_2,\ldots,X_n=x_n) =
P(X_1=x_1) P(X_2=x_2) \ldots P(X_n=x_n)
</math>

לכן:
:<math>
P(X=x) =
{e^{-\lambda} \lambda^{x_1} \over x_1 !} \cdot
{e^{-\lambda} \lambda^{x_2} \over x_2 !} \cdots
{e^{-\lambda} \lambda^{x_n} \over x_n !} \, =
e^{-n\lambda} \lambda^{(x_1+x_2+\cdots+x_n)} \cdot
{1 \over x_1 ! x_2 !\cdots x_n ! } \, =
e^{-n\lambda} \lambda^{T(x)} \cdot
{1 \over x_1 ! x_2 !\cdots x_n ! } \,
</math>

אם נסמן: <math>g(T(x),\lambda) = e^{-n\lambda} \lambda^{T(x)} \,</math>, וכן: <math>h(X) = {1 \over x_1 ! x_2 !\cdots x_n ! } \,</math>,

נקבל את תנאי משפט הפירוק, ושאכן <math>T(X) = X_1+...+X_n</math> הוא סטטיסטי מספיק.

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

קטגוריה:סטטיסטיקה