Categorical-Data

如何解釋 Cochran-Mantel-Haenszel 檢驗?

  • February 28, 2014

我正在測試由 C 分層的兩個變量 A 和 B 的獨立性。A 和 B 是二元變量,C 是分類變量(5 個值)。運行 Fisher 對 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

最後,使用 A、B 和 C 運行Cochran-Mantel-Haenszel (CMH) 測試,我得到:

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

CMH 檢驗的結果表明,A 和 B在每個層次上**並不獨立(p < 0.05);**但是,大多數層內測試都不顯著,這表明我們沒有足夠的證據來否定 A 和 B 在每個層都是獨立的。

那麼,什麼結論是正確的呢?鑑於這些結果,如何報告結論?C是否可以被認為是一個混雜變量?

編輯:我對零假設進行了 Breslow-Day 檢驗,即各層的優勢比相同,p 值為 0.1424。

第一個測試告訴您,忽略C的 A 和 B 之間的優勢比與1 不同。查看分層分析有助於您確定忽略 C 是否可以。

CMH 測試告訴您 A 和 B 之間的優勢比(針對 C 進行調整)與 1 不同。它返回特定層的優勢比的加權平均值,所以如果這些是在某些階層和在其他情況下,它們可能會抵消並錯誤地告訴您 A 和 B 之間沒有關聯。因此,我們必須測試假設 C 的所有級別的優勢比(在總體級別)相等是否合理。交互作用的 Breslow-Day 檢驗正是這樣做的,假設所有層都具有相同的優勢比,而優勢比不必等於 1。該測試在 EpiR R 包中實現。Breslow-Day p 值為 0.14 意味著我們可以做出這個假設,因此調整後的優勢比是合理的。

但這並不能幫助我們在 CMH 和 Fisher 精確(或 Pearson) 測試。如果 Breslow-Day 檢驗顯著,則需要報告特定層的優勢比。既然不是,那你就需要問是否需要針對 C 進行調整。C 是否“混淆”了 A 和 B 之間的關聯?我學到的啟發式方法(不是統計測試)是檢查未調整和調整後的優勢比之間的比例差異是否超過 10%。這裡,所以CMH是合適的。

引用自:https://stats.stackexchange.com/questions/88315

comments powered by Disqus