מבחן גריינג'ר

ממוזער|348x348px|כאשר סדרה עתית X גריינג'ר-מסבירה את סדרה עתית Y, כיוון ועוצמת השינוי בערכי X, יופיעו בקירוב גם בסדרה Y, לעיתים בפיגור.
**מבחן גריינג'ר** הוא מבחן סטטיסטי הבודק את כוחה של סדרה עתית אחת בניבוי סדרה עתית אחרת. מבחן זה עושה שימוש ברגרסיה ליניארית עם משתנים בפיגור ובמבחני מובהקות שונים כגון מבחן F, קריטריון האינפורמציה של Akaike וקריטריון בייס. המבחן בודק את ההשערה כי תוספת של ערכים בפיגור של משתנה מסביר X למשוואה המכילה ערכים בפיגור של משתנה מוסבר Y, תגדיל את כוח ההסבר של המודל. במקרה שהשערה זו מתקבלת, נהוג לומר כי "X גריינג'ר-מסביר את Y".

המבחן פותח על ידי החוקר קלייב גריינג'ר בשנת 1969, והפך לרווח בעיקר בתחומי הכלכלה. עם השנים הפך המבחן לשימושי ומקובל גם בתחומי המדעים המדויקים כגון מדעי המוח, הנדסה ומדעי המחשב.

## הגדרה

### קשר סיבתי חד כיווני פשוט
תהיה Y_t סדרה עתית סטאציונריות בה מתקיים קשר אוטו-רגרסיבי פשוט: y_t=a+\sum_{j=1}^{K_y}b_jy_{t-j}, עבור K_y כלשהו, ותהיה X_t סדרה עתית סטציונרית נוספת.

עבור U ,מרחב האינפורמציה שנצברה עבור כלל המשתנים במודל עד לתקופה t, ומרחבי האינפורמציה של הסדרות X_t ו- Y_t, נסמן את שני תתי-המרחבים הבאים:
* מרחב האינפורמציה שנצברה אודות המשתנה בעבר: \overline{A_i}:=\{A_{t-j}: j=1,2,3,...\}
* מרחב האינפורמציה שנצברה אודות המשתנה בעבר ובנקודה הנוכחית: \overline{\overline{A_i}}:=\{A_{t-j}: j=0,1,2,...\}
נאמר כי X גריינג'ר-מסביר את Y אם ורק אם: \sigma^2(Y|{U}).

כלומר, שימוש בכלל האינפורמציה במרחב U עדיף על שימוש בכלל האינפורמציה אותו המרחב, *למעט* האינפורמציה אודות X.

מבחינת משוואות האמידה, נאמר כי המשוואה: y_t=a+\sum_{j=1}^{K_y}b_jy_{t-j}+\sum_{j=1}^{K_x}b_jx_{t-j} בעלת כוח הסבר רב יותר מן המשוואה: y_t=a+\sum_{j=1}^{K_y}b_jy_{t-j}.

## הרחבת המודל

### קשרים חד כיווניים

#### קשר סיבתי מיידי
ההבדל העיקרי בין קשר סיבתי פשוט לקשר סיבתי מיידי, הוא בבחינת כוח ההסבר של ערך הסדרה X_t בנקודת הזמן בה אומדים את Y_t. על כן, נאמר כי X גריינג'ר-מסביר באופן מיידי את Y אם ורק אם: \sigma^2(Y|{\overline{U},\overline{\overline{X}}}).

כלומר, שימוש בכלל האינפורמציה במרחב \overline{U}, עדיף על שימוש בכלל האינפורמציה אותו המרחב, *למעט* האינפורמציה אודות \overline{\overline{X}}.

מבחינת משוואות האמידה, נאמר כי המשוואה: y_t=a+\sum_{j=1}^{K_y}b_jy_{t-j}+\sum_{j=0}^{K_x}b_jx_{t-j} בעלת כוח הסבר רב יותר מן המשוואה: y_t=a+\sum_{j=1}^{K_y}b_jy_{t-j}.

#### קשר סיבתי בפיגור
עבור הסדרה X_t נגדיר את הסדרה X(m):=\{X_{t-m}, X_{t-(m+1)}, X_{t-(m-2)},...\}, כאשר: m\in \N\backslash{\{0,1\}}.

נאמר כי X גריינג'ר-מסביר בפיגור m את Y, אם ורק אם m הוא המספר הטבעי הקטן ביותר עבורו: \sigma^2(Y|{U\backslash X(m)}).

כלומר, שימוש בערכים: X_t,X_{t-1},...,X_{t-(m-1)} לא תורם לכח ההסבר של X את Y.

מבחינת משוואת האמידה, נאמר כי המשוואה: y_t=a+\sum_{j=1}^{K_y}b_jy_{t-j}+\sum_{j=m}^{K_x}b_jx_{t-j} בעלת כוח הסבר רב יותר מהמשוואה: y_t=a+\sum_{j=1}^{K_y}b_jy_{t-j}+\sum_{j=0}^{K_x}b_jx_{t-j}, וגם מן המשוואה: y_t=a+\sum_{j=1}^{K_y}b_jy_{t-j}.

### קשרים דו-כיווניים
מבחן גריינג'ר מתיר גם הרחבה של המודל, לבדיקה של קשרים דו-כיוונים של הסדרות  X_t ו- Y_t זו על זו. לשם כך נניח כעת כי עבור שתי הסדרות העתיות, X_t,Y_t, קיים מודל אוטו-רגרסיבי מובהק. דהיינו, קיים צמד מודלים, בעלי כח הסבר מתקבל על הדעת, מן הצורה:

y_t=a+\sum_{j=1}^{K_y}b_jy_{t-j} 

x_t=a+\sum_{j=1}^{K_x}b_jx_{t-j}

### היזון חוזר
תחת אותן ההגדרות ואותם הסימונים, נאמר כי קיים היזון חוזר בין X ל- Y אם ורק אם מתקיימים יחדיו שני התנאים הבאים:

\sigma^2(Y|{U})

\sigma^2(X|{U})

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

y_t=a+\sum_{j=1}^{K_y}b_jy_{t-j}+\sum_{j=1}^{K_x}b_jx_{t-j} 

x_t=a+\sum_{j=1}^{K_x}b_jx_{t-j}+\sum_{j=1}^{K_y}b_jy_{t-j}

יהיה בעל כוח הסבר רב יותר מזוג המשוואות:

y_t=a+\sum_{j=1}^{K_y}b_jy_{t-j}

x_t=a+\sum_{j=1}^{K_x}b_jx_{t-j}

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

## דוגמה ליישום מבחן גרייג'ר
חיזוי צמיחת תוצר במשק כלכלי. 

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

Y_t=\alpha+\sum_{i=1}^N\beta_iY_{t-i}

כאשר Y הוא שיעור צמיחת התוצר השנתי, עבור N שנים.

הכלכלן כורה את נתוני צמיחת התוצר של המשק ב-20 השנים האחרונות ומגיע למשוואת האמידה הבאה:

y_t=a+b_1y_{t-1}+b_2y_{t-2}+\varepsilon_t

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

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

Y_t=\alpha+\sum_{i=1}^N\beta_iY_{t-i}+\sum_{j=1}^M\gamma_jX_{t-j}

כאשר X הוא שיעור המועסקים השנתי במשך, עבור M שנים.

מכיוון שהכלכלן מעוניין לבדוק את כוחה של סדרת נתוני התעסוקה בניבוי צמיחת התוצר, המודל מצטמצם לכדי המודל הבא:

Y_t=\alpha+\sum_{i=1}^2\beta_iY_{t-i}+\sum_{j=1}^M\gamma_jX_{t-j}

הכלכלן כורה את נתוני התעסוקה של המשק ב-20 השנים האחרונות ומגיע למשוואת האמידה הבאה:

y_t=\tilde{a}+\tilde{b_1}y_{t-1}+\tilde{b_2}y_{t-2}+\tilde{b_3}x_{t_1}+\varepsilon_t

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

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

## ראו גם
* אקונומטריקה
* רגרסיה ליניארית

## לקריאה נוספת

* Granger, C. W. (1969). Investigating causal relations by econometric models and cross-spectral methods. *Econometrica: Journal of the Econometric Society*, 424-438.


קטגוריה:מבחנים סטטיסטיים