R软件gemtc程序包在网状Meta分析中的应用

网状Meta分析

张超－董圣杰ｚ

曾宪涛－，，＋

１．湖北医药学院附属太和医院循证医学中心（湖北十堰４４２０００）；２．苏州大学附属第一医院骨科（江苏苏州２１５００６）；３．湖北医药学院附属太和医院口腔科（湖北十堰“２０００）

摘要Ｒ软件相关程序包对于网状Ｍｅｔａ分析的理论原理大多是基于贝叶斯统计模型，也有少部分使用广义ＩＡＧＳ、０ｐｅｎＢｕＧＳ及ｗｉｎＢｕＧＳ）来执行网状Ｍｅｔａ分析。同时，还可以为ＧｅＭＴＣ软件生成相关的数据存储文件。

其本质是调用基于ＭｃＭｃ法的软件进行网状Ｍｅｔａ分析。本文简要介绍Ｒ软件ｇｅｍｔｃ程序包调用ｏｐｅｎＢｕＧｓ软

件进行网状Ｍｅｔａ分析。

关键词Ｒ软件；网状Ｍｅｔａ分析；ｇｅｍｔｃ程序包

Ａｐｐｌｉｃａｔｉｏｎ

ｏｆ

ＧｅＭＴＣ

ＲＰａｃｋａｇｅｏｆＲ

Ｓｏ觚随ｒｅｆｏｒＮｅ咖ｒｋＭｅｔａ－Ａｎａｌｙｓｉｓ

ＺＨＡＮＧＣｈａ０１，ＤＯＮＧＳｈｅｎｇ－ｊｉｅ２，ＺＥＮＧ

Ｘｉａｎ－ｔａ０１，３’

Ｊ．＆ｎ把咖ｒ西砘Ｈ卯一Ｂ船ｅｄ胁成ｃｆｎｅｎ积Ｃｆｆ”ｆ叩ｆＲ酷阳砌，诒他ＨＤ妒ｆ矾Ｈ“搬咖眈圩砂矿胁成咖Ｂｓ＾咖ｎ４４２０００，∞ｆｎ研

２．Ｄ叩口ｒ所ｌ阴ｆ矿ｏｒ咖Ｄｐ口鲥ｆｃ，而ＰＲ圩ｆ４历“口把ｄＨｏ印ｆ缸ｆ矿ｓ００幽ｏｗｕｎ浙倦ｆ劬ｓ“拍Ｄ“２ｊ５００６，ｃＡｆ”口ｊ

３．Ｄ印ｎｒｆｍｅｎｆ

Ｄ厂ｓ幻，，ｌ咖ｆ哪而舭Ｈ０印ｆ砒ｍ蛔Ｕｈｆｖｅ圩砂Ｄ，胁加∞跏咖”４４２００以∞ｆＨ日

ａ

Ａｂｓｔｒａｃｔ

ｏｎ

ＴｈｅｔｌｌｅｏｒｅｔｉｃａｌｆｏｕｎｄａｔｉｏｎｏｆｒｅｌｅｖａｎｔｐａｃｋａｇｅｓｏｆＲｓｏｆｈｖａｒｅｆｏｒｎｅｈｒｏｒｋｍｅｔａ—ａｎａｌｙＳｉｓｉｓｍａｉｌｌｌｙｂａｓｅｄ

ｆｅｗｏｆｔｈｅｍ

ｕｓｅ

Ｂａｙｅｓｉａｎｓｔａｔｉｓｔｉｃａｌｍｏｄｅｌａｎｄｇｅｎｅｒａｌｉｚｅｄｌｉｎｅａｒｍｏｄｅｌ．Ｎｅｔｗｏｒｋｍｅｔａ—ａｎａｌｙｓｉｓ

ｃａｎ

ｉｓｐｅｒｆｏｒｍｅｄｕｓｉｎｇ

ＧｅＭＴｃＲＰａｃｋａｇｅｔｌｌｒｏｕｇｈｃａｌｌｉｎｇｍｅｃｏｒｒｅｓｐｏｎｄｉｎｇ

ｒｊａｇｓ

ｐａｃｋａｇｅ，ＢＲｕｇｓＰａｃｋａｇｅ，ｏｒ

Ｒ２ｗｉｎＢｕＧｓｐａｃｋａｇｅ（ｎａｍｅｌｙ，

ｇｅｎｅｒａｔｅｄａｔａｓｔｏｒａｇｅ丘ｌｅｓ

ｏｎ

ＩＡＧｓ，ＯｐｅｎＢｕＧｓ，ａｎｄｗｉｎＢｕＧｓｓｏ觚ａｒｅ，ｒｅｓｐｅｃｔｉＶｅｌｙ）．Ｍｅａｎｗｈｉｌｅ，ＧｅＭＴｃＲｐａｃｋａｇｅ

ｕｓｅ

ｆｏｒＧｅＭＴＣｓｏｆｈⅣａｒｅ．Ｔｅｃｈｏｎｉｃａｌｌｙ，ｎｅｔｗｏｒｋｍｅｔａ—ａｎａｌｙｓｉｓｉｓｐｅｒｆｏｒｍｅｄｔｈｒｏｕｇｈｃａｌｌｉｎｇｔｈｅｓｏｆｈｖａｒｅｂａｓｅｄＣｈａｉｎＭｏｎｔｅＣａｒｌｏｍｅｔｌｌｏｄ．Ｉｎｔｈｉｓａｎｉｃｌｅ，ｗｅｂｒｉｅｎｙｉｎｔｒｏｄｕｃｅｈｏｗｔｏａｎａｌｙＳｉｓ

Ｍａｒｋｏｖ

ＧｅＭＴＣＲｐａｃｋａｇｅ

ｔｏ

ｐｅｒｆｏｒｍｎｅｔｗｏｒｋｍｅｔａ—

ｔｌｌｒｏｕｇｈｃａｌｌｉＩｌｇｔ１１ｅＯｐｅｎＢＵＧＳｓｏｆｎｖａｒｅ．

ＲｓｏＲｗａｒｅ；Ｎｅｔｗｏｒｋｍｅｔａ—ａｎａｌｙＳｉｓ；ＧｅＭＴＣｐａｃｋａｇｅ

Ｋ吖ｗｏｒｄｓ

随着网状Ｍｅｔａ分析…的发展，关于网状Ｍｅｔａ分析的软件日益增多。Ｒ软件因其开放性和灵活性使其拥有广泛的程序包，如ｎｌｍｅ程序包、Ｒ２ｗｉｎ—ＢｕＧｓ程序包、ｇｅｍｔｃ程序包等。ｇｅｍｔｃ程序包是通过相关指令来调用基于马尔科夫链．蒙特卡罗

（ＭａｒｋｏｖｃｈａｉｎＭｏｎｔｅ

ｌ安装及加载

１．１软件安装

Ｒ软件当前最新版本为Ｒ一３．ｏ．１，用户可从官方网站ｈｔｔｐ：／／ｗⅥＭｒ．ｐｒｏｊｅｃｔ．ｏｒｇ上获取最新的版本及相应的统计包【６】。ＯｐｅｎＢｕＧｓ当前最新版本为３．２．２，

ｃａｒｌｏ，ＭｃＭｃ砝的相关软件，

如ＩＡＧｓ软件、ＯｐｅｎＢｕＧｓ软件及ｗｉｎＢｕＧｓ软件来实现网状Ｍｅｔａ分析。现已有介绍ｎｌｍｅ程序包【２】和Ｒ２ｗｉｎＢｕＧＳ程序包【３】的文献，本文对Ｒ软件ｇｅｍｔｃ程序包【４，５】在网状Ｍｅｔａ分析中的应用进行简介，调用的是ＯｐｅｎＢＵＧｓ软件。

ＤＯＩ：ｌＯ．７５０７／１６７２－２５３１．２０１３０２１５

下载网址为ｈｔｔｐ：／／Ｍ州．ｏｐｅｎｂｕｇｓ．ｉｎｆｏ／ｗ．ｃｇｉ／Ｄｏｗｎ一

１０ａｄｓ。从官方网站下载软件后，安装到本地计算机。

安装完毕后，Ｒ软件需要进一步安装和加载所需的

程序包。

１．２程序包加载

软件安装完毕后，双击桌面的Ｒ软件图标，可启动Ｒ的交互式窗口（Ｒ—ＧｕＩ）。于命令提示符“＞”后输人命令ｉｎｓｔａｌｌ．Ｐａｃｋａｇｅｓ（“ｇｅｍｔｃ”），选择镜像安装（ｃＲＡＮ），安装完成后可由１ｉｂｒａｒｙ（“ｇｅｍｔｃ”）命

令完成加载。注意在加载ｇｅｍｔｃ包之前，还需安装

基金项目：湖北省教育科学“十二五”规划２０１２年度重点课题（编号２０１２Ａ０５０），湖北医药学院２０１１年度优秀中青年科技创新团队项目（编号：２０１１ｃｚｘ０１），湖北医药学院附属太和医院２０１３年循证医学“苗圃基金”项目（编号：ＥＢＭ２０１３００４）。

作者简介：张超，男（１９８８年～），本科，住院医师，以循证医学、系统评价与Ｍｅｔａ分析为主要研究方向。Ｅｍ日ｆｆ：拍ｎ，≈ａ幽ｎ００８０３＠ｊ２矗ｃｏ“４通讯作者，ＥⅢ日ｆｆ：ｚＰＨ麟ｆｎｎ缸ＤｊＪ２８＠Ｊ６３．ｃｏ卅

１ａｔｔｉｃｅ及ｃｏｄａ程序包，具体方法如上。

ｃｌｔＩ讯

１２５８

＠２０１３Ｅｄｉｆｏｒ｛ｎｆＢＤａｒｄ

ｏ，Ｃｈ｛ｎ，．Ｅｖ｛ｄ—ｂＤｓｅｄＭｅｄ

ｗｗｗ巧ｅ扫卅．Ｄ，哥ｃＨ

网状Meta分析

２数据处理

我们先前介绍了命令读取和手动输入两种数据录入方式【３】。但ｇｅｍｔｃ程序包的数据录入方式有所不同。本文同样引用‘Ｒ软件Ｒ２ｗｉｎＢｕＧｓ程序包在网状Ｍｅｔａ分析中的应用”文中【３】的实例为例，具体操作方法如下：２．１治疗方案标识

首先，我们需要对治疗方案进行相关标记，此

处我们需要用到ｒｅａｄ．ｔａｂｌｅ（）命令和ｔｅ】【ｔｃｏｎｎｅｃ—

ｔｉｏｎ（）选择项，对数据进行读取并进行关联。我们所引用的实例中的１４种药物的对应关系如下：

Ａ—ｐｌａｃｅｂｏ，Ｂ－ｂｕｐｒｏｐｉｏｎ，Ｃ—ｃｉｔａｌｏｐｒａｍ，Ｄ－ｄｅｓｖｅｎｌａｆａｘｉｎｅ，Ｅ．ｄｕｌｏｘｅｔｉｎｅ，Ｆ－ｅｓｃｉｔａｌｏｐｒａｍ，Ｇ．ｆｌｕｏｘｅｔｉｎｅ，Ｈ—ｆｌｕｖｏｘａｍｉｎｅ，Ｉ—ｍｉｒｔａｚａｐｉｎｅ，７一ｎｅｆａｚｏｄｏｎｅ，Ｋ ｐａｒｏｘｅｔｉｎｅ，Ｌ—ｓｅｒｔｒａｌｉｎｅ，Ｍ－ｔｒａｚｏｄｏｎｅ，Ｎ－

ｖｅｌａｆａ】【ｉｎｅ

其在Ｒ软件运用ｇｅｍｔｃ程序包的录入标记命令如下：

ｔｒｅａｔｍｅｎｔｓ＜一ｒｅａｄ．ｔａｂｌｅ（ｔｅ】（ｔＣｏｎｎｅｃｔｉｏｎ（

’ｉｄｄｅｓｃｒｉＰｔｉｏｎＡｐｌａｃｅｂｏＢｂｕＰｒｏＰｉｏｎＣｃｉｔａｌｏｐｒａｍＤ

ｄｅｓｖｅｌｌｌａｆ技ｉＩｌｅ

ＥｄｕｌｏｘｅｔｉｎｅＦｅｓｃｉｔａｌｏｐｒａｍＧｎｕｏｘｅｔｉｎｅＨｎｕｖｏｘａｍｉｎｅＩｍｉｒｔａｚａｐｉｎｅＴｎｅｆａｚｏｄｏｎｅＫＰａｒｏｘｅｔｉｎｅＬｓｅｒｔｒａｌｉｎｅ】、，ＩｔｒａｚｏｄｏｎｅＮｖｅｌａｆａｘｉＩｌｅ’

），ｈｅａｄｅｒ＝ＴＲＵＥ）

需要注意的是：ｔｅ】【ｔＣｏｎｎｅｃｔｉｏｎ（”）里面为英文

模式的单引号，ｉｄ标识列与ｄｅｓｃｒｉｐｔｉｏｎ标识列之间均以空格键一一隔开，如Ａ与ｐｌａｃｅｂｏ间空一格。

２．２数据录入

同样，以上述相同方式进行录人。其所需录

入的列依次为‘Ｓｔｕ酊’、‘．Ｉ＇ｒｅａｔｍｅｎｔ，’、‘仅ｅｓｐｏｎｄｅｒＳ’’、

‘Ｓａｍｐｌｅ

ｓｉｚｅ’’，分别代表研究组、治疗方案、阳性组

数及总例数。具体数据及格式见表１。

２．３

ｎｅ咖ｒｋ数据创建

重季亳丝矽：，：篡弓三嚣三医学杂志编辑部

与常用的数据存储不同，ｇｅｍｔｃ程序包其数据的存储方式较为严格，其采用的方式为ＧｅＭＴＣ软件【４】的数据识别存储模式。对于ｎｅｔｗｏｒｋ数据创建需要运用ｍｔｃ．ｎｅｔｗｏｒｋ（）命令来进行执行，命令如

下：

ｎｅｔｗｏｒｋ＜－ｍｔｃ．ｎｅｔｗｏｒｋ（ｄａｔａ，ｄｅｓｃｒｉＰｔｉｏｎ＝”Ｅｘ—ａｍｐｌｅ”，ｔｒｅａｔｍｅｎｔｓ＝ｔｒｅａｔｍｅｎｔｓ）

ｄａｔａ为录入的研究数据，ｄｅｓｃｒｉｐｔｉｏｎ为描述研究的名称，ｔｒｅａｔⅡｌｅｎｔｓ为录人的治疗方案标识。如此，一个ｎｅｔｗｏｒｋ数据已成功创建。２．４生成ｇｅｍｔｃ格式数据

Ｒ软件还具有另一项功能，就是可以为那些数据录入繁琐的软件生成相关格式数据存储文件，从而为这些软件的数据录人提供方便。ｇｅｍｔｃ程序包可以为ＧｅＭＴｃ软件生成相关的‘：ｇｅｍｔＣ’格式存储文件，具体命令如下：

ｗｒｉｔｅ．ｍｔｃ．ｎｅｔｗｏｒｋ（ｎｅｔｗｏｒｋ，”Ｃ：／Ｕｓｅｒｓ／Ａｄｍｉｎｉｓ—

ｔｒａｔｏｒ／Ｄｅｓｋｔｏｐ／Ｒｗｏｒｌ【／ｍｅ．ｇｅｍｔｃ”）

其中：ｎｅ似ｏｒｋ为事先临时存放的数据集。因默认保存路径在软件深层，难以找到，建议以‘℃：／ｕｓｅｒｓ／Ａｄｍｉｎｉｓｔｒａｔｏｒ，Ｄｅｓｋｔｏｐ／Ｒｗｏｒｌ（，ｆｉｌｅ．ｇｅｍｔｃ，’路径保存，‘矗１ｅ．ｇｅｍｔＣ’为储存数据文件。２．５ｎ曲帅ｒｋ数据读取

ｎｅ№ｒｋ数据读取与以往的读取不同，其需要

命令‘‘ｒｅａｄ．ｍｔｃ．ｎｅｔｗｏｒｋ（）’’来读取，具体命令代码如

下：

ｎｅｔｗｏｒｋ＜－ｒｅａｄ．ｍｔｃ．ｎｅｔｗｏｒｋｒＣ：／Ｕｓｅｒｓ／Ａｄｍｉｎ—

ｉｓｔｒａｔｏｒ／Ｄｅｓｋｔｏｐ／Ｒｗｏｒｋ／ｆｉｌｅ．ｇｅｍｔｃ”）“ｒｅａｄ．ｍｔｃ．ｎｅｔ－ｗｏｒｋ（）”

为使用ｎｅｔ、∞ｒｋ数据格式来进行读取。需要说明的是，本例至２．４已经完成。２．４中介绍的数据保存，若再次仍采用Ｒ读取，还可以使用‘：戗ｔ’’格式。２．５是介绍再次使用Ｒ读取采用２．４保存的文件源的方法。

３实现分析

ｇｅｍｔｃ程序包自身并没有运算程序功能，需要ｍｏｄｅｌ进行相关设置及需要相关指令调用其他网状ｍｅｔａ分析软件来进行相关的迭代运算。

３．１

设置ｎｅ咖ｒｋｍｏｄｅｌ

在选取具有网状ｍｅｔａ分析功能软件之前，我

们需要对ｍｏｄｅｌ进行相关设置，需要‘‘ｍｔｃ．ｍｏｄｅｌ（）”命令，具体命令如下：

ｍｏｄｅｌ＜一ｍｔｃ．ｍｏｄｅｌ（ｎｅｔｗｏｒｋ，ｔｙｐｅ＝”Ｃｏｎｓｉｓｔ—ｅｎｃｙ”，ｎ．ｃｈａｉｎ＝３）

Ｃ睢矗朋¨

．１２５９．¨

网状Meta分析

表１本文使用的数据及排列方式

ｓｍｄｙ

Ｔ掀吐ｍ∞ｔ

Ｒ髂ｐｏｎｄ哪ｓａｍ砖ｅｓｉ臻

ｓｔ嘶Ｔ鼢乜ｌｌｅｎｔＲｅｓＰ∞ｄｅ娼Ｓ姗ｐｌｅｓｉｚｅ

Ｓｔｕｄｙ

Ｔ嗽衄１ｅｎｔ

Ｒ唧ｏｎｄ鹏

ｓａｍｐｌｅｓＩｚｅ

１Ａ７３１５２１９Ｇ５２１０３４２Ｉ３５６６ｌＢ７６１５０１９Ｎ５４１００４３Ｇ２７６ｌ１Ｇ８３１５４２０Ａ５１８４３Ｉ２９６４２Ａ６６１２４２０Ｈ９１８４４Ｇ６７１０１２Ｂ７８１２２２１Ａ１５４２４４Ｋ６７１０２２Ｌ６６１１８２１Ｔ２５３９４５Ｇ２７４５３Ａ５５１２ｌ２２Ａ１４４５４５Ｋ３０４５３Ｂ７７１２０２２Ｉ４１９０４６Ｇ２６５０３Ｌ７９１１９２３Ａ１２５６４６Ｋ２５５０４Ａ４０１６ｌ２３Ｋ２４５５４７Ｇ５７９２４Ｄ２０５３２４２４Ａ４５１２９４７Ｋ６４９６５Ａ３９１２２２４Ｌ７０１２９４７Ｌ７０９６５Ｄ５２１２５２５Ａ１６４９４８Ｇ３５１２０６Ａ４８１２６２５Ｌ１９４９４８Ｌ４８１１８６Ｄ１４２２４９２６Ａ４９１５０４９Ｇ６３１４４７Ａ３６１２ｌ２６Ｌ７７１４９４９Ｌ７３１４２７Ｄ４６１２３２７Ａ４３１２９５０Ｇ３ｌ５４８Ａ６１１６４２７Ｌ６５１３２５０Ｎ３６５５８Ｄ１３２３１５２８Ａ１３１１６５ｌＧ３５４７８Ｅ７４１５９２８Ｌ２６ｌｌｌ５ｌＮ３５４０９Ａ５４１４ｌ２９Ａ２９１０２５２Ｇ９８１７０９Ｅ５５１４１２９Ｎ５３９５５２Ｎ８ｌ１７１１０Ａ４９１３９３０Ｂ３７６１５３Ｇ１５３１８６１０Ｅ６４１２８３０Ｇ３５６２５３Ｎ１７０１９６１ｌＡ２６１２２３ｌＢ８ｌ１２２５４Ｇ９５１６１１１Ｅ５４１２３３１Ｌ９３１２６５４Ｎ１０７１５３１２Ａ４４１３７３２Ｂ

３３６３５５Ｇ３４７３１２Ｅ１１７２７３３２Ｍ

２ｌ６１５５Ｎ４８７３１２Ｆ１１２２７４３３Ｃ８７１２０５６Ｉ７４１３９１３Ａ２４７０３３Ｆ８３１２０５６Ｋ６６１３６１３Ｅ３２７０３４Ｃ３３１０８５７Ｉ

６１１００１３Ｇ１５３３３４Ｈ３１１０９５７Ｍ

５１１００１４Ａ５ｌ９９３５Ｅ６６１３８５８Ｉ１１２０１４Ｅ１２９１９６３５Ｆ８３１４０５８Ｋ１６２０１４Ｋ５９９７３６Ｅ８ｌ１５１

５９Ｉ４２７８１５Ａ４１９３３６Ｆ９４

１４４

５９Ｌ４ｌ８２１５Ｅ１２６１８８３７Ｅ１４４

２３８６０Ｋ４８５３１５Ｋ６３８６３７Ｋ１５７２４０６０Ｍ４８５５１６Ａ１８７８３８Ｆ９４１２３６１Ｌ

３７６０１６Ｇ１３２２８５３８Ｇ８９１１７６ｌＭ

４６６２１７ＡｌＯ１９３９Ｆ１７５２３２６２Ｌ４１７２１７Ｇ３ｌ５４３９Ｋ１４６２２７６２Ｎ４９７５１７Ｋ３２５５４０Ｆ７５１０７６３Ｌ５６７９１８Ａ３７１０２４０Ｌ７４１０８６３Ｎ５６８４１８Ｇ４５１０４４ｌ

Ｆ５９９８６４Ｌ４５８２１８Ｎ５１１０２４ｌ

Ｎ４７１００６４

Ｎ

４９

７８

１９

Ａ

４ｌ

９８

４２

Ｇ

３０

６６

Ａ：ｐｌａｃｅｂｏ；Ｂ：ｂｕＰｒｏｐｉｏｎ；Ｃ：ｃｉｔａｌｏｐｒａｍ，；Ｄ：ｄｅｓＶｅｎｌａｆａｘｉｎｅ：Ｅ：ｄｕｌｏｘｅｔｉｎｅ：Ｆ：ｅｓｃｉｔａｌｏＰｒａｍ；Ｇ：ｎｕｏｘｅｔｉｎｅ；Ｈ：ｎｕＶｏｘａｍｉｎｅ；Ｉ：ｍｉｒｔａｚａｐｉｎｅ

Ｊ：ｎｅｆａｚｏｄｏｎｅ：Ｋ：Ｐａｒｏｘｅｔｉｎｅ：Ｌ：ｓｅｒｔｒａｌｉｎｅ；Ｍ：ｔｒａｚｏｄｏｎｅ；Ｎ：Ｖｅｌａｆ奴ｉｎｅ

。２叫３副打”抽旧∞州叻ｍ，篇；嬲嚣凰回

网状Meta分析

其中，ｎｅｔｗｏｒｋ为ｎｅＭｏｒｋ数据，帅ｅ为是否选

取一致性模型，ｎ．ｃｈａｉｎ为迭代运算中链的条数。需要注意的是，目前该程序包通过Ｒ软件只能使用一致性模型，非一致性模型（ｉｎｃｏｎｓｉｓｔｅｎｃｙｍｏｄｅｌ）尚不能使用。

３．２软件选取

‘‘ｒｊａ￥’、‘ＢＲｕ∥’、‘砣ｗｉｎＢｕＧＳ，’依次为调用，ＡＧｓ、

ＯｐｅｎＢｕＧＳ及ｗｉｎＢｕＧＳ软件。当‘‘Ｓａｍｐｌｅ，’设为ＮＡ时，则为依次调用上述三种软件进行迭代运算。上述操作计算机将自动完成建模、数据加载及初始化等过程。

需要注意的是：调用三种软件的前提是计算机

系统盘上必须安装相应的软件及在Ｒ软件拥有相

ｇｅｍｔｃ程序包提供了３种软件的调用方式，依次为ＪＡＧｓ、ＯｐｅｎＢｕＧｓ及ｗｉｎＢｕＧｓ软件，通过‘‘ｍｔｃ．ｍｎ（）＇’命令选取并执行。对于调用软件的选取可通过设置参数‘‘Ｓａｍｐｌｅｆ’值来实现，其取值

关的程序包。由于ｗｉｎＢｕＧｓ软件已停止更新，对于有些复杂数据可能存在无法处理的情况，相比之下，笔者建议最好能够选用ＯｐｅｎＢｕＧＳ软件。当然，

＞ｒｅｓｕｈｓ＜一ｍｔｃ．ｒｕｎ（ｍｏｄｅｌ，ｓａｍｐｌｅｒ二”ＢＲｕｇｓ”，ｎ．ａｄａｐｔ＝５０００，ｎ．ｉｔｅｒ＝２００００，ｔｈｉｎ＝１）ＷｅｌｃｏｍｅｔｏＢＲｕｇｓｃｏｎｎｅｃｔｅｄ

ｍｏｄｅｌｉｓｓｙｎｔａｃｔｉｃａｌｌｙｃｏｒｒｅｃｔｄａｔａｌｏａｄｅｄｍｏｄｅｌｃｏｍＰｉｌｅｄ

ｔｏ

０ｐｅｎＢＵＧＳＶｅｒｓｉｏｎ３．２．２

【１］”Ｄ：＼＼ｕｓｅｒｓ＼＼ＡＤＭＩＮＩ～１＼＼ＡｐｐＤａｔａ＼＼Ｌｏｃａｌ＼＼Ｔｅｍｐ＼＼ＲｔｍｐｕＣｑ９５ｗ／ｉｎｉｔｓｌ．伍ｔ”［２］”Ｄ：＼＼Ｕｓｅｒｓ＼＼ＡＤＭＩＮＩ～１＼＼ＡｐｐＤａｔａ＼＼Ｌｏｃａｌ＼＼Ｔｅｍｐ＼＼Ｒｔｍｐｕｃｑ９５ｗ／ｉｎｉｔｓ２．ｔ】【ｔ”【３］¨Ｄ：＼＼ＵＳｅｒｓ＼＼ＡＤＭＩＮＩ～ｌ＼＼ＡｐｐＤａｔａ＼＼Ｌｏｃａｌ＼＼ＴｅｍＰ、＼Ｒｔｍｐｕｃｑ９５ｗ／ｉｎｉｔｓ３．戗ｔ”

Ｉｎｉｔｉａｌｉｚｉｎｇｃｈａｉｎ１：

ｉｎｉｔｉａｌｖａｌｕｅｓ１０ａｄｅｄａｎｄｃｈａｉｎｉｎｉｔｉａｌｉｚｅｄｂｕｔａｎｏｔｈｅｒｃｈａｉｎｃｏｎｔａｉｎｕｎｉｎｉｔｉａｌｉｚｅｄｖａｒｉａｂｌｅｓＩｎｉｔｉａｌｉｚｉｎｇｃｈａｉｎ２：

ｉｎｉｔｉａｌｖａｌｕｅｓｌｏａｄｅｄａｎｄｃｈａｉｎｉｎｉｔｉａｌｉｚｅｄｂｕｔａｎｏｔｈｅｒｃｈａｉｎｃｏｎｔａｉｎｕｎｉｎｉｔｉａｈｚｅｄｖａｒｉａｂｌｅｓＩｎｉｔｉａｌｉｚｉｎｇｃｈａｉｎ３：ｍｏｄｅｌｉｓｉｎｉｔｉａｌｉｚｅｄ

ｍｏｄｅｌｉｓａｌｒｅａｄＶｉｎｉｔｉａｌｉｚｅｄ

５０００

ｕｐｄａｔｅｓｔｏｏｋ１２

ｓ

ｄｅ、ｒｉａｎｃｅｓｅｔ

ｎｌｏｎｉｔｏｒｓｅｔ

ｆｂｒｖａｒｉａｂｌｅ’ｄ．Ａ．Ｂ’

ｍｏｎｉｔｏｒｓｅｔｆｏｒｖａｒｉａｂｌｅ’ｄ．Ａ．Ｄ’ｍｏｎｉｔｏｒｓｅｔｆｏｒｖａｒｉａｂｌｅ‘ｄ．Ａ．Ｅ‘ｍｏｎｉｔｏｒｓｅｔｆｏｒｖａｒｉａｂｌｅ’ｄ．Ａ．Ｆ’ｍｏｎｉｔｏｒｓｅｔｆｂｒｖａｒｉａｂｌｅ’ｄ．Ａ．Ｇ’ｍｏｎｉｔｏｒｓｅｔｆｏｒｖａｒｉａｂｌｅ’ｄ．Ａ．Ｈ’ｍｏｎｉｔｏｒｓｅｔｆｏｒｖａｒｉａｂｌｅ’ｄＡ．Ｔ’ｍｏｎｉｔｏｒｓｅｔｆｏｒｖａｒｉａｂｌｅ’ｄ．Ａ．Ｋ’ｍｏｎｉｔｏｒｓｅｔｆ’ｏｒｖａｒｉａｂｌｅ１ｄ．Ａ．Ｌ１ｍｏｎｉｔｏｒｓｅｔｆｏｒｖａｒｉａｂｌｅ’ｄ．Ａ．Ｎ’ｍｏｎｉｔｏｒｓｅｔｆｏｒｖａｒｉａｂｌｅ’ｄ．Ｂ．Ｍ’ｍｏｎｉｔｏｒｓｅｔｆｏｒｖａｒｉａｂｌｅ’ｄ．ＦＣ’ｍｏｎｉｔｏｒｓｅｔｆｏｒｖａｒｉａｂｌｅ’ｄ．Ｇ．１１ｍｏｎｉｔｏｒｓｅｔｆ－ｏｒｖａｒｉａｂｌｅ’ｓｄ．ｄ’

２００００

ｕｐｄａｔｅｓｔｏｏｋ５２

ｓ

图１Ａ

ＯｐｅｎＢＵＧＳ行网状Ｍｅｔａ分析的迭代进程

＠２０１３中国循证医学杂志编辑部ｗｗ砒巧Ｐ扫卅．Ｄ７＿ｇｃｎ

ＣＪＥｔＭ

１２６ｌ

网状Meta分析

选择ｗｉｎＢｕＧｓ软件仍是可以的。

３．３实现分析

数，‘协ｉｄ’为步长。命令运行的过程及结果见图１。

本例以调用ＯｐｅｎＢｕＧｓ为例作为演示，具体命

令如下：

ｒｅｓｕｌｔｓ＜一ｍｔｃ．ｒｕｎ（ｍｏｄｅｌ，ｓａｍｐｌｅｒ＝”ＢＲｕｇｓ”，ｎ．ａｄａｐｔ＝５０００，ｎ．ｉｔｅｒ＝２００００，ｔ１１ｉｎ＝１）

４绘制图形

在完成上述命令操作后，我们基于Ｒ软件绘图功能，对数据结果进行相关图形的绘制。需要说明的是，以下截取的图形均可采用Ｒ软件保存图形的‘‘ｎ．ａｄａｐｔ，’为预迭代次数，‘‘ｎ．ｉｔｅ，’为迭代运算次

ＣＪＥＯＭ

１２６２

方式选择对应的位置及格式进行保存。

ＩｔｅｒａｔｉＯｎｓ＝５００１：２５０００１１１ｊｎｎｉｎｇｍｔｅｒｖａｌ＝ｌＮｕｍｂｅｒｏｆｃｈａｉｎｓ＝３

ＳａｍＰｌｅｓｉｚｅｐｅｒｃｈａｍ＝２００００

１．Ｅｍｐｍｃａｌｍｅａｎａｎｄｓｔａｎｄａｒｄｄｅｖｉａｔｉｏｎｆｏｒｅａｃｈｖａｒｉａｂｌｅ，

Ｐ１ｕｓｓｔａｎｄａｒｄｅＨｏｒｏｆｔｌｌｅｍｅａｎ：

Ｍｅａｎ

ＳＤＮ撕ｅ

ＳＥ

Ｔｉｍｅ—ｓｅｒｉｅｓＳＥｄ．Ａ．Ｂ０．５８８４０．１６０６５Ｏ．０００６５５８０．００１７７０１ｄ．Ａ．Ｄ０．６８７７０．１５１６３０．０００６１９０．００１２２ｌｌｄ．Ａ．Ｅ０．５８０２０．１１３３９０．０００４６２９０．００１１４３４ｄ．Ａ．Ｆ０．８７５６０．１４９３２Ｏ．０００６０９６０．００１５９６４ｄ．Ａ．Ｇ０．５００３０．１１０７６０．０００４５２２０．００１５１９７ｄ．Ａ．Ｈ０．９５ｌ０．４５８０８０．００１８７０ｌ０．００６６８８２ｄ．Ａ．Ｔ０．６５９４０．２２３４４０．０００９１２２Ｏ．００２６１５５ｄ．Ａ．ＫＯ．６４７６Ｏ．１３４６６０．０００５４９７Ｏ．００１６６１ｄ．Ａ．Ｌ０．６７７１Ｏ．１０８８ｌＯ．０００４４４２０．００１２５１４ｄ．Ａ．Ｎ

０．７６１７０．１３２５４０．０００５４１ｌ０．００１５９７６ｄ．Ｂ．Ｍ

．Ｏ．１２６３Ｏ．２７０２５Ｏ．００１１０３３０．００３２１１２ｄ．Ｆ．Ｃ０．１６８９Ｏ＿３５２７８Ｏ．００１４４０２０．００４３６３５ｄ．Ｇ．Ｉ０．３４２８０．２５１７７０．００１０２７９０．００２４９５３ｓｄ．ｄ

０．２５６５

Ｏ．０４９４３

０．０００２０１８

０．０００７４７９

２．Ｑｕａｎｔｉｌｅｓｆｏｒｅａｃｈｖａｒｉａｂｌｅ：

２．５０％

２５％５０％７５％９７．５０％ｄ．Ａ．Ｂ０．２７３７４Ｏ．４８０７５０．５８７６Ｏ．６９５４０．９０８２ｄ．Ａ．ＤＯ．３８５３６０．５８８６８０．６８８６Ｏ．７８７７０．９８６３ｄ．ＡＥＯ．３５７６９０．５０４６８Ｏ．５８０６Ｏ．６５６４０．８０３３ｄ．Ａ．Ｆ０．５８４３４０．７７４９２Ｏ．８７４５Ｏ．９７５２１．１７３２ｄ．Ａ．Ｇ０．２８２２０．４２６０３０．５００９０．５７４２０．７１７３ｄ．Ａ．ＨＯ．０５１５７０．６４７０３Ｏ．９５０３１．２５６ｌ１．８５４６ｄ．Ａ．Ｔ０．２２６２７０．５０８５６０．６５７２０．８０９ｌ１＿０９８ｌｄ．Ａ．Ｋ０－３８６４７０．５５７０ｌ０．６４６９０．７３５９０．９１７２ｄ．Ａ．Ｌ０．４６３２８０．６０３８７Ｏ．６７７Ｏ．７５０３Ｏ．８９１４ｄ．Ａ．Ｎ

０．５０２４１０．６７２４２０．７６１８Ｏ．８４９８１．０２１９ｄ．Ｂ．Ｍ

．Ｏ．６５９７７—０．３０７２－Ｏ．１２６２Ｏ．０５６３Ｏ．４０６ｄ．Ｆ．Ｃ．０．５２３８９ Ｏ．０６６０４０．１６８２Ｏ．４０３９０．８６０２ｄ．Ｇ．Ｉ一０．１４９９２０．１７４４６０．３４２１０．５１０８０．８４０８ｓｄ．ｄ

Ｏ．１６５１３

０．２２２９１

０．２５４２

０．２８８４

Ｏ．３５９

图１Ｂ网状Ｍｅｔａ分析结果

＠２ＤＪ３尉泐ｒｆｎｆＢＤ４耐矿劭ｆｎ，西珏６琊鲥胁ｄ

ｗｗｗＥｆ曲小．Ｄｒｇ．ｃｎ