הצגת המקור של הדף "מבחן קוקראן-מנטל-הנזל"

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

וריאציות שונות של המבחן הוצעו על ידי [[ויליאם קוקראן]], נתן מנטל וויליאם הנזל{{הערה|{{צ-מאמר|
מחבר=Cochran, W. G.|
שם=Some methods for strengthening the common chi-square tests|
כתב עת=Biometrics|
כרך=10|
סדרה=4|
עמ=417-451|
שנת הוצאה=1954|
doi=10.2307/3001616
}}}}{{הערה|{{צ-מאמר|
מחבר=Mantel, N. and Haensze, W.|
שם=Statistical aspects of the analysis of data from retrospective studies of disease|
כתב עת=Journal of the national cancer institute|
כרך=22|
סדרה=4|
עמ=719-748|
שנת הוצאה=1959|
doi=10.1093/jnci/22.4.7191|
קישור=http://www.med.mcgill.ca/epidemiology/hanley/c634/stratified/Mantel_Haenszel_1.pdf
}}}}{{הערה|{{צ-מאמר|
מחבר=Mantel, N.|
שם=Chi-square tests with one degree of freedom; extensions of the Mantel-Haenszel procedure|
כתב עת=Journal of the American Statistical Association|
כרך=58|
סדרה=303|
עמ=690-700|
שנת הוצאה=1963|
doi=10.1080/01621459.1963.10500879|
קישור=http://www.med.mcgill.ca/epidemiology/hanley/c634/stratified/Mantel_Haenszel_1.pdf
}}}}.

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

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

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

==הגדרה==
===מוטיבציה===
בניסוי קליני שנערך בבית חולים מסוים, חלק מהחולים קיבלו טיפול נסיוני למחלתם, בעוד חלק אחר קיבלו את הטיפול המקובל. עבור כל חולה, נבדק האם מצבו השתפר לאחר שבוע של טיפול. ניתן לרכז את הנתונים בטבלה:
{| class="wikitable"
|סטטוס
|שיפור
|ללא שיפור
|סך הכל
|-
|טיפול ניסיוני
|A
|B
|A+B
|-
|טיפול מקובל
|C
|D
|C+D
|-
|סך הכל
|A+C
|B+D
|N
|}

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

[[יחס הסיכויים]] הוא לכן <math>OR= \frac {A \cdot D} {B \cdot C}</math>.

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

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

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

<math>OR_1= \frac {A_1 \cdot D_1} {B_1 \cdot C_1}</math>,

<math>OR_2= \frac {A_2 \cdot D_2} {B_2 \cdot C_1}</math>.

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

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

מבחן קוקראן-מנטל-הנזל מוגדר על ידי<ref>
{{צ-ספר|
מחבר=Agresti, Alan|
שם=An Introduction to Categorical Data Analysis, 2nd edition|
מקום הוצאה=Hoboken, New Jersey|
מו"ל=JohnWiley & Sons, Inc.|
שנת הוצאה=2007|
עמ=114-115|
פרק=4.3.4|
ISBN=978-0-471-22618-5}}
</ref>:

<math>CMH = {{\sum_{j=1}^J {{A_j \cdot D_j} \over N_j}} \over {\sum_{j=1}^J {{B_j \cdot C_j} \over N_j}}},</math>.

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

==דוגמה: אפליה על רקע מגדרי באוניברסיטת ברקלי==
[[אוניברסיטת קליפורניה בברקלי|באוניברסיטת ברקלי]] בקליפורניה, נתוני הקבלה לסתיו 1973 הראו כי לגברים יש סיכויים גבוהים יותר להתקבל ללימודים בהשוואה לנשים{{הערה|[http://homepage.stat.uiowa.edu/~mbognar/1030/Bickel-Berkeley.pdf .Sex Bias in Graduate Admissions: Data from Berkeley, Bickel et al., Science Vol. 187, 1975]|כיוון=שמאל}}. מתוך 8442 גברים שביקשו להתקבל ללימודים מתקדמים, התקבלו בערך 44% מהפונים. לעומת זאת, מתוך 4321 נשים שביקשו להתקבל ללימודים, התקבלו רק 35%. לא ניתן להסביר פער זה כצרוף מקרים סטטיסטי (במילים אחרות, השערת אי התלות נדחתה).

ניתן לראות זאת גם בנתוני הקבלה של ששת החוגים הגדולים ביותר:
{| class="wikitable"
|
|התקבלו
|נדחו
|סך הכל
|-
|גברים
|1198
|1493
|2691
|-
|נשים
|557
|1278
|1835
|-
|סך הכל
|1755
|2771
|4526
|}
יחס הסיכויים מראה בבירור כי שיעור הקבלה של הגברים גבוה בהרבה מזה של הנשים: <math>OR= \frac {1198 \cdot 1278} {1493\cdot 557} = 1.84</math>

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

{| class="wikitable"
|חוג
|מספר הגברים שהתקבלו
|מספר הגברים שנדחו
|מספר הנשים שהתקבלו
|מספר הנשים שנדחו
|סך הכל
|-
|1
|512
|313
|89
|19
|933
|-
|2
|353
|207
|17
|8
|585
|-
|3
|120
|205
|202
|391
|918
|-
|4
|138
|279
|131
|244
|792
|-
|5
|53
|138
|94
|299
|584
|-
|6
|22
|351
|24
|317
|714
|-
|סך הכל
|1198
|1493
|557
|1278
|4526
|}

ערכו של מבחן קוקראן-מנטל-הנזל הוא <math>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</math>. חישוב [[ערך-p|ערך ה-p]] מעלה כי <math>p\mbox{-value}=0.2323</math>. בהנחה שקבענו מראש כי רמת המובהקות של המבחן היא <math>\alpha=0.05</math> לא נדחה את השערת האפס ונסיק כי לא נצפו הבדלים בהעדפות בין גברים ונשים.

==הערות שוליים==
{{ltr|
{{הערות שוליים}}
}}

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