Stata軟件在Meta系統分析中的應用*

江西中醫藥大學 馬艷嬌 陳曉凡

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

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

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

1 Stata界面窗口功能

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

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

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

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

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

2 異質性評價分析

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

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

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

3.1 二分類變量

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

3.2 連續性變量

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

4 實例分析

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

表1 Meta中納入的文獻及相關數據Tab.1 Included literature and related data in Meta

常規Meta分析(森林圖)：數據發現，I2=59.5%，p=0.003，具有異質性，隨機效應模型被選中，0.97（0.78,1.22）為可信區間，無效線相交于菱形（代表無統計學意義），因此回歸與亞組分析應展開（如圖1所示）。

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

Meta回歸：異質性與年齡存在相關性（如表2所示）。

表2 Meta中回歸數據Tab.2 Regression data in Meta

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

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

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

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

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

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

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

圖5 Meta分析發表偏倚圖Fig.5 Publication bias of Meta analysis

5 展望

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

引用

[1] 孫振球,徐勇勇.醫學統計學[M].北京:人民衛生出版社,2005:740-788.

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

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

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

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

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

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