近段時(shí)間，有不少小伙伴給我留言，為什么二分類變量的meta分析，得到的OR值是負(fù)數(shù)，這個(gè)結(jié)果怎么解釋？?我們都知道，OR是“Odd ratio”的縮寫，中文名稱是“比值比”，是兩組計(jì)數(shù)資料頻數(shù)分布的比值差異，其取值范圍是0到+∞?？墒?，為什么軟件會(huì)出現(xiàn)OR<0的結(jié)果呢？

?
這個(gè)事情，看上去是軟件的問題，但仔細(xì)一想，問題出在我們自己身上。?先說軟件的問題，現(xiàn)在很多meta分析教程都是基于低版本軟件的，例如Stata 11/12，對于這兩個(gè)版本，二分類變量的meta分析，若使用OR值為效應(yīng)值，計(jì)算得到的是OR值。?以下面這些數(shù)據(jù)為例，用OR (95%CI)評(píng)估兩組死亡風(fēng)險(xiǎn)的差異

?
通過菜單或命令行操作，得到隨機(jī)效應(yīng)模型的合并結(jié)果，森林圖以O(shè)R (95%CI)展示結(jié)果。

?
數(shù)據(jù)窗口，軟件也生成了兩個(gè)臨時(shí)變量_ES和_selogES，也就是OR和selogOR。

?
如果接下來我們要做敏感性分析或發(fā)表偏倚檢驗(yàn)，需要計(jì)算logOR和selogOR，那就要利用以下代碼將_ES轉(zhuǎn)換為logOR：gen logOR=log(_ES)?可是，當(dāng)我們用Stata14.0操作的時(shí)候，出來的數(shù)據(jù)是這樣的！兩個(gè)臨時(shí)變量，_selogOR變成了_seES，但數(shù)據(jù)完全一樣；而_ES的變量名沒變，可數(shù)據(jù)截然不同，甚至出現(xiàn)了負(fù)數(shù)。

?
這個(gè)時(shí)候再進(jìn)行對數(shù)轉(zhuǎn)換，就會(huì)出現(xiàn)部分?jǐn)?shù)據(jù)缺失的情況。為什么有缺失？負(fù)數(shù)沒法進(jìn)行自然對數(shù)轉(zhuǎn)換，這是中學(xué)數(shù)學(xué)的知識(shí)！

?
那是軟件出錯(cuò)了嗎？軟件也沒錯(cuò)，它只是在版本升級(jí)的同時(shí)，做了一些更新，將二分類變量的meta分析結(jié)果以對數(shù)值的形式展示，也就是此時(shí)的_ES實(shí)際上就是11/12版本的logOR。
我們可以驗(yàn)證一下，下面的logOR是否與上圖的_ES (保留小數(shù)點(diǎn)后7位)一模一樣？

?
OK，第一個(gè)情況解決了。如果你是使用Stata14.0完成了二分類變量的meta分析，得到_ES后就不要再進(jìn)行對數(shù)轉(zhuǎn)換了，因?yàn)樗褪荗R的自然對數(shù)。?還有一種更有趣的情況，使用Stata16的新菜單“meta”實(shí)現(xiàn)二分類變量的森林圖，如果菜單選擇有誤，得到的森林圖都是以對數(shù)值顯示結(jié)果。