ViewWiki

העולם, קצת יותר מקרוב

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

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

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

01 / הניסוי

מצאו את הישר שבתוך הענן.

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

תמונה של הסצנה התלת־ממדית „סטטיסטיקה ורגרסיה”, שעוסקת בשאלה: איך מוצאים את הישר שמסתתר בענן של נקודות, ואיך אומרים כמה טוב הוא מתאים?
סצנה אינטראקטיבית

הגרסה האינטראקטיבית נטענת…

ההסבר שלמטה זמין כבר עכשיו, בזמן שהסצנה נטענת.

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

02 / שלושה ענני נקודות

תלול זה לא אותו דבר כמו חזק.

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

תלול ורועש01

תלול אבל מפוזר

y עולה ביחידה שלמה על כל יחידה של x, אבל דחיפות אקראיות גדולות מפזרות את הנקודות. הישר תלול; ההתאמה רק סבירה.

כלל
שיפוע 1, רעש 2
הישר הטוב ביותר
ŷ = −0.26 + 1.05x
r
0.79
מה r אומר
התאמה סבירה: הישר מוריד בערך שלוש חמישיות מהריבועים
מתון וצמוד02

מתון אבל צמוד

y עולה רק ברבע מהקצב, אבל כל נקודה יושבת כמעט על הישר. גם שיפוע מתון יכול לבוא עם r קרוב ל־1.

כלל
שיפוע 0.25, רעש 0.15
הישר הטוב ביותר
ŷ = 3.73 + 0.25x
r
0.97
מה r אומר
כמעט התאמה מושלמת, אף שהישר כמעט שטוח
עקום03

קשת

y עולה ויורד שוב, כמעט בלי רעש. הישר הטוב ביותר כמעט אופקי ו־r קרוב ל־0, אף ש־x קובע את y כמעט לגמרי.

כלל
y עולה ואז יורד
הישר הטוב ביותר
ŷ = 7.51 − 0.01x
r
−0.02
מה r אומר
אין קשר של ישר, אף ש־y תלוי ב־x כמעט לגמרי
הרעיון הגדול

מקדם המתאם מודד כמה הנקודות צמודות לישר: לא כמה הישר תלול, ולא אם יש בין x ל־y קשר מסוג אחר. הריבוע שלו הוא החלק מהפיזור שהישר מוריד.

r² = 1 − sₑ² ⁄ s²

ענן של מדידות

מדדו שתי תכונות של אותם עצמים, ושרטטו אחת מול השנייה: הגובה ומוטת הידיים של עשרים אנשים, הטמפרטורה ומכירות הגלידה של עשרים ימים, המשקל שתלוי על קפיץ וכמה הוא נמתח. כל עצם הופך לנקודה אחת, (x,y)(x, y)⁠, ויחד הן יוצרות דיאגרמת פיזור. אם הנקודות עולות משמאל לימין, xx⁠ גדול יותר נוטה לבוא עם yy⁠ גדול יותר; אם הן יורדות, עם yy⁠ קטן יותר. אבל מדידות אמיתיות אף פעם לא יושבות בדיוק על ישר. משהו תמיד דוחף כל אחת מהן קצת: סיבות אחרות, מקרה, הגבולות של סרגל המדידה.

הלוח יוצר ענן כזה מכלל שאתם קובעים. כל נקודה מתחילה על הישר y=5+m(x−5)y = 5 + m(x - 5)⁠, כש־mm⁠ הוא השיפוע, ואז נדחפת למעלה או למטה בשיעור אקראי. שיפוע קובע כמה תלול הכלל עולה או יורד; רעש קובע את גודל הדחיפות האקראיות, את סטיית התקן של ההגרלות. כשהרעש הוא 00⁠ כל הנקודות יושבות על הכלל. הסטטיסטיקה שואלת בכיוון ההפוך: כשיש בידינו רק את הנקודות, איך מוצאים את הכלל ואומרים כמה אפשר לסמוך עליו?

האמצע והפיזור

התחילו מערכי yy⁠ לבדם. הממוצע שלהם, yˉ\bar{y}⁠, הוא הסכום שלהם חלקי מספרם. קדמו את התאמת הישר אל ממוצע ופיזור, וקו אופקי נדלק בגובה yˉ\bar{y}⁠. הממוצע הוא נקודת האיזון של הנתונים: סכום המרחקים של הנקודות שמעליו שווה בדיוק לסכום המרחקים של הנקודות שמתחתיו, כך שאילו כל נקודה הייתה משקולת על סרגל, הסרגל היה מתאזן שם. הטבעת מסמנת את נקודת האיזון של הענן כולו, (xˉ,yˉ)(\bar{x}, \bar{y})⁠.

כמה רחוק הנקודות סוטות מהממוצע? המרחק של כל אחת מהקו, הסטייה שלה, תלוי עכשיו ממנו כמקל גיר, ועליו מצויר ריבוע. סטיות למעלה ולמטה מבטלות זו את זו כשמחברים אותן, ולכן הסטטיסטיקה מעלה אותן קודם בריבוע. הממוצע של הריבועים הוא השונות, s2s^2⁠, והשורש שלה הוא סטיית התקן, ss⁠ (בספרי לימוד רבים σ\sigma⁠):

s=∑(y−yˉ)2ns = \sqrt{\frac{\sum (y - \bar{y})^2}{n}}⁠

כלומר ss⁠ היא הצלע של הריבוע הממוצע, שעומד ליד הגרף. הרצועה המוצללת מ־yˉ−s\bar{y} - s⁠ עד yˉ+s\bar{y} + s⁠ מכילה את רוב הנקודות. בערכים שמתפלגים לפי עקומת הפעמון היא מכילה בערך שתיים מכל שלוש, ורצועה רחבה פי שניים מכילה בערך 95 מכל מאה.

הישר שמחטיא הכי מעט

עכשיו קדמו אל ישר הרגרסיה. הקו האופקי מסתובב סביב נקודת האיזון, וכל מקל עוקב אחריו, כך שכל מקל הופך לשארית: בכמה הישר מחטיא את הנקודה, כשמודדים ישר למעלה או למטה. ככל שהישר נוטה לכיוון המגמה הריבועים מתכווצים, והוא נעצר במקום שבו השטח הכולל שלהם הכי קטן שאפשר. זה ישר הריבועים הפחותים, או ישר הרגרסיה, שכותבים y^=a+bx\hat{y} = a + bx⁠: ישר, הפונקציה הקווית של האלגברה בבית הספר, עם שיפוע bb⁠ וחיתוך aa⁠.

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

b=∑(x−xˉ)(y−yˉ)∑(x−xˉ)2,a=yˉ−bxˉb = \frac{\sum (x - \bar{x})(y - \bar{y})}{\sum (x - \bar{x})^2}, \qquad a = \bar{y} - b\bar{x}⁠

ולכן ישר הרגרסיה עובר תמיד דרך (xˉ,yˉ)(\bar{x}, \bar{y})⁠. בערכי ההתחלה הלוח כותב y^=1.90+0.62x\hat{y} = 1.90 + 0.62x⁠, קרוב לכלל הנסתר, שהשיפוע שלו 0.60.6⁠. זה לא בדיוק הכלל, כי עשרים ההגרלות האלה של הרעש מטות את הענן במקרה מעט. מדדו עוד עשרים, והישר ייצא קצת אחר; כמה הוא מתנדנד, ואיך הנדנוד קטן ככל שהמדגם גדל, זה הנושא של מודלים ואי־ודאות.

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

מקדם המתאם r

הריבועים מודדים גם כמה טוב הישר מתאים. לפני הישר, השטח הממוצע של הריבועים סביב הממוצע הוא s2s^2⁠. אחריו, לריבועים של השאריות יש שטח ממוצע קטן יותר, se2s_e^2⁠. החלק מהשטח שהישר מוריד נקרא r2r^2⁠:

r2=1−se2s2r^2 = 1 - \frac{s_e^2}{s^2}⁠

והשורש שלו, עם הסימן של השיפוע, הוא מקדם המתאם rr⁠. הוא נע בין −1-1⁠ ל־+1+1⁠. ב־+1+1⁠ או ב־−1-1⁠ כל הנקודות על ישר אחד ולא נשארים ריבועים; קרוב ל־00⁠ הישר הטוב ביותר כמעט אופקי ומוריד כמעט כלום. בערכי ההתחלה r=0.89r = 0.89⁠, ולכן הישר מוריד 78% מהשטח. אותו מספר מתקבל ישירות מהנתונים:

r=∑(x−xˉ)(y−yˉ)∑(x−xˉ)2⋅∑(y−yˉ)2r = \frac{\sum (x - \bar{x})(y - \bar{y})}{\sqrt{\sum (x - \bar{x})^2 \cdot \sum (y - \bar{y})^2}}⁠

ככלל אצבע, מתאם שגודלו מעל 0.70.7⁠ נקרא חזק, ומתחת ל־0.30.3⁠ חלש. קרל פירסון נתן לנוסחה את הצורה הזאת ב־1896, על בסיס העבודה של פרנסיס גלטון על תורשה; גלטון גם נתן לרגרסיה את שמה, אחרי שהבחין שלהורים גבוהים נוטים להיוולד ילדים גבוהים, אבל פחות מהם, קרובים יותר לממוצע.

נסו: הזיזו את רעש אל 00⁠, וכל הנקודות נופלות על הישר: r=1r = 1⁠. עכשיו הזיזו אותו אל 22⁠. השיפוע של הישר הטוב ביותר כמעט לא משתנה, אבל הריבועים גדלים ו־rr⁠ יורד ל־0.630.63⁠. אחר כך קבעו את שיפוע על 0.20.2⁠ ואת רעש על 00⁠: ישר מתון, ובכל זאת rr⁠ שוב 11⁠. המקדם rr⁠ מודד כמה הנקודות צמודות לישר, ולא כמה הישר תלול.

חזרה לניסוי

ערכים חריגים, עקומות וסיבות

הפעילו את ערך חריג, ומופיעה מדידה שגויה אחת, רחוקה מתחת לכל השאר. השארית שלה גדולה, ולכן הריבוע שלה ענק, ומכיוון ששיטת הריבועים הפחותים שוקלת כל ריבוע באותה מידה, הנקודה הבודדת הזאת מושכת את הישר אליה: השיפוע יורד מ־0.620.62⁠ ל־0.470.47⁠, ו־rr⁠ מ־0.890.89⁠ ל־0.690.69⁠. ניתוח אמיתי היה בוחן נקודה כזאת היטב לפני שהוא סומך על הישר, כי טעות בהעתקה או חיישן מקולקל יכולים לגבור על תשע־עשרה מדידות טובות.

המקדם rr⁠ מודד רק קשר של ישר. נקודות על קשת מושלמת, כמו הגובה של כדור שנזרק כשמודדים אותו במרווחים שווים מהזריקה ועד התפיסה, נותנות rr⁠ קרוב ל־00⁠, אף ש־yy⁠ תלוי ב־xx⁠ לגמרי. הסתכלו תמיד בדיאגרמת הפיזור, ולא רק במספר.

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

מה המודל הזה משמיט

הרעש הוא קבוצה קבועה אחת של עשרים הגרלות אקראיות מעקומת הפעמון, שהותאמו כך שהממוצע שלהן יהיה בדיוק 00⁠ וסטיית התקן שלהן בדיוק 11⁠, ולכן רעש הוא בדיוק הפיזור של הנקודות סביב הכלל. שאר התכונות שלהן נשארו כפי שהמקרה קבע, ולכן הישר המותאם קרוב לכלל אבל לא עליו. רק ערכי yy⁠ רועשים: כל xx⁠ נחשב כאילו נמדד בדיוק מושלם. בנתונים אמיתיים יש לעתים קרובות שגיאות בשניהם. והשיטה לא סימטרית בכל מקרה: חיזוי של xx⁠ מתוך yy⁠ נותן ישר אחר מחיזוי של yy⁠ מתוך xx⁠, אלא אם rr⁠ שווה ±1.

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

מותר לעשות בהם שימוש חוזר, עם קרדיט ל־ViewWiki וקישור לעמוד הזה.

04 / ללמוד את זה בקורס

תרגול, צעד אחר צעד.

5 השיעורים האלה מלמדים את הנושא עם שאלות שנבדקות, לפי סדר הקורס.

  1. ממוצע, חציון ושכיח מתמטיקה · סטטיסטיקה
  2. טווח וסטיית תקן מתמטיקה · סטטיסטיקה
  3. מתאם מתמטיקה · סטטיסטיקה
  4. ישר הרגרסיה מתמטיקה · סטטיסטיקה
  5. ללמוד ישר מדעי המחשב · למידת מכונה