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

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

הגדרה מתמטית

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

משפט הפירוק של פישר וניימן

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

$${\displaystyle f_{\theta }(x)=g_{\theta }(T(x))h(x)}$$

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

קל לראות שאם ${\displaystyle f(t)}$ היא פונקציה חד-חד-ערכית ו-${\displaystyle T}$ הוא סטטיסטי מספיק, אזי ${\displaystyle f(T)}$ הוא סטטיסטי מספיק גם כן. בפרט, ניתן לכפול סטטיסטי מספיק בקבוע שאינו אפס, ולקבל סטטיסטי מספיק.

משמעות לגבי הסקה סטטיסטית

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

ניסוח מדויק

יהיו ${\displaystyle X_{1},X_{2},\ldots ,X_{n}}$ תצפיות של מדגם מקרי מהתפלגות עם פונקציית צפיפות ${\displaystyle f(x,\theta )}$ עבור ${\displaystyle \theta \in \Theta }$, ויהי ${\displaystyle Y_{1}=u_{1}(X_{1},X_{2},\ldots ,X_{n})}$ סטטיסטי שפונקציית הצפיפות שלו היא ${\displaystyle g_{1}(y_{1};\theta )}$. אז ${\displaystyle Y_{1}=u_{1}(X_{1},X_{2},\ldots ,X_{n})}$ הוא סטטיסטי מספיק עבור ${\displaystyle \theta }$ אם ורק אם קיימת פונקציה ${\displaystyle H}$ כך ש:

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

הוכחה

ההוכחה להלן ניתנה על ידי רוברט הוג ואלן קרייג.^[1]

כיוון אחד

נניח ש:

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

ונוכיח ש-${\displaystyle Y_{1}}$ הוא סטטיסטי מספיק.

יהיו ${\displaystyle Y_{i}=u_{i}(X_{1},X_{2},\ldots ,X_{n})}$, ‏ ${\displaystyle i=2,3,\ldots ,n}$ סטטיסטים כלשהם שפונקציות הצפיפות שלהם הן ${\displaystyle g_{i}(y_{i};\theta )}$ בהתאמה.

נבצע את הטרנספורמציה: ${\displaystyle y_{i}=u_{i}(x_{1},x_{2},\ldots ,x_{n})}$ לכל ${\displaystyle i=1,2,\ldots ,n}$, כשהפונקציות ההפוכות הן: ${\displaystyle x_{j}=w_{j}(y_{1},y_{2},\ldots ,y_{n})}$ לכל ${\displaystyle j=1,2,\ldots ,n}$, ויעקוביאן ${\displaystyle J=[w_{j}/y_{i}]}$. אז:

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

הביטוי בצד שמאל של המשוואה הוא הצפיפות המשותפת: ${\displaystyle g(y_{1},y_{2},\dots ,y_{n};\theta )}$ של המשתנים (הסטטיסטים): ${\displaystyle Y_{1},Y_{2},\ldots ,Y_{n}}$.

בצד ימין של המשוואה, הביטוי ${\displaystyle g_{1}(y_{1};\theta )}$ הוא פונקציית הצפיפות של ${\displaystyle Y_{1}}$, וכך יוצא שהביטוי: ${\displaystyle H[w_{1},\dots ,w_{n}]\cdot |J|}$ שווה ל- ${\displaystyle g(y_{1},\dots ,y_{n};\theta )}$ חלקי ${\displaystyle g_{1}(y_{1};\theta )}$, כלומר, הוא שווה לפונקציית הצפיפות המותנה ${\displaystyle h(y_{2},\dots ,y_{n}\mid y_{1};\theta )}$ של ${\displaystyle Y_{2},\dots ,Y_{n}}$ בהינתן ${\displaystyle Y_{1}=y_{1}}$.

אבל לפי ההנחה, ${\displaystyle H(x_{1},x_{2},\dots ,x_{n})}$, וממילא גם ${\displaystyle H\left[w_{1}(y_{1},\dots ,y_{n}),\dots ,w_{n}(y_{1},\dots ,y_{n}))\right]}$ אינו תלוי ב-${\displaystyle \theta }$. מכיוון ש-${\displaystyle \theta }$ לא הוכנס בתוך הטרנספורמציה שביצענו, וכן לא בא לידי ביטוי ביעקוביאן ${\displaystyle J}$, יוצא ש-${\displaystyle h(y_{2},\dots ,y_{n}\mid y_{1};\theta )}$ אינו תלוי ב-${\displaystyle \theta }$ (אלא רק ב-${\displaystyle y_{1}}$, וזאת לכל סטטיסטים ${\displaystyle Y_{2},Y_{3},\ldots ,Y_{n}}$ אפשריים), ולכן ${\displaystyle Y_{1}}$ הוא סטטיסטי מספיק עבור ${\displaystyle \theta }$.

כיוון שני

נניח ש-${\displaystyle Y_{1}}$ הוא סטטיסטי מספיק, ונוכיח ש:

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

כיוון ש-ש-${\displaystyle Y_{1}}$ הוא סטטיסטי מספיק, ניתן לרשום:

${\displaystyle g(y_{1},\dots ,y_{n};\theta )=g_{1}(y_{1};\theta )h(y_{2},\dots ,y_{n}\mid y_{1})\,}$

כש-${\displaystyle h(y_{2},\dots ,y_{n}\mid y_{1})}$ אינו תלוי ב-${\displaystyle \theta }$ כיוון ש-${\displaystyle Y_{2}...Y_{n}}$ תלוי רק ב-${\displaystyle X_{1}...X_{n}}$, אשר הם בלתי תלויים ב-${\displaystyle \theta }$ בהינתן שידוע ${\displaystyle Y_{1}}$ שהוא סטטיסטי מספיק.

באמצעות חלוקת שני צידי המשוואה בערך המוחלט של היעקוביאן ${\displaystyle J}$, והחלפת ${\displaystyle y_{1},\dots ,y_{n}}$ בפונקציות: ${\displaystyle u_{1}(x_{1},\dots ,x_{n}),\dots ,u_{n}(x_{1},\dots ,x_{n})}$, מתקבלת המשוואה:

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

כש-${\displaystyle J*}$ הוא היעקוביאן עם ${\displaystyle y_{1},\dots ,y_{n}}$ שהוחלפו על ידי: ${\displaystyle x_{1},\dots ,x_{n}}$.

הצד הימני הוא בהכרח הצפיפות המשותפת ${\displaystyle \prod _{i=1}^{n}f(x_{i};\theta )}$ של ${\displaystyle X_{1},\dots ,X_{n}}$. נסמן: ${\displaystyle 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*|}}}$, ונקבל:

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

כיוון ש-${\displaystyle h(y_{2},\dots ,y_{n}\mid y_{1})}$ (ולכן גם ${\displaystyle h(u_{2},\dots ,u_{n}\mid u_{1})}$) אינו תלוי ב-${\displaystyle \theta }$ (לפי ההנחה), מתקבל ש:

${\displaystyle H(x_{1},\dots ,x_{2})={\frac {h(u_{2},\dots ,u_{n}\mid u_{1})}{|J*|}}}$

היא פונקציה שאינה תלויה ב-${\displaystyle \theta }$, כלומר:

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

כנדרש.

סטטיסטי מספיק מינימלי

סטטיסטי מספיק הוא מינימלי אם הוא יכול להיות מוצג כפונקציה של כל סטטיסטי מספיק אחר. במילים אחרות, ${\displaystyle S(X)}$ הוא סטטיסטי מספיק מינימלי אם ורק אם:

1. ${\displaystyle S(X)}$ הוא סטטיסטי מספיק
2. לכל סטטיסטי מספיק ${\displaystyle T(X)}$, קיימת פונקציה ${\displaystyle f}$ כך ש- ${\displaystyle S(X)=f(T(X))}$

באופן אינטואיטיבי, סטטיסטי מספיק מינימלי תופס באופן היעיל ביותר את כל המידע בנוגע לפרמטר ${\displaystyle \theta }$.

מאפיין שימושי של סטטיסטי מספיק מינימלי הוא שכאשר פונקציית הצפיפות ${\displaystyle f_{\theta }}$ קיימת, ${\displaystyle S(X)}$ הוא סטטיסטי מספיק מינימלי אם ורק אם:

${\displaystyle {\frac {f_{\theta }(x)}{f_{\theta }(y)}}}$ בלתי תלוי ב-${\displaystyle \theta }$ ${\displaystyle \Longleftrightarrow }$ ${\displaystyle S(x)=S(y)}$

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

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

בפרט, במרחב אוקלידי, התנאים הללו תמיד מתקיימים אם המשתנים המקריים (הקשורים ל- ${\displaystyle P_{\theta }}$) הם כולם בדידים או כולם רציפים.

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

אוסף היחסים של פונקציות נראות ${\displaystyle \left\{{\frac {L(\theta _{1}|X)}{L(\theta _{2}|X)}}\right\}}$ הוא סטטיסטי מספיק מינימלי אם ${\displaystyle P(X|\theta )}$ הוא בדיד או בעל פונקציית צפיפות.

דוגמאות

התפלגות ברנולי

אם ${\displaystyle X_{1},...,X_{n}}$ משתנים מקריים בלתי תלויים בעלי התפלגות ברנולי עם פרמטר ${\displaystyle p}$, אזי הסכום ${\displaystyle T(X)=X_{1}+...+X_{n}}$ הוא סטטיסטי מספיק עבור ${\displaystyle p}$.

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

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

ועל ידי קיבוץ חזקות של ${\displaystyle p}$ ושל ${\displaystyle 1-p}$, מתקבל:

${\displaystyle f_{p}(X)=p^{\sum x_{i}}(1-p)^{n-\sum x_{i}}=p^{T(x)}(1-p)^{n-T(x)}\,\!}$

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

${\displaystyle h(X)=1;\quad g(T(x),p)=p^{\sum x_{i}}(1-p)^{n-\sum x_{i}})}$

הפרמטר הבלתי ידוע, ${\displaystyle \ p}$, תלוי בנתונים (${\displaystyle X_{1},...,X_{n}}$) רק דרך הסטטיסטי שלהם ${\displaystyle T(x)={\sum _{i}x_{i}}}$

התפלגות אחידה

אם ${\displaystyle X_{1},...,X_{n}}$ משתנים מקריים בלתי תלויים המתפלגים באופן אחיד בטווח ${\displaystyle [0,\theta ]}$ אזי ${\displaystyle T(X)=\max(X_{1},...,X_{n})}$ הוא סטטיסטי מספיק ל-${\displaystyle \theta }$.

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

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

כאשר ${\displaystyle \mathbf {1} _{\{\cdots \}}}$ היא פונקציית האינדיקטור. לכן, אם נסמן: ${\displaystyle 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 \}}}$

וכן ${\displaystyle h(X)=\mathbf {1} _{\{0\leq \min(X_{i})\}}}$,

נקבל את תנאי משפט הפירוק, ושאכן ${\displaystyle T(X)=\max(X_{i})}$ הוא סטטיסטי מספיק.

התפלגות פואסון

אם ${\displaystyle X_{1},...,X_{n}}$ משתנים מקריים בלתי תלויים בעלי התפלגות פואסון עם פרמטר ${\displaystyle \lambda }$, אזי הסכום ${\displaystyle T(X)=X_{1}+...+X_{n}}$ הוא סטטיסטי מספיק עבור ${\displaystyle \lambda }$.

על מנת להיווכח בכך, נתבונן בפונקציית ההסתברות המשותפת:

${\displaystyle \Pr(X=x)=P(X_{1}=x_{1},X_{2}=x_{2}\ldots ,X_{n}=x_{n})\,}$

ומכיוון שהתצפיות בלתי תלויות:

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

לכן:

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

אם נסמן: ${\displaystyle g(T(x),\lambda )=e^{-n\lambda }\lambda ^{T(x)}\,}$, וכן: ${\displaystyle h(X)={1 \over x_{1}!x_{2}!\cdots x_{n}!}\,}$,

נקבל את תנאי משפט הפירוק, ושאכן ${\displaystyle T(X)=X_{1}+...+X_{n}}$ הוא סטטיסטי מספיק.

הערות שוליים

1. Hogg, Robert V.; Craig, Allen T. (1995). Introduction to Mathematical Statistics. Prentice Hall. ISBN 978-0-02-355722-4.