一、問題與數(shù)據(jù)
Franik S等人擬探討芳香化酶抑制劑治療多囊卵巢綜合征患者,其臨床活產(chǎn)兒率是否高于使用克羅米酚檸檬酸鹽治療的患者。研究者通過檢索數(shù)據(jù)庫,對最終納入的文獻(xiàn),提取了如下信息:作者、發(fā)表年代、試驗組活產(chǎn)兒數(shù)量、試驗組樣本量、對照組活產(chǎn)兒數(shù)量、對照組樣本量。
最終納入6篇研究,合計1353名婦女。各研究提取得到的數(shù)據(jù)如表1。試問:對于多囊卵巢綜合征,使用芳香化酶抑制劑治療的患者,其臨床活產(chǎn)兒率是否高于使用克羅米酚檸檬酸鹽治療的患者?
表1 文獻(xiàn)信息提取表
二、對數(shù)據(jù)結(jié)構(gòu)的分析
該Meta分析納入了6個研究,每個研究均給出了試驗組活產(chǎn)兒數(shù)量、試驗組樣本量、對照組活產(chǎn)兒數(shù)量和對照組樣本量4個數(shù)據(jù)。Meta分析的結(jié)局指標(biāo)為兩組活產(chǎn)兒率比較的OR值,因此屬于兩組率的Meta分析。
三、Stata操作方法
1. 創(chuàng)建數(shù)據(jù)集并錄入數(shù)據(jù)
(1) 打開數(shù)據(jù)錄入窗口:Data → Data Editor → Data Editor (Edit)
關(guān)于數(shù)據(jù)錄入細(xì)節(jié),可參考文章:Stata做Meta分析:以診斷試驗準(zhǔn)確性為例
(2) 錄入數(shù)據(jù)后,點擊右上角關(guān)閉按鈕
各變量名代表含義如下表所示:
(注:非活產(chǎn)兒數(shù)量=樣本量-活產(chǎn)兒數(shù)量)
2. 兩組率的Meta分析
(1) 操作方法:在Command窗口中進(jìn)行編程。
在“Command”欄中輸入編程語句→點擊“Enter”鍵
注:在進(jìn)行分析前,需要先在stata中安裝meta分析模塊。
在Command欄中輸入語句:
ssc install metan
(2) 基本原理
進(jìn)行兩組率的meta分析,實際上為對原始研究中試驗組和對照組活產(chǎn)兒率的比值進(jìn)行合并。
(3) 具體編程步驟:使用固定效應(yīng)模型進(jìn)行兩組率的meta分析。
程序為:metan aevent anonevent bevent bnonevent, or fixed
程序解讀:
兩組率的meta分析的程序為:
Metan 試驗組發(fā)生數(shù) 試驗組未發(fā)生數(shù) 對照組發(fā)生數(shù) 對照組未發(fā)生數(shù), 效應(yīng)指標(biāo) 使用的模型
metan代表對數(shù)據(jù)進(jìn)行meta分析;or代表選擇的效應(yīng)指標(biāo)為比值比(病例對照比較研究中,兩組率的meta分析的效應(yīng)指標(biāo)均為or);fixed代表固定效應(yīng)模型;(若要選擇隨機(jī)效應(yīng)模型,應(yīng)使用random)。
四、結(jié)果解讀
1. 首先看異質(zhì)性檢驗(Heterogeneity)的結(jié)果:I2=0%,說明6個原始研究的數(shù)據(jù)之間不存在異致性,應(yīng)選用固定效應(yīng)模型。
2. 再看合并的OR值其95%CI(Model Results):合并的OR=1.805,95%CI:1.396 -2.333。
下圖所示為Meta分析的森林圖。結(jié)果解讀同上。
五、撰寫結(jié)論
固定效應(yīng)模型Meta分析結(jié)果顯示:對于多囊卵巢綜合征,使用芳香化酶抑制劑治療的患者,其臨床活產(chǎn)兒率高于使用克羅米酚檸檬酸鹽治療的患者,OR值為1. 805(95%CI:1.396 -2.333)。
本文案例文獻(xiàn)來源:
聯(lián)系客服