江西中醫(yī)藥大學(xué) 馬艷嬌 陳曉凡

Stata軟件是一款統(tǒng)計分析軟件，不但功能強大，而且具有極為全面的Meta系統(tǒng)分析作用，所以一直受到影響力、認可度極高的國外雜志認可。以Meta分析RevMan與Stata軟件利弊、Meta中Stata軟件運行與分析、Meta分析中常用命令、Stata界面窗口功能等為切入點來描述Stata軟件。

循證醫(yī)學(xué)證據(jù)數(shù)量眾多，其中主要方法之一Meta分析，可以數(shù)個研究結(jié)果整合為基礎(chǔ)將量化聯(lián)系或平均效果提供出來，通過樣本量加大這一基本途徑來實現(xiàn)結(jié)論把握度提升，就此將研究結(jié)果偏差有效化解是其典型特征和基本優(yōu)勢。隨機效應(yīng)模型或者固定效應(yīng)模型是Meta采用兩種基本模型，因為真實效應(yīng)相同（納入研究間）是固定效應(yīng)模型假設(shè)，因此其應(yīng)用會受到某種約束；但隨機效應(yīng)模型極為普遍，其應(yīng)用最多[1,2]。

Stata軟件實質(zhì)上屬于整合性統(tǒng)計分析軟件，能將圖表繪制、統(tǒng)計分析、數(shù)據(jù)管理功能提供出來，所以其普及運用于生物統(tǒng)計學(xué)、經(jīng)濟學(xué)及社會學(xué)。美國Computer Resource Center（計算機資源中心）為其原始研發(fā)單位，后為Stata公司收購并成為其旗下產(chǎn)品[3]。1.0版于1985年發(fā)布，自此，十種主要版本陸續(xù)發(fā)布，且自4.0版開始以來全面迎來Win時代，在擴展、升級、創(chuàng)新持續(xù)進行過程中逐步完善軟件功能，同時具有下列各項典型特征，即簡單、靈活、操作方便，程序語言、矩陣計算軟件、繪圖軟件、統(tǒng)計分析軟件及數(shù)據(jù)管理軟件[4]。2015年，可在Stata 14.0版本中直接進行Meta分析，它可分步提供詳細的結(jié)果和理想的圖形，對那些對于統(tǒng)計學(xué)了解不多的研究者，Stata 14.0應(yīng)該是更好的選擇[5,6]。

1 Stata界面窗口功能

對Stata軟件圖標(biāo)執(zhí)行單擊操作即可完成Stata界面登錄。菜單欄包括“File”“Edit”“Data”“Graph ics”“Statistics”“User”“Windows”和“Help”等（Win版軟件通用狀態(tài)欄、工具欄、菜單欄除外），而Stata運行包括四個窗口，分別為結(jié)果窗口（Stata Results）、命令窗口（Stata Command）、命令回顧窗口（Review）和變量名窗口（Variables）。

(1)結(jié)果窗口（Stata Results）：界面右上部即為結(jié)果窗口，此處主要展示出錯信息、執(zhí)行結(jié)果和所執(zhí)行命令等軟件運行信息。

(2)命令窗口（Stata Command）：結(jié)果窗口右下部即為命令窗口，和命令行（DOS中）比較相似，需執(zhí)行指令在此鍵入，回車完成，執(zhí)行開始，且結(jié)果窗口內(nèi)會顯示有關(guān)結(jié)果。

(3)命令回顧窗口（Review）：界面左上部即為命令回顧窗口，該窗口將會依次列出所有命令（已經(jīng)執(zhí)行過），執(zhí)行單擊操作即可自動復(fù)制命令至命令窗口；若有重復(fù)執(zhí)行需要，只要在Review窗口中用鼠標(biāo)對有關(guān)命令執(zhí)行雙擊操作就能達到目標(biāo)。

(4)變量名窗口（Variables）：界面左下部即為變量名窗口，各種變量名稱。

2 異質(zhì)性評價分析

用于各類效應(yīng)量變異程度描述，能將研究結(jié)果的差異性（不包括可預(yù)見偶然機會）表達出來。方法學(xué)異質(zhì)性、隨機誤差、臨床異質(zhì)性均為異質(zhì)性的主要來源，H統(tǒng)計量、I2統(tǒng)計量和Q值統(tǒng)計量等檢驗法均為統(tǒng)計學(xué)檢驗評價異質(zhì)性的基本方法，敏感性分析、Meta回歸和亞組分析則為異質(zhì)性探索技術(shù)。

3 Stata軟件在Meta中的運行及分析

Metan是Meta Analysis命令，可以合并二分類及連續(xù)數(shù)據(jù)；且Meta分析效應(yīng)量、可信區(qū)間或相應(yīng)標(biāo)準誤差；且可將P值、I2統(tǒng)計量獲得并運用于檢驗異質(zhì)性，并將Forest Plot（森林圖）同步給出。其命令行格式為：Metan變量名[, 選擇項]。

3.1 二分類變量

Meta命令有四個變量。分別是a、b、c、d，也就是試驗組發(fā)生事件數(shù)、試驗組沒有發(fā)生事件數(shù)、對照組發(fā)生事件數(shù)、對照組沒有發(fā)生事件數(shù)。導(dǎo)入Excel，檢驗數(shù)據(jù)是否導(dǎo)入，進行數(shù)據(jù)分析：User窗口：(1)常規(guī)Meta分析（森林圖）；(2)Meta回歸；(3)Meta回歸圖；(4)常規(guī)Meta亞組分析；(5)Meta漏斗圖；(6)發(fā)表偏倚；(7)Meta敏感分析。

3.2 連續(xù)性變量

Meta命令有對照組的樣本量、對照組標(biāo)準差、對照組均數(shù)，試驗組的樣本量、試驗組標(biāo)準差、試驗組均數(shù)6個變量，如Metan n1 m1 sd1 n2 m2 sd2。導(dǎo)入Excel，檢驗數(shù)據(jù)是否導(dǎo)入，進行數(shù)據(jù)分析：User窗口：(1)常規(guī)Meta分析（森林圖）；(2)Meta回歸；(3)Meta回歸圖；(4)常規(guī)Meta亞組分析；(5)Meta漏斗圖；(6)發(fā)表偏倚；(7)Meta敏感分析。

4 實例分析

(1)二分類變量：實例演示基于系統(tǒng)數(shù)據(jù)（Stata軟件Metan命令）展開。例如：納入樣本年齡（3、2、1分別代表老、中、青年組）用Age表示；納入研究用Study表示；納入病例的時期（1表示早期，2表示中期和3表示晚期）用TMN表示；試驗組發(fā)生數(shù)和未發(fā)生數(shù)分別用a和b表示；對照組發(fā)生數(shù)用c表示；對照組沒有發(fā)生數(shù)用d表示(如表1所示)。

表1 Meta中納入的文獻及相關(guān)數(shù)據(jù)Tab.1 Included literature and related data in Meta

常規(guī)Meta分析(森林圖)：數(shù)據(jù)發(fā)現(xiàn)，I2=59.5%，p=0.003，具有異質(zhì)性，隨機效應(yīng)模型被選中，0.97（0.78,1.22）為可信區(qū)間，無效線相交于菱形（代表無統(tǒng)計學(xué)意義），因此回歸與亞組分析應(yīng)展開（如圖1所示）。

圖1 Meta分析森林圖Fig.1 Forest map of Meta-analysis

Meta回歸：異質(zhì)性與年齡存在相關(guān)性（如表2所示）。

表2 Meta中回歸數(shù)據(jù)Tab.2 Regression data in Meta

Meta回歸圖：異質(zhì)性與年齡存在相關(guān)性(如圖2所示)。

圖2 Meta分析回歸圖Fig.2 Regression diagram of Meta-analysis

常規(guī)Meta亞組分析：中年組沒有任何意義，老年組與青年組有意義(如圖3所示)。

圖3 Meta分析亞組圖Fig.3 Subsets of Meta-analysis

Meta漏斗圖：文獻并不對稱，所以具有發(fā)表偏倚(如圖4所示)。

圖4 Meta分析漏斗圖Fig.4 Funnel diagram of Meta analysis

發(fā)表偏倚：文獻并不對稱，具有發(fā)表偏倚(如圖5所示)。(2)連續(xù)性變量：與二分類變量相似。

圖5 Meta分析發(fā)表偏倚圖Fig.5 Publication bias of Meta analysis

5 展望

Meta分析目前更傾向于選擇Stata軟件，其更是高質(zhì)量國外雜志進行Meta系統(tǒng)分析的首選，其能迅速、高效的識別發(fā)表性偏倚（Meta分析），同時可以圖文并茂、生動有趣的完成Stata軟件技術(shù)描述，并且方便、快捷、準確、可靠做出Meta分析的圖[7]。

引用

[1] 孫振球,徐勇勇.醫(yī)學(xué)統(tǒng)計學(xué)[M].北京:人民衛(wèi)生出版社,2005:740-788.

[2] 陳新林,胡月,郎建英,等. RevMan、Stata和R軟件在生存資料Meta分析的應(yīng)用[J].中國循證醫(yī)學(xué)雜志編輯部,2016(6):736-740.

[3] 張?zhí)鞎?Stata和R軟件在生存資料Meta分析中的正確使用[J].中國循證醫(yī)學(xué)雜志,2017(10):1237-1230.

[4] 莫傳偉,陳群,徐志偉.多中心臨床隨機對照試驗四格表資料的綜合分析及Stata軟件應(yīng)用[J].中醫(yī)藥學(xué)刊,2006(24):1013-1015.

[5] 趙景波.累積Meta分析方法及其在臨床醫(yī)學(xué)研究中的應(yīng)用[J].循證醫(yī)學(xué),2002(2):167-171.

[6] 王靜,莫傳偉,陳群,等.應(yīng)用Stata軟件的Meta分析完成醫(yī)學(xué)科研定量評價[J].中華中醫(yī)藥學(xué)刊,2008(5):947-949.

[7] 康德英,洪旗,劉關(guān)鍵,等.Meta分析中發(fā)表性偏倚的識別與處理[J].中國循證醫(yī)學(xué)雜志,2003(3):45-49.