תיקון בסל

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

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

: <math>\frac n {n-1}</math>

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

: <math>(x_1-\overline{x},\,\dots,\,x_n-\overline{x}),</math>

כאשר <math>\overline{x}</math> הוא ממוצע המדגם. בעוד שיש n תצפיות בלתי תלויות במדגם, ישנם רק n-1 שארים בלתי תלויים במדגם, שכן השארים נסכמים ל-0.

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

== מינוח ==

: μ הוא הערך האמיתי של ממוצע האוכלוסייה (שאינו ידוע),

: <math>\overline{x}</math> הוא ממוצע המדגם,

: σ<sup>2</sup> הוא הערך האמיתי של שונות האוכלוסייה (שאינו ידוע),

: s<sub>n</sub><sup>2</sup> היא השונות המדגמית המוטה (ללא תיקון בסל),

: s<sup>2</sup> היא השונות המדגמית הלא מוטה (עם תיקון בסל),

=== תיקון בסל ===
ממוצע המדגם ניתן על ידי

:<math>\overline{x}=\frac{1}{n}\sum_{i=1}^n x_i.</math>

השונות המדגמית המוטה היא

:<math>s_n^2 = \frac {1}{n} \sum_{i=1}^n \left(x_i - \overline{x} \right)^ 2 = \frac{\sum_{i=1}^n \left(x_i^2\right)}{n} - \frac{\left(\sum_{i=1}^n x_i\right)^2}{n^2}</math>

והשונות המדגמית הלא מוטה (לאחר תיקון בסל) היא:

:<math>s^2 = \frac {1}{n-1} \sum_{i=1}^n \left(x_i - \overline{x} \right)^ 2 = \frac{\sum_{i=1}^n \left(x_i^2\right)}{n-1} - \frac{\left(\sum_{i=1}^n x_i\right)^2}{(n-1)n} = \left(\frac{n}{n-1}\right)\,s_n^2.</math>

== הוכחה ==
התוחלת של אי-ההתאמה בין האומד המוטה של השונות המדגמית לבין השונות האמיתית היא:

:<math>
\begin{align}
\operatorname{E} \left[ \sigma^2 - s_n^2 \right] &= \operatorname{E}\left[ \frac{1}{n} \sum_{i=1}^n(x_i - \mu)^2 - \frac{1}{n}\sum_{i=1}^n (x_i - \overline{x})^2 \right] \\
&= \operatorname{E}\left[ \frac{1}{n} \sum_{i=1}^n\left((x_i^2 - 2 x_i \mu + \mu^2) - (x_i^2 - 2 x_i \overline{x} + \overline{x}^2)\right) \right] \\
&= \operatorname{E}\left[ \frac{1}{n} \sum_{i=1}^n\left(\mu^2 - \overline{x}^2 + 2 x_i (\overline{x}-\mu) \right) \right] \\
&= \operatorname{E}\left[ \mu^2 - \overline{x}^2 + \frac{1}{n} \sum_{i=1}^n 2 x_i (\overline{x} - \mu) \right] \\
&= \operatorname{E}\left[ \mu^2 - \overline{x}^2 + 2(\overline{x} - \mu) \overline{x} \right] \\
&= \operatorname{E}\left[ \mu^2 - 2 \overline{x} \mu + \overline{x}^2 \right] \\
&= \operatorname{E}\left[ (\overline{x} - \mu)^2 \right] \\
&= \operatorname{Var} (\overline{x}) \\
&= \frac{\sigma^2}{n}
\end{align}
</math>

כאשר המעבר האחרון נובע מכך שהשונות של סכום המשתנים המקריים <math>x_1,x_2,...,x_n</math> היא סכום השונויות שלהם (שכן המשתנים הללו בלתי תלויים ולכן גם בלתי מתואמים), כך ש-<math>\operatorname{Var}(\overline{x})=\frac{n\cdot \sigma^2}{n^2} = \frac{\sigma^2}{n}</math>.

לפיכך, התוחלת של האומד המוטה של שונות המדגם היא:

:<math> \operatorname{E} \left[ s^2_n \right] = \sigma^2 - \frac{\sigma^2}{n} = \frac{n-1}{n} \sigma^2 </math>

ולכן, האומד הבלתי מוטה ניתן בנוסחה:

:<math> s^2 = \frac{n}{n-1} s_n^2 </math>.

=== אינטואיציה ===
עבור האומד המוטה, כאשר עושים שימוש בממוצע המדגם במקום בממוצע האמיתי (ממוצע האוכלוסייה), לוקים בהערכת חסר בשיעור <span style="text-decoration: overline">x</span>&nbsp;&minus;&nbsp;µ עבור כל איבר x<sub>i</sub>&nbsp;&minus;&nbsp;µ. מכיוון שהשונות של סכום היא סכום השונויות (בעבור משתנים בלתי מתואמים), אז, כדי למצוא את אי ההתאמה בין האומד המוטה לבין השונות האמיתית, יש למצוא את התוחלת של <math>(\bar{x}-\mu)^2</math>.

גודל זה הוא בדיוק השונות של ממוצע המדגם, שהיא σ<sup>2</sup>/n. לכן, אנו מצפים שהאומד המוטה לוקה בהערכת חסר של σ<sup>2</sup> בשיעור של σ<sup>2</sup>/n. כדי לפצות על הערכת חסר זאת, יש להפעיל את תיקון בסל, כלומר לכפול את השונות המוטה ב-<math>\frac{n}{n-1}</math>.

== מקורות ==
* Wolfram Mathworld - Bessel's Correction

== ראו גם ==
* שונות

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

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