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.