Meta簡明教程目錄

5.Meta分析數(shù)據(jù)提取| Meta簡明教程（5）

6.一文學(xué)會revman軟件| Meta簡明教程（6）

Meta簡明教程（7）

上一期介紹了Revman 軟件對二分類數(shù)據(jù)、連續(xù)型數(shù)據(jù)、 診斷性試驗(yàn)數(shù)據(jù)、 生存-時間數(shù)據(jù)進(jìn)行meta分析，本期將利用Stata對以上數(shù)據(jù)進(jìn)行meta分析。大家可以到本公眾號下載Stata軟件（重磅推薦：分類最全的統(tǒng)計(jì)分析相關(guān)軟件，了解一下？請關(guān)注、收藏以備用）

Stata12.0 界面

一、二分類數(shù)據(jù)分析

數(shù)據(jù)形式

例：研究阿司匹林（aspirin）預(yù)防心肌梗死（MI）7個臨床隨機(jī)對照試驗(yàn)，觀察死亡率，數(shù)據(jù)提取如下：

操作步驟

1.構(gòu)建數(shù)據(jù)

1)啟動Stata 12.0 軟件后，可以直接點(diǎn)擊工具欄中DataEditor (edit)按鈕。

也可在在菜單欄中點(diǎn)擊Data→Data Editor→ DataEditor (edit)，出現(xiàn)以下界面。

2）點(diǎn)擊變量名位置，依次輸入研究名稱（research），阿司匹林組死亡數(shù)(a)，阿司匹林組存活數(shù)(b)，安慰劑組死亡數(shù)(c)，安慰劑組存活數(shù)(d)

3）錄入數(shù)據(jù)：在變量值區(qū)域輸入數(shù)據(jù)

2. 數(shù)據(jù)分析

1）導(dǎo)入meta模塊：

在Command窗口中進(jìn)行編程，首先需要在Stata中安裝meta模塊:在Command窗口輸入“ssc install metan”，選中點(diǎn)回車。結(jié)果窗口中出現(xiàn)下面的結(jié)果，說明已經(jīng)安裝了meta模塊。

2）輸入meta分析代碼：

在Command窗口輸入 “Command窗口輸入 “metan a b c d, or fixed”,點(diǎn)回車，完成結(jié)果分析。a,b,c,d分別代表阿司匹林組死亡數(shù)，阿司匹林組存活數(shù)，安慰劑組死亡數(shù)，安慰劑組存活數(shù)，or代表效應(yīng)指標(biāo)為OR值，fixed代表選用固定效應(yīng)模型。相應(yīng)的可以用rr 及random替換代碼。

結(jié)果

1．合并效應(yīng)

2.森林圖

二、連續(xù)型數(shù)據(jù)分析

數(shù)據(jù)形式

例：研究小劑量潑尼松和安慰劑治療類風(fēng)濕性關(guān)節(jié)炎的7個臨床隨機(jī)對照試驗(yàn)，觀察類風(fēng)濕性關(guān)節(jié)炎患者的關(guān)節(jié)壓痛指數(shù)，數(shù)據(jù)提取如下：

操作步驟

1.構(gòu)建數(shù)據(jù)

數(shù)據(jù)的構(gòu)建和二分類的步驟相同，需要設(shè)置的變量名有研究名稱（research）、試驗(yàn)組樣本量(tn)，試驗(yàn)組均數(shù)(tmean)，試驗(yàn)組標(biāo)準(zhǔn)差(tsd)，對照組均數(shù)(cmean)，對照組標(biāo)準(zhǔn)差(csd)。

2. 數(shù)據(jù)分析

1）導(dǎo)入meta模塊（與二分類數(shù)據(jù)一致，如果已經(jīng)導(dǎo)入可跳過此步驟）

2）輸入meta分析代碼：

在Command窗口輸入 “metan tn tmean tsd cn cmean csd, fixed”,點(diǎn)回車，完成結(jié)果分析。tn tmean tsd分別代表試驗(yàn)組樣本量、均數(shù)、標(biāo)準(zhǔn)差，cncmean csd分別代表對照組樣本量、均數(shù)、標(biāo)準(zhǔn)差，fixed代表選用固定效應(yīng)模型。如果計(jì)算合并均差（WMD）需要在fixed后加nostandard,如果不加，默認(rèn)合并標(biāo)化均差（SMD）。

結(jié)果

1．合并效應(yīng)

2. 森林圖

三、診斷性試驗(yàn)數(shù)據(jù)的分析

數(shù)據(jù)形式

例：研究ProGRP和NSE對小細(xì)胞肺癌診斷10個診斷試驗(yàn)，觀察ProGRP對小細(xì)胞肺癌的診斷價值，數(shù)據(jù)提取如下：

操作步驟及結(jié)果

1.構(gòu)建數(shù)據(jù)

數(shù)據(jù)的構(gòu)建和二分類的步驟相同，需要設(shè)置的變量名有研究名稱（research）、真陽性（tp）、假陽性(fp)、假陰性(fn)、真陰性(tn)。

2.數(shù)據(jù)分析

1）導(dǎo)入診斷性試驗(yàn)meta分析模塊：

在Command窗口輸入 “sscinstall midas” ，點(diǎn)回車，結(jié)果窗口顯示

2）效應(yīng)值合并：

在Command窗口輸入 “midas tp fp fn tn, res(sum)”,點(diǎn)回車，完成結(jié)果分析。tp fp fn tn真陽性、假陽性、假陰性、真陰性例數(shù)，res(sum)代表計(jì)算合并點(diǎn)值，包括靈敏度、特異度、陽性似然比、陰性似然比、診斷比值比。

3）森林圖繪制

在Command窗口輸入 “midas tp fp fn tn, ford fors bfor(dss)”,點(diǎn)回車, 完成結(jié)果分析。tp fp fn tn真陽性、假陽性、假陰性、真陰性例數(shù). Ford代表研究的置信區(qū)間，fors代表異質(zhì)性結(jié)果的分析，bfor(dss)列出每一項(xiàng)研究。

4）SROC曲線繪制

在Command窗口輸入 “midas tp fp fn tn, plot sroc(both)”,點(diǎn)回車, 完成結(jié)果分析。Plot sroc(both)代表繪制SROC曲線。

四、生存-時間數(shù)據(jù)分析

數(shù)據(jù)形式

例：結(jié)直腸癌肝轉(zhuǎn)移生存風(fēng)險的3個臨床研究，觀察肝轉(zhuǎn)移對結(jié)直腸癌生存風(fēng)險比的影響，數(shù)據(jù)提取如下：

操作步驟

1.構(gòu)建數(shù)據(jù)

數(shù)據(jù)的構(gòu)建和二分類的步驟相同，需要設(shè)置的變量名有研究（research）,風(fēng)險比（hr）,風(fēng)險比置信區(qū)間下限（lci）,風(fēng)險比置信區(qū)間上限（uci）。

2.數(shù)據(jù)分析

1）導(dǎo)入meta模塊（與二分類數(shù)據(jù)一致，如果已經(jīng)導(dǎo)入可跳過此步驟）

2）輸入meta分析代碼：

在Command窗口輸入 “metan hr lci uci, effect(HR) fixed”,點(diǎn)回車，完成結(jié)果分析。hr、lci、uci分別代表各研究風(fēng)險比,風(fēng)險比置信區(qū)間下限,風(fēng)險比置信區(qū)間上限，effect(HR)代表效應(yīng)量，fixed代表選用固定效應(yīng)模型。

結(jié)果

1．合并效應(yīng)

2.森林圖

五、總結(jié)

相比于Revman，Stata操作步驟更為簡單，而且靈活多樣，功能也比Revman強(qiáng)大，但Stata主要是代碼界面，需要進(jìn)行簡單編程，上手比Revman難。

掃碼關(guān)注我們

了解更多醫(yī)學(xué)論文寫作方法、真實(shí)世界研究統(tǒng)計(jì)分析方法,請關(guān)注公眾號：“醫(yī)學(xué)論文與統(tǒng)計(jì)分析”或“medicalstat”