二分類變量的meta分析,兩組發(fā)生率為0,沒法分析,那100%呢?
二分類變量的meta分析,兩組的發(fā)生率都是100%,能不能做meta分析?

以下面的數(shù)據(jù)為例,我們用Stata和RevMan驗證一下,兩組發(fā)生率為100%時,用RR為效應值,是否能得到結(jié)果。

首先,看看Stata11的結(jié)果,隨機效應模型(I-V法)的結(jié)果是RR (95%CI) = 0.60 (0.50, 0.73)。然而,兩組發(fā)生率為100%的A5,并沒有結(jié)果,在森林圖中被排除(Excluded)。

接下來,是Stata 14的隨機效應模型I-V法的結(jié)果,RR (95%CI) = 0.62 (0.52, 0.75)。A5被納入分析,該研究的RR (95%CI) = 1.00 (0.94, 1.06)。

為什么Stata11中,A5會被排除呢?
因為二分類變量的meta分析,Stata使用的數(shù)據(jù)是“a, b, c, d”,而A5對應的數(shù)據(jù)是“35, 0, 30, 0”,同樣符合兩組為0,自動排除的規(guī)則。
當然,這是一個Bug,因為RR的計算公式是:
RR = (a/(a+b))/(c/(c+d))
將數(shù)據(jù)代入公式,RR = [35/(35+0)]/[30/(30+0)] = 1
因此,Stata后來更新了計算的代碼,更新后的版本可以將兩組發(fā)生率為100%的研究納入RR的森林圖。
我們再用RevMan 5.3.5驗證一下,隨機效應模型I-V法的結(jié)果為RR (95%CI) = 0.62 (0.52, 0.75),與Stata14一致。

二分類變量的meta分析,以RR值為合并效應值,兩組發(fā)生率為100%時,你知道怎么分析了嗎?
以OR值為合并效應值呢?歡迎留言分享答案。
