|
言歸正傳，下面筆者將以實(shí)例演示的形式展示如何使用Stata軟件畫出一張森林圖。本文所舉例子的數(shù)據(jù)是來自于若干個鏈激酶干預(yù)急性心肌梗死的隨機(jī)對照試驗(yàn)，試驗(yàn)組是鏈激酶+標(biāo)準(zhǔn)治療，對照組為標(biāo)準(zhǔn)治療。觀察終點(diǎn)為：心肌梗死所致病死率。 第一步，在Excel中整理各納入研究的原始數(shù)據(jù)如下，從左至右每列依次為：試驗(yàn)名稱，研究發(fā)表的年代，實(shí)驗(yàn)組的死亡人數(shù)，實(shí)驗(yàn)組的存活人數(shù)，對照組的死亡人數(shù)，對照組的存活人數(shù)。對于大多數(shù)臨床醫(yī)生來說Excel管理數(shù)據(jù)已經(jīng)綽綽有余，建議讀者采用Excel作為通用工具管理數(shù)據(jù)，再通過導(dǎo)入或者拷貝的方式把已經(jīng)整理好的數(shù)據(jù)錄入統(tǒng)計(jì)軟件。如下圖所示： 第二步，拷貝上述數(shù)據(jù)并粘貼至Stata軟件，具體操作如所示：拷貝目標(biāo)數(shù)據(jù)→在Stata數(shù)據(jù)編輯窗口，單擊左鍵選中左上角的單元格→點(diǎn)擊右鍵，粘貼。這時(shí)會彈出一個對話框，選擇“treat first rowas variable names”。如下所示。 至此數(shù)據(jù)錄入完畢。數(shù)據(jù)錄入要點(diǎn)，對于二分類數(shù)據(jù)的meta分析，數(shù)據(jù)錄入的順序最好采用如下順序：試驗(yàn)組事件數(shù)、實(shí)驗(yàn)組未發(fā)生事件數(shù)、對照組事件數(shù)、對照組未發(fā)生事件數(shù)。這一順序也是后續(xù)meta分析命令強(qiáng)制要求的順序，如果錄入數(shù)據(jù)就采用這樣的順序，便于后續(xù)操作。 輸入如下命令：
“metan”是上述這段命令的主語，需要提前在stata軟件中安裝metan程序包；“ntdeath ntalivencdeath ncalive”分別是實(shí)驗(yàn)組的死亡人數(shù)，實(shí)驗(yàn)組的存活人數(shù)，對照組的死亡人數(shù)，對照組的存活人數(shù)，這個順序不能改變；“l(fā)cols(trialname year)”表示在森林圖的左邊分別顯示trialname和year兩列變量；“random”指定隨機(jī)效應(yīng)模型；“or”指定效應(yīng)量為Odds Ratio（優(yōu)勢比）；“counts group1(Test) group2(Control)”在森林圖的右邊顯示實(shí)驗(yàn)組與對照組的事件數(shù)與總例數(shù)；“xlabel(0.05,1,10)force”強(qiáng)制森林圖僅顯示0.05,1,10三個橫坐標(biāo)，且森林圖只能在這個范圍內(nèi)顯示，超出這個范圍則用箭頭代替；“boxsca(40)”顯示權(quán)重的方塊大小為40個單位。森林圖如下： 數(shù)字化結(jié)果如下： 第四步，結(jié)果的解讀。 合并的OR，0.782（0.693, 0.884）；異質(zhì)性檢驗(yàn)結(jié)果：p=0.066<0.1，存在異質(zhì)性；異質(zhì)性的程度I2=31.5%，輕度異質(zhì)性。異質(zhì)性可接受。研究間的方差Tau2=0.0168，研究間的變異可接受。如下圖所示。專業(yè)解釋：鏈激酶干預(yù)組較對照組顯著降低因急性心肌梗死病死率，OR=0.782（0.693, 0.884）。
作者| 周支瑞，復(fù)旦大學(xué)腫瘤醫(yī)院放療科在讀博士。AME學(xué)術(shù)沙龍委員，喜讀書，廣交友。業(yè)余時(shí)間擔(dān)任丁香園循證醫(yī)學(xué)討論版版主，自學(xué)臨床流行病學(xué)與循證醫(yī)學(xué)五年余。主要研究方向：惡性腫瘤放射治療與放射生物學(xué)研究、循證醫(yī)學(xué)與Meta分析方法學(xué)研究。目前以第一作者及共同作者發(fā)表SCI論文十余篇，參編學(xué)術(shù)著作兩部。
|
|
|
