מבחן קוקראן-מנטל-הנזל

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

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

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

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

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

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

כאשר A הוא מספר החולים שקיבלו טיפול ניסיוני ומצבם השתפר, B הוא מספר החולים שקיבלו טיפול ניסיוני ומצבם לא השתפר, וכן הלאה. כמו כן, N הוא המספר הכולל של החולים שהשתתפו בניסוי.

יחס הסיכויים הוא לכן OR= \frac {A \cdot D} {B \cdot C}.

בהנחה כי אין קשר בין הטיפול ומצב החולה, נצפה כי מספר החולים שמצבם השתפר לאחר קבלת הטיפול הניסיוני הוא E= \frac {(A+B) \cdot (A+C)} {N}.

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

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

OR_1= \frac {A_1 \cdot D_1} {B_1 \cdot C_1},

OR_2= \frac {A_2 \cdot D_2} {B_2 \cdot C_1}.

מבחן קוקראן-מנטל-הנזל משקלל את יחסי הסיכויים של שני בתי החולים לפי מספר החולים בכל אחד מהם. באופן אינטואיטיבי, השקלול צריך להיות ממוצע משוקלל של יחסי הסיכויים משני בתי החולים. אולם, אפשרות אחרת לבדיקת ההשערה כייחס הסיכויים הכולל שווה ל-1 היא לבדק את הערך של S=(A_1-E_1)^2 + (A_2-E_2)^2. ביטוי זה מהווה את הבסיס להגדרת מבחן קוקראן-מנטל-הנזל, והוא מתקבל לאחר שקלול מתאים של S בשונויות של A_1 ו-A_2.

### הגדרה פורמלית
נתון לוח שכיחות תלת־ממדי מסדר 2 \times 2 \times J כאשר לכל j, לוח השכיחות הדו-ממדי הוא \begin{bmatrix} A_j & B_j \\ C_j & D_j \end{bmatrix}. נסמן
N_j = A_j + B_j + C_j + D_j.

מבחן קוקראן-מנטל-הנזל מוגדר על ידי:

CMH = {{\sum_{j=1}^J  \over {\sum_{j=1}^J },.

תחת השערת האפס כי יחס הסיכויים הכולל שווה ל-1, ל-CMH יש התפלגות אסימפטוטית \chi^2 עם דרגת חופש אחת.

## דוגמה: אפליה על רקע מגדרי באוניברסיטת ברקלי
באוניברסיטת ברקלי בקליפורניה, נתוני הקבלה לסתיו 1973 הראו כי לגברים יש סיכויים גבוהים יותר להתקבל ללימודים בהשוואה לנשים. מתוך 8442 גברים שביקשו להתקבל ללימודים מתקדמים, התקבלו בערך 44% מהפונים. לעומת זאת, מתוך 4321 נשים שביקשו להתקבל ללימודים, התקבלו רק 35%. לא ניתן להסביר פער זה כצרוף מקרים סטטיסטי (במילים אחרות, השערת אי התלות נדחתה).

ניתן לראות זאת גם בנתוני הקבלה של ששת החוגים הגדולים ביותר:

יחס הסיכויים מראה בבירור כי שיעור הקבלה של הגברים גבוה בהרבה מזה של הנשים: OR= \frac {1198 \cdot 1278} {1493\cdot 557} = 1.84

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

ערכו של מבחן קוקראן-מנטל-הנזל הוא CMH = \frac { \frac {512 \cdot 19} {933 } + \frac {353 \cdot 8} {585} + \ldots}{ \frac {313 \cdot 89} {933 } + \frac {207 \cdot 17} {585} + \ldots } = 1.4269. חישוב ערך ה-p מעלה כי p\mbox{-value}=0.2323. בהנחה שקבענו מראש כי רמת המובהקות של המבחן היא \alpha=0.05 לא נדחה את השערת האפס ונסיק כי לא נצפו הבדלים בהעדפות בין גברים ונשים.

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

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