兩組率的Meta分析，手把手教Stata操作！

一、問題與數據

Franik S等人擬探討芳香化酶抑制劑治療多囊卵巢綜合征患者，其臨床活產兒率是否高於使用克羅米酚檸檬酸鹽治療的患者。研究者通過檢索資料庫，對最終納入的文獻，提取了如下信息：作者、發表年代、試驗組活產兒數量、試驗組樣本量、對照組活產兒數量、對照組樣本量。

最終納入6篇研究，合計1353名婦女。各研究提取得到的數據如表1。試問：對於多囊卵巢綜合征，使用芳香化酶抑制劑治療的患者，其臨床活產兒率是否高於使用克羅米酚檸檬酸鹽治療的患者？

表1 文獻信息提取表

二、對數據結構的分析

該Meta分析納入了6個研究，每個研究均給出了試驗組活產兒數量、試驗組樣本量、對照組活產兒數量和對照組樣本量4個數據。Meta分析的結局指標為兩組活產兒率比較的OR值，因此屬於兩組率的Meta分析。

三、Stata操作方法

1. 創建數據集並錄入數據

(1)打開數據錄入窗口：Data Data Editor Data Editor (Edit)

關於數據錄入細節，可參考文章：Stata做Meta分析：以診斷試驗準確性為例

各變數名代表含義如下表所示：

（註：非活產兒數量=樣本量-活產兒數量）

2. 兩組率的Meta分析

(1)操作方法：在Command窗口中進行編程。

在「Command」欄中輸入編程語句點擊「Enter」鍵

註：在進行分析前，需要先在stata中安裝meta分析模塊。

在Command欄中輸入語句：

ssc install metan

(2)基本原理

進行兩組率的meta分析，實際上為對原始研究中試驗組和對照組活產兒率的比值進行合并。

(3)具體編程步驟：使用固定效應模型進行兩組率的meta分析。

程序為：metan aevent anonevent bevent bnonevent, or fixed

程序解讀：

兩組率的meta分析的程序為：

Metan 試驗組發生數 試驗組未發生數 對照組發生數 對照組未發生數, 效應指標 使用的模型

metan代表對數據進行meta分析；or代表選擇的效應指標為比值比（病例對照比較研究中，兩組率的meta分析的效應指標均為or）；fixed代表固定效應模型；（若要選擇隨機效應模型，應使用random）。

四、結果解讀

1.首先看異質性檢驗（Heterogeneity）的結果：I2=0%，說明6個原始研究的數據之間不存在異致性，應選用固定效應模型。

2.再看合并的OR值其95%CI（Model Results）：合并的OR=1.805，95%CI：1.396 -2.333。

下圖所示為Meta分析的森林圖。結果解讀同上。

五、撰寫結論

固定效應模型Meta分析結果顯示：對於多囊卵巢綜合征，使用芳香化酶抑制劑治療的患者，其臨床活產兒率高於使用克羅米酚檸檬酸鹽治療的患者，OR值為1. 805（95%CI：1.396 -2.333）。

喜歡這篇文章嗎？立刻分享出去讓更多人知道吧！