จะตีความการทดสอบ Cochran-Mantel-Haenszel ได้อย่างไร


10

ฉันกำลังทดสอบความเป็นอิสระของตัวแปรสองตัวคือ A และ B แบ่งเป็นชั้น ๆ โดย C. A และ B เป็นตัวแปรไบนารีและ C คือหมวดหมู่ (5 ค่า) ทำการทดสอบที่แม่นยำของฟิชเชอร์สำหรับ A และ B (ชั้นทั้งหมดรวมกัน) ฉันได้รับ:

##          (B)
##      (A) FALSE TRUE
##    FALSE  1841   85
##    TRUE    915   74

OR: 1.75 (1.25 --  2.44), p = 0.0007 *

โดยที่ OR เป็นอัตราต่อรอง (ประมาณและช่วงความมั่นใจ 95%) และ*หมายความว่า p <0.05

ทำแบบทดสอบเดียวกันสำหรับแต่ละสแตรท (C) ฉันได้รับ:

C=1, OR: 2.31 (0.78 --  6.13), p = 0.0815
C=2, OR: 2.75 (1.21 --  6.15), p = 0.0088 *
C=3, OR: 0.94 (0.50 --  1.74), p = 0.8839
C=4, OR: 1.48 (0.77 --  2.89), p = 0.2196
C=5, OR: 3.38 (0.62 -- 34.11), p = 0.1731

ในที่สุดก็ทำการทดสอบCochran-Mantel-Haenszel (CMH) โดยใช้ A, B และ C ฉันจะได้รับ:

OR: 1.56 (1.12 --  2.18), p = 0.0089 *

ผลจากการทดสอบ CMH แสดงให้เห็นว่า A และ B ไม่เป็นอิสระในแต่ละชั้น (p <0.05) อย่างไรก็ตามการทดสอบภายในสตราตัมส่วนใหญ่ไม่มีนัยสำคัญซึ่งจะแนะนำว่าเราไม่มีหลักฐานเพียงพอที่จะละทิ้งว่า A และ B มีความเป็นอิสระในแต่ละสแตรท

ดังนั้นข้อสรุปที่ถูกต้องคืออะไร? จะรายงานข้อสรุปอย่างไรเมื่อได้ผลลัพธ์เหล่านั้น? C ถือเป็นตัวแปรที่รบกวนหรือไม่?

แก้ไข: ฉันทำการทดสอบ Breslow-Day สำหรับสมมติฐานว่างว่าอัตราต่อรองนั้นเหมือนกันทั่วทั้งชั้นและค่า p คือ 0.1424


3
คุณไม่ได้ทำการทดสอบ Cochran-Mantel-Haenszel อย่างแม่นยำเพราะหลักฐานของอัตราต่อรองที่แตกต่างจากที่อาจจะอ่อนแอสำหรับแต่ละชั้นถือว่าเป็นรายบุคคล แต่แข็งแรงสำหรับการพิจารณาร่วมกัน?
Scortchi - Reinstate Monica

ฉันแสดง CMH เพราะฉันต้องการคำตอบเดียวที่ครบวงจรและฉันต้องการให้แน่ใจว่าผลที่สังเกตระหว่าง A และ B ไม่ได้เกิดจาก C. ฉันอยู่ในเส้นทางที่ถูกต้องหรือไม่? ฉันควรรายงานสถิติของแต่ละชั้นหรือไม่
rodrigorgs

คำตอบ:


10

การทดสอบครั้งแรกจะบอกคุณว่าอัตราต่อรองระหว่าง A และ B โดยไม่สนใจ C นั้นแตกต่างจาก 1 การดูการวิเคราะห์แบบแบ่งชั้นจะช่วยให้คุณตัดสินใจได้ว่าการเพิกเฉยต่อ C นั้นถูกต้องหรือไม่

<1>1ในคนอื่น ๆ พวกเขาสามารถยกเลิกและบอกคุณอย่างผิด ๆ ว่าไม่มีการเชื่อมโยงระหว่าง A และ B ดังนั้นเราต้องทดสอบว่ามีเหตุผลหรือไม่ที่จะสมมติว่าอัตราส่วนอัตราต่อรองเท่ากัน (ในระดับประชากร) ในทุกระดับของ C. การทดสอบการมีปฏิสัมพันธ์ระหว่างวันของ Breslow ทำอย่างนี้โดยมีสมมติฐานว่างว่าชั้นทั้งหมดมีอัตราส่วนอัตราต่อรองเท่ากันซึ่งไม่จำเป็นต้องเท่ากับหนึ่ง การทดสอบนี้ถูกนำไปใช้ในแพ็คเกจ EpiR R ค่า Breslow-Day ที่ 0.14 หมายถึงเราสามารถตั้งสมมติฐานนี้ได้ดังนั้นอัตราต่อรองที่ปรับแล้วนั้นถูกต้องตามกฎหมาย

χ21.751.561.75=0.108


ฉันแก้ไขคำถามของฉันเพื่อเพิ่มผลลัพธ์จากการทดสอบ Breslow-Day (0.14) ดังนั้นฉันสามารถพูดได้ว่ามันมีเหตุผลที่จะสมมติว่าอัตราต่อรองที่เท่ากัน? ในกรณีนี้ฉันควรรายงานอัตราต่อรองของฟิชเชอร์หรือ CMH หรือไม่
rodrigorgs

2
สมมติฐานว่างของ Breslow-Day คือ "อัตราต่อรองที่เป็นเนื้อเดียวกันทั่วทั้งชั้น" เนื่องจากค่า ap> 0.05 ไม่ได้หมายความว่า null เป็นจริงคุณจึงไม่สามารถสันนิษฐานได้ว่าอัตราต่อรองเท่ากัน
Michael M

@MichaelMayer: ฉันคิดว่าคุณหมายถึงการพูดว่า "สมมติฐานของอัตราต่อรองที่เป็นเนื้อเดียวกันไม่น่าอดสู แต่คุณไม่ควรสับสนที่จะปฏิเสธโมฆะด้วยการพิสูจน์โมฆะ"
Scortchi - Reinstate Monica

@ vafisher: มีสิ่งหนึ่งที่ผิด - ประโยคที่ 3: การทดสอบของชาวประมงยังคงไม่เหมาะสมเมื่ออัตราต่อรองที่แตกต่างกันในระดับที่แตกต่างกันของ C.
Scortchi - Reinstate Monica

@Scortchi: จุดดี!
vafisher
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.