网络Meta分析在R软件中的实现
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
Ｊｉｎｇ’ａｎ Ｄｉｓｔｒｉｃｔ Ｃｅｎｔｒａｌ Ｈｏｓｐｉｔａｌ，Ｓｈａｎｇｈａｉ ２０００４０，Ｃｈｉｎａ）
Ａｂｓｔｒａｃｔ： Ｏｂｊｅｃｔｉｖｅ Ｔｏ ｉｎｔｒｏｄｕｃｅ ｍｅｔｈｏｄｏｌｏｇｙ ｆｏｒ ｕｎｄｅｒｔａｋｉｎｇ ａ ｎｅｔｗｏｒｋ ｍｅｔａ—ａｎａｌｙｓｉｓ ａｎｄ ｉｔｓ ｒｅａｌｉｚａｔｉｏｎ ｉｎ Ｒ． Ｍｅｔｈｏｄｓ Ｂｙ ｕｓｉｎｇ ｌｍｅ（）ｆｕｎｃｔｉｏｎ ｉｎ Ｒ，ｎｅｔｗｏｒｋ ｍｅｔａ—ａｎａｌｙｓｉｓ ｏｆ ｒａｎｄｏｍｉｚｅｄ ｃｏｎｔｒｏｌｌｅｄ ｔｒｉａｌｓ（ＲＣＴｓ）ｗａｓ ｕｓｅｄ ｔｏ ｃｏｍｂｉｎｅ ｄｉｒｅｃｔ ａｎｄ ｉｎｄｉｒｅｃｔ ｅｓｔｉｍａｔｅｓ ｏｆ ｔｈｅ ｅｆｆｅｃｔ ｏｆ ｔｈｒｅｅ ｄｒｕｇｓ ａｎｄ ｐｌａｃｅｂｏ ｆｏｒ ｍｙｃｏｂａｅｔｅｒｉｕｍ ａｖｉｕｍ ｃｏｍｐｌｅｘ（ＭＡＣ）ｉｎ ｐａｔｉｅｎｔｓ ｉｎｆｅｃｔｅｄ ｗｉｔｈ ｔｈｅ ｈｕｍａｎ ｉｍｍｕｎｏｄｅ胁ｉｅｎｃｙ ｖｉｒｕｓ（ＨＩＶ）ｆｒｏｍ ５ ＲＣＴｓ．Ｒｅｓｕｌｔｓ Ａｌｌ ｄｒｕｇｓ ｄｉｆｆｅｒｅｄ ｆｒｏｍ ｐｌａｃｅｂｏ ｉｎ ｒｅｄｕｃｉｎｇ ＭＡＣ ｅｖｅｎｔｓ ｆｒｏｍ ｔｈｅ ｒｅｓｕｌｔｓ ｏｆ ｄｉｒｅｃｔ ａｎｄ ｉｎｄｉｒｅｃｔ ｃｏｍｐａｒｉｓｏｎｓ．Ｃｌａｒｉｔｈｒｏｍｙｅｉｎ ａｎｄ Ａｚｉｔｈｒｏｍｙｃｉｎ ｄｉｆｉｅｒｅｄ ｆｒｏｍ Ｒｉｆａｂｕｔｉｎ ｉｎ ｒｅｄｕｃｉｎｇ ＭＡＣ ｅｖｅｎｔｓ ｆｒｏｍ ｔｈｅ ｒｅｓｕｌｔｓ ｏｆ ｄｉｒｅｃｔ ａｎｄ ｉｎｄｉｒｅｃｔ ｃｏｍｐａｒｉｓｏｎｓ ｒｅｓｐｅｃｔｉｖｅｌｙ．Ｔｈｅ ｉｎｄｉｒｅｃｔ ｅｆｆｉｃａｃｙ ｏｆ Ｃｌａｒｉｔｈｒｏｍｙｃｉｎ ａｎｄ Ａｚｉｔｈｒｏｍｙｃｉｎ ｈａｄ ｎｏ ｓｉｇｎｉｆｉｃａｎｃｅ ｉｎ ｓｔａｔｉｓｔｉｃａｌｌｙ． Ｃｏｎｃｌｕｓｉｏｎｓ Ａ ｎｅｔｗｏｒｋ ｍｅｔａ—ａｎａｌｙｓｉｓ ｃａｎ ｂｅ ｕｓｅｄ ｔｏ ｃｏｍｂｉｎｅ ｄｉｒｅｃｔ ａｎｄ ｉｎｄｉｒｅｃｔ ｔｒｅａｔｍｅｎｔ ｅｆｆｅｃｔｓ ｉｎ ａ ｆｏｒｍａｌ ｗｈｅｎ ｓｅｖｅｒａｌ ｔｒｅａｔｍｅｎｔｓ ａｒｅ ａｖａｉｌａｂｌｅ ａｎｄ ｎｏｔ ａｌｌ ｔｒｅａｔｍｅｎｔｓ ｈａｖｅ ｂｅｅｎ ｃｏｍｐａｒｅｄ ｄｉｒｅｃｔｌｙ ｉｎ ｓｏｍｅ ｔｗｏ—ａｒｍｅｄ ＲＣＴｓ．
一 效应估计部分，以“Ａｚｉｔｈｒｏｍｙｅｉｎ”为参照药物，其他 三种药物分另４与其相对比较的ｌｏｇｏｒ及其标准误、相 应Ｐ值，从而可以求得其间接比较的比数比（ｏｄｄｓ
一一 ｒａｔｉｏ．ＯＲ）及其９５％可信区间（ｃｏｎｆｉｄｅｎｃｅ ｉｎｔｅｒｖａｌ，
呲一 ＣＩ）．如得到“Ｃｌａｒｉｔｈｒｏｍｙｃｉｎ ＶＳ．Ａｚｉｔｈｒｏｍｙｃｉｎ”间接比
比较的结果（纵向排列药物为参照）。直接比较和
间接比较均提示三种药物均较之于“Ｐｌａｃｅｂｏ”、
“Ａｚｉｔｈｒｏｍｙｃｉｎ” 和 “Ｃｌａｒｉｔｈｒｏｍｙｅｉｎ” 较之于
一一一～一一一
“Ｒｉｆａｂｕｔｉｎ”可降低发生ＭＡＣ事件的风险，而且直 接比较和间接比较结果是一致的：而 “Ａｚｉｔｈｒｏｍｙｃｉｎ”和“Ｃｌａｒｉｔｈｒｏｍｙｃｉｎ’’在降低ＭＡＣ发 生风险方而无明显差异；同时还发现间接比较的
表２预处理后的数据
３网络Ｍｅｔａ分析的具体步骤
以“ｎｌｍｅ”扩展包中ｌｍｅ（）函数来实现网络 Ｍｅｔａ分析。具体过程如下：
第二步：读人数据，将表中的数据输入Ｅｘｃｅｌ表 格中，输人完成后，将其另存为“文本文件（制表符 分隔）（术．ｔｘｔ）”，并命名为“ｈｉｖ．ｍａ．ｔｘｔ”，假设存储在 ｃ盘根目录下。在Ｒ软件命令输入窗口分别键入：
表３直接比较与间接比较获得结果比较
预措施与安慰剂比较和Ｂ干预措施与安慰剂比较 的随机对照研究，由此可能会获得Ａ干预措施和 Ｂ干预措施效应量的间接比较结果，但其证据级别 低于直接比较Ａ干预措施和Ｂ干预措施的证据级 别，因此当两种证据都存在时，首先考虑使用直接 比较的证据［５。。虽然间接比较的证据力度存在争议， 但众多学者不断探讨间接比较Ｍｅｔａ分析方法学， 目前常用的有“网络Ｍｅｔａ分析（ｎｅｔｗｏｒｋ ｍｅｔａ—
１ 数据来源
该方法需要数量较多的干预措施比较研究，
１８６
循证医学
２０１２年第１２卷第３期
但限于篇幅和简便说明起见，我们采用一简单数 据为例，其源于公开发表的文献㈦，含有五个研究， 主要观察药物预防人类免疫缺陷病毒感染者出现 复合鸟分支杆菌（ｍｙｅｏｂａｃｔｅｒｉｕｍ ａｖｉｕｍ ｃｏｍｐｌｅｘ， ＭＡＣ）事件例数，其用药物、患者数、出现ＭＡＣ事 件人数具体见表１。
２０１２年６月 第１２卷第３期
循证医学
Ｔｈｅ Ｊｏｕｒｎａｌ ｏｆ Ｅｖｉｄｅｎｃｅ—Ｂａｓｅｄ Ｍｅｄｉｅｉｎｅ
·循证医学中的医学统计学问题·
Ｊｕｎ．２０１２ Ｖ０１．１２ Ｎｏ．３
网络Ｍｅｔａ分析在Ｒ软件中的实现
张天嵩８一， 熊 茜８
（上海市静安区中心医院ａ．科教部；ｂ．中医科，上海２０００４０）
第二步：进行网络Ｍｅｔａ分析，在Ｒ软件命令 输入窗口分别键人：
＞ｌｉｂｒａｒｙ（ｎｌｍｅ） ＞ｌｍｅ １＜－ｌｍｅ（１０９０ｒ～０＋Ｒｉｆａｂｕｔｉｎ＋Ｃｌａｒｉｔｈｒｏｍｙｃｉｎ
＋Ｐｌａｃｅｂｏ， ｒａｎｄｏｍ＝～１ ｔｒｔｐａｉｒ，ｄａｔａ＝ｍｙｄａｔａ，
ｗｅｉｇｈｔｓ＝ｖａｒＣｏｎｓｔＰｏｗｅｒ（ｆｏｒｍ＝－ｓｅｌｏｇｏｒ，ｆｉｘｅｄ＝ｌｉｓｔ （ｐｏｗｅｒ＝１）））
为了进行网络Ｍｅｔａ分析．必须对表中的数据 进行预处理，使之适宜于线性混合模型分析的要 求。将数据整理成两两比较形式，共得到５个比较 组，按经典方法，分别计算出两比较组的对数比数 比（１０９０ｒ）及其标准误（ｓｅｌｏｇｏｒ）；将配对比较的药物 分别记为１与一１，无关的药物记为０。预处理后的 数据见表２。
［中图分类号］Ｒ１９５．１
［文献标识码］Ａ
［文章编号］ １６７１—５１４４（２０１２）０３—０１８５—０４
Ｎｅｔｗｏｒｋ Ｍｅｔａ－Ａｎａｌｙｓｉｓ ａｎｄ Ｉｔｓ Ｓｏｌｕｔｉｏｎ ｂｙ Ｒ
ＺＨＡＮＧ Ｔｉａｎ—ｓｏｎ９８一．ＸＩＯＮＧ Ｑｉａｎ８ （ａ．Ｄｅｐａｒｔｍｅｎｔ ｏｆＳｃｉｅｎｃｅ＆Ｅｄｕｃａｔｉｏｎ；ｂ．Ｄｅｐａｒｔｍｅｎｔ ｏｆ ＴＣＭ，
［２］Ｌｕｍｌｅｙ Ｔ－
Ｎｅｔｗｏｒｋ ｍｅｔａ—ａｎａｌｙｓｉｓ ｆｏｒ ｉｎｄｉｒｅｍⅡｅａ“ｎ”‘
源软件，是一套完整的数据处理、计算和绘图软件 …ｃｃｏ。ｍ。ｐＪａｒ，ｉ８蔷篡．Ｃ言？‘Ｉ：；。ｅｄ。，ｄｂ２。０昭０２Ｋ，Ａ２１，（１：？：ｉ３１篡惫。。。，。。
系统．大多数经典的统计方法和最新的技术都可ｄｉｒｅｃｔ ａｎｄ ａｄｊｕｓｔｅｄ ｉｎｄｉｒｅｃｔ ｃ。ｍｐａｒｉｓｏｎ ｔ。ｐｌａｃｅｂｏ： Ａｎ
上述命令是载入“ｎｌｍｅ”扩展包，采用ｌｍｅ（）函 数拟合线性混合效应模型。其中“ｌｏｇｏｒ～０＋ Ｒｉｆａｂｕｔｉｎ＋Ｃｌａｒｉｔｈｒｏｍｙｅｉｎ＋Ｐｌａｃｅｂｏ”为固定效应：
张天嵩，等．网络Ｍｅｔａ分析在Ｒ软件中的实现
１８７
在结果的随机效应估计部分．其截距标准差代 表网络不相关性，其数值越小表明网络越稳定；固定
网络Ｍｅｔａ分析是由Ｌｕｍｌｅｙ建立的一种用于 间接比较的方法。可以获得直接和间接比较的结 果。该方法基于线性混合效应模型．包括两个随机
１８８
循汪医学
２０１２年第１２卷第３期
效应模型：研究水平的异质性、直接比较与间接比
［参考文献］
较结果比较的非一致性（即不相关性），只有当网 络中的比较对象形成一个闭合环（如图１）Ｈｅ，才
Ｋｅｙ ｗｏｒｄｓ：ｎｅｔｗｏｒｋ ｍｅｔａ—ａｎａｌｙｓｉｓ；ｌｉｎｅａｒ ｍｉｘｅｄ—ｅｆｆｅｃｔｓ ｍｏｄｅｌ；ｉｎｄｉｒｅｃｔ ｃｏｍｐａｒｉｓｏｎｓ；Ｒ
在理想情况下，直接比较某一干预措施与对 照或阳性药物的随机对照试验是评价治疗效应的 金标准，但在实际情况下，许多干预措施缺乏或没 有直接比较的证据…。随着医学研究的进展，在没 有直接证据时．越来越需要干预措施间接比较的 证据。采用网络Ｍｅｔａ分析可以同时进行直接与间 接比较．即使是相比较的两种治疗药物从未进行
人类免疫缺陷病毒感染者出现复合鸟分支杆菌事件发生：而对无直接比较证据的克拉霉素与阿奇霉素，间接比
较结果显示两者差异无统计学意义。结论合并多个干预措施效果和没有所有干预措施直接比较证据时，采用
网络Ｍｅｔａ分析方法合并直接和间接比较结果是一种正规的方法。
［关键词］ 网络Ｍｅｔａ分析；线性混合效应模型；间接比较；Ｒ软件
Ａｚｉｔｈｒｏｍｙｃｉｎ”直接比较的证据．但可以通过网络 Ｍｅｔａ分析来计算出间接比较的证据。
表１实例数据
图１药物直接比较来自百度文库间接比较网络图
２数据预处理
从表１中可以看出，有“Ｒｉｆａｂｕｔｉｎ与Ｐｌａｃｅｂｏ”、 “Ｃｌａｒｉｔｈｒｏｍｙｃｉｎ与Ｐｌａｃｅｂｏ”、“Ａｚｉｔｈｒｏｍｙｅｉｎ与 Ｐｌａｃｅｂｏ”、“Ｒｉｆａｂｕｔｉｎ与Ｃｌａｒｉｔｈｒｏｍｙｃｉｎ”直接比较的 证据．图１中以实线表示，但无“Ｃｌａｒｉｔｈｒｏｍｙｅｉｎ与
［１］Ｔｈｉｊｓ Ｖ，Ｌｅｍｍｅｎｓ Ｒ，ＦｉｅｕＷｓ Ｓ·Ｎｅｔｗｏｒｋ ｍｅｔａ－ａｎａｌｙｓｉｓ：
慕１：：＝：。ｒｅｉ。ｃ。ｔ‰ａ－ａｎ血ａｌｙｓ叭ｉｓ。兰：；＝：。８焉节矾繁淼。；ｉ
能估算不相关性１２］，但未提供不相关性大小的切
２００８．２９（９）：１０８６—１０９２．
割点。Ｒ软件属于ＧＮＵ系统的一个自由、免费、开
万方数据
ａｎａｌｙｓｅｓ，ＮＭＡ）”、“混合治疗比较（ｍｉｘｅｄ ｔｒｅａｔｍｅｎｔ ｃｏｍｐａｒｉｓｏｎ，ＭＴＣ）”、“多种干预措施Ｍｅｔａ分析 （ｍｕｌｔｉｐｌｅ ｔｒｅａｔｍｅｎｔｓ ｍｅｔａ—ａｎａｌｙｓｉｓ，ＭＴＣ）”等‘５。，这 些Ｍｅｔａ分析方法可以通过概率统计或贝叶斯方法 来实现［２，６Ｉ。
［摘要］ 目的介绍网络Ｍｅｔａ分析方法及其在Ｒ软件的实现。方法 以实例说明，采用Ｒ软件ｌｍｅ（）函
数，通过对来自５个预防人类免疫缺陷病毒感染者发生复合鸟分支杆菌感染事件随机对照研究的３种药物和安
慰剂进行直接和间接比较，进行网络Ｍｅｔａ分析。结果直接和间接比较结果表明，所有药物较之于安慰剂均可
以减少人类免疫缺陷病毒感染者出现复合鸟分支杆菌事件发生，克拉霉素和阿奇霉素较之于利福平均可以减少
＞ｌｉｂｒａｒｙ（ｆｏｒｅｉｇｎ） ＞ｍｙｄａｔａ＜．ｒｅａｄ．ｔａｂｌｅ（”ｅ：／ｈｉｖ—ｍａ．ｔｘｔ”，ｈｅａｄｅｒ＝ ＴＲＵＥ） ＞ｍｙｄａｔａ＜－ｔｒａｎｓｆｏｒｍ（ｍｙｄａｔａ，ｔｒｔｐａｉｒ ２ ｐａｓｔｅ （，ｄｒｕｇｌ，ｄｎｌ９２））
万方数据
上述命令是将数据读人，建立一个名为 “ｍｙｄａｔａ”的数据集，并经过处理，将两两比较的药 物结成对子，以便于进一步分析。
［作者简介］
张天嵩（１９７０一），男，山东昌邑人，医学博士， 主任医师，副教授，主要从事呼吸系统疾病 的中西医结合治疗及循证医学方法学研究。
万方数据
过直接对比．这种分析方法可以将一系列不同治 疗方法随机临床试验数据汇总，然后就给定的治 疗终点进行点估计及可信区间估计，同时对非相 关性进行评估，Ｌｕｍｌｅｙ教授首先提出了网络Ｍｅｔａ 分析（ｎｅｔｗｏｒｋ ｍｅｔａ．ａｎａｌｙｓｉｓ）的模型和方法１ ２『，可以 在多种统计软件中通过拟合线性混合效应模型而 实现，本文以实例说明以Ｒ１４．１的“ｎｌｍｅ”软件包来 进行网络Ｍｅｔａ分析的具体过程。
较的ＯＲ及９５％ＣＩ为１．０３（０．７６，１．３９）。
面一． ～一
两种药物预防ＭＡＣ事件出现直接和间接比 较结果（ＯＲ及９５％Ｃ１）见表３。表中右上部分是直 接比较的结果（横向排列药物为参照）， “Ｃｌａｒｉｔｈｒｏｍｙｃｉｎ”与“Ａｚｉｔｈｒｏｍｙｃｉｎ”没有直接比较的 结果：左下部分以浅灰色为底标记的数据是间接
结果更为精确些。
４讨论
～酡一 渤
～也一 一～一～一‰～一 一～一一５
在没有直接比较的证据或者在比较两个以上 干预措施有效性时，必然会涉及间接比较．“，间接 比较结果在理论上可能比直接比较结果准确，但 其证据级别低于直接比较的证据ｃ引．如比较Ａ和Ｂ 两种干预措施的有效性，尽管没有Ａ干预措施和 Ｂ干预措施直接比较的随机对照研究，但有Ａ干
Ａｂｓｔｒａｃｔ： Ｏｂｊｅｃｔｉｖｅ Ｔｏ ｉｎｔｒｏｄｕｃｅ ｍｅｔｈｏｄｏｌｏｇｙ ｆｏｒ ｕｎｄｅｒｔａｋｉｎｇ ａ ｎｅｔｗｏｒｋ ｍｅｔａ—ａｎａｌｙｓｉｓ ａｎｄ ｉｔｓ ｒｅａｌｉｚａｔｉｏｎ ｉｎ Ｒ． Ｍｅｔｈｏｄｓ Ｂｙ ｕｓｉｎｇ ｌｍｅ（）ｆｕｎｃｔｉｏｎ ｉｎ Ｒ，ｎｅｔｗｏｒｋ ｍｅｔａ—ａｎａｌｙｓｉｓ ｏｆ ｒａｎｄｏｍｉｚｅｄ ｃｏｎｔｒｏｌｌｅｄ ｔｒｉａｌｓ（ＲＣＴｓ）ｗａｓ ｕｓｅｄ ｔｏ ｃｏｍｂｉｎｅ ｄｉｒｅｃｔ ａｎｄ ｉｎｄｉｒｅｃｔ ｅｓｔｉｍａｔｅｓ ｏｆ ｔｈｅ ｅｆｆｅｃｔ ｏｆ ｔｈｒｅｅ ｄｒｕｇｓ ａｎｄ ｐｌａｃｅｂｏ ｆｏｒ ｍｙｃｏｂａｅｔｅｒｉｕｍ ａｖｉｕｍ ｃｏｍｐｌｅｘ（ＭＡＣ）ｉｎ ｐａｔｉｅｎｔｓ ｉｎｆｅｃｔｅｄ ｗｉｔｈ ｔｈｅ ｈｕｍａｎ ｉｍｍｕｎｏｄｅ胁ｉｅｎｃｙ ｖｉｒｕｓ（ＨＩＶ）ｆｒｏｍ ５ ＲＣＴｓ．Ｒｅｓｕｌｔｓ Ａｌｌ ｄｒｕｇｓ ｄｉｆｆｅｒｅｄ ｆｒｏｍ ｐｌａｃｅｂｏ ｉｎ ｒｅｄｕｃｉｎｇ ＭＡＣ ｅｖｅｎｔｓ ｆｒｏｍ ｔｈｅ ｒｅｓｕｌｔｓ ｏｆ ｄｉｒｅｃｔ ａｎｄ ｉｎｄｉｒｅｃｔ ｃｏｍｐａｒｉｓｏｎｓ．Ｃｌａｒｉｔｈｒｏｍｙｅｉｎ ａｎｄ Ａｚｉｔｈｒｏｍｙｃｉｎ ｄｉｆｉｅｒｅｄ ｆｒｏｍ Ｒｉｆａｂｕｔｉｎ ｉｎ ｒｅｄｕｃｉｎｇ ＭＡＣ ｅｖｅｎｔｓ ｆｒｏｍ ｔｈｅ ｒｅｓｕｌｔｓ ｏｆ ｄｉｒｅｃｔ ａｎｄ ｉｎｄｉｒｅｃｔ ｃｏｍｐａｒｉｓｏｎｓ ｒｅｓｐｅｃｔｉｖｅｌｙ．Ｔｈｅ ｉｎｄｉｒｅｃｔ ｅｆｆｉｃａｃｙ ｏｆ Ｃｌａｒｉｔｈｒｏｍｙｃｉｎ ａｎｄ Ａｚｉｔｈｒｏｍｙｃｉｎ ｈａｄ ｎｏ ｓｉｇｎｉｆｉｃａｎｃｅ ｉｎ ｓｔａｔｉｓｔｉｃａｌｌｙ． Ｃｏｎｃｌｕｓｉｏｎｓ Ａ ｎｅｔｗｏｒｋ ｍｅｔａ—ａｎａｌｙｓｉｓ ｃａｎ ｂｅ ｕｓｅｄ ｔｏ ｃｏｍｂｉｎｅ ｄｉｒｅｃｔ ａｎｄ ｉｎｄｉｒｅｃｔ ｔｒｅａｔｍｅｎｔ ｅｆｆｅｃｔｓ ｉｎ ａ ｆｏｒｍａｌ ｗｈｅｎ ｓｅｖｅｒａｌ ｔｒｅａｔｍｅｎｔｓ ａｒｅ ａｖａｉｌａｂｌｅ ａｎｄ ｎｏｔ ａｌｌ ｔｒｅａｔｍｅｎｔｓ ｈａｖｅ ｂｅｅｎ ｃｏｍｐａｒｅｄ ｄｉｒｅｃｔｌｙ ｉｎ ｓｏｍｅ ｔｗｏ—ａｒｍｅｄ ＲＣＴｓ．
一 效应估计部分，以“Ａｚｉｔｈｒｏｍｙｅｉｎ”为参照药物，其他 三种药物分另４与其相对比较的ｌｏｇｏｒ及其标准误、相 应Ｐ值，从而可以求得其间接比较的比数比（ｏｄｄｓ
一一 ｒａｔｉｏ．ＯＲ）及其９５％可信区间（ｃｏｎｆｉｄｅｎｃｅ ｉｎｔｅｒｖａｌ，
呲一 ＣＩ）．如得到“Ｃｌａｒｉｔｈｒｏｍｙｃｉｎ ＶＳ．Ａｚｉｔｈｒｏｍｙｃｉｎ”间接比
比较的结果（纵向排列药物为参照）。直接比较和
间接比较均提示三种药物均较之于“Ｐｌａｃｅｂｏ”、
“Ａｚｉｔｈｒｏｍｙｃｉｎ” 和 “Ｃｌａｒｉｔｈｒｏｍｙｅｉｎ” 较之于
一一一～一一一
“Ｒｉｆａｂｕｔｉｎ”可降低发生ＭＡＣ事件的风险，而且直 接比较和间接比较结果是一致的：而 “Ａｚｉｔｈｒｏｍｙｃｉｎ”和“Ｃｌａｒｉｔｈｒｏｍｙｃｉｎ’’在降低ＭＡＣ发 生风险方而无明显差异；同时还发现间接比较的
表２预处理后的数据
３网络Ｍｅｔａ分析的具体步骤
以“ｎｌｍｅ”扩展包中ｌｍｅ（）函数来实现网络 Ｍｅｔａ分析。具体过程如下：
第二步：读人数据，将表中的数据输入Ｅｘｃｅｌ表 格中，输人完成后，将其另存为“文本文件（制表符 分隔）（术．ｔｘｔ）”，并命名为“ｈｉｖ．ｍａ．ｔｘｔ”，假设存储在 ｃ盘根目录下。在Ｒ软件命令输入窗口分别键入：
表３直接比较与间接比较获得结果比较
预措施与安慰剂比较和Ｂ干预措施与安慰剂比较 的随机对照研究，由此可能会获得Ａ干预措施和 Ｂ干预措施效应量的间接比较结果，但其证据级别 低于直接比较Ａ干预措施和Ｂ干预措施的证据级 别，因此当两种证据都存在时，首先考虑使用直接 比较的证据［５。。虽然间接比较的证据力度存在争议， 但众多学者不断探讨间接比较Ｍｅｔａ分析方法学， 目前常用的有“网络Ｍｅｔａ分析（ｎｅｔｗｏｒｋ ｍｅｔａ—
１ 数据来源
该方法需要数量较多的干预措施比较研究，
１８６
循证医学
２０１２年第１２卷第３期
但限于篇幅和简便说明起见，我们采用一简单数 据为例，其源于公开发表的文献㈦，含有五个研究， 主要观察药物预防人类免疫缺陷病毒感染者出现 复合鸟分支杆菌（ｍｙｅｏｂａｃｔｅｒｉｕｍ ａｖｉｕｍ ｃｏｍｐｌｅｘ， ＭＡＣ）事件例数，其用药物、患者数、出现ＭＡＣ事 件人数具体见表１。
２０１２年６月 第１２卷第３期
循证医学
Ｔｈｅ Ｊｏｕｒｎａｌ ｏｆ Ｅｖｉｄｅｎｃｅ—Ｂａｓｅｄ Ｍｅｄｉｅｉｎｅ
·循证医学中的医学统计学问题·
Ｊｕｎ．２０１２ Ｖ０１．１２ Ｎｏ．３
网络Ｍｅｔａ分析在Ｒ软件中的实现
张天嵩８一， 熊 茜８
（上海市静安区中心医院ａ．科教部；ｂ．中医科，上海２０００４０）
第二步：进行网络Ｍｅｔａ分析，在Ｒ软件命令 输入窗口分别键人：
＞ｌｉｂｒａｒｙ（ｎｌｍｅ） ＞ｌｍｅ １＜－ｌｍｅ（１０９０ｒ～０＋Ｒｉｆａｂｕｔｉｎ＋Ｃｌａｒｉｔｈｒｏｍｙｃｉｎ
＋Ｐｌａｃｅｂｏ， ｒａｎｄｏｍ＝～１ ｔｒｔｐａｉｒ，ｄａｔａ＝ｍｙｄａｔａ，
ｗｅｉｇｈｔｓ＝ｖａｒＣｏｎｓｔＰｏｗｅｒ（ｆｏｒｍ＝－ｓｅｌｏｇｏｒ，ｆｉｘｅｄ＝ｌｉｓｔ （ｐｏｗｅｒ＝１）））
为了进行网络Ｍｅｔａ分析．必须对表中的数据 进行预处理，使之适宜于线性混合模型分析的要 求。将数据整理成两两比较形式，共得到５个比较 组，按经典方法，分别计算出两比较组的对数比数 比（１０９０ｒ）及其标准误（ｓｅｌｏｇｏｒ）；将配对比较的药物 分别记为１与一１，无关的药物记为０。预处理后的 数据见表２。
［中图分类号］Ｒ１９５．１
［文献标识码］Ａ
［文章编号］ １６７１—５１４４（２０１２）０３—０１８５—０４
Ｎｅｔｗｏｒｋ Ｍｅｔａ－Ａｎａｌｙｓｉｓ ａｎｄ Ｉｔｓ Ｓｏｌｕｔｉｏｎ ｂｙ Ｒ
ＺＨＡＮＧ Ｔｉａｎ—ｓｏｎ９８一．ＸＩＯＮＧ Ｑｉａｎ８ （ａ．Ｄｅｐａｒｔｍｅｎｔ ｏｆＳｃｉｅｎｃｅ＆Ｅｄｕｃａｔｉｏｎ；ｂ．Ｄｅｐａｒｔｍｅｎｔ ｏｆ ＴＣＭ，
［２］Ｌｕｍｌｅｙ Ｔ－
Ｎｅｔｗｏｒｋ ｍｅｔａ—ａｎａｌｙｓｉｓ ｆｏｒ ｉｎｄｉｒｅｍⅡｅａ“ｎ”‘
源软件，是一套完整的数据处理、计算和绘图软件 …ｃｃｏ。ｍ。ｐＪａｒ，ｉ８蔷篡．Ｃ言？‘Ｉ：；。ｅｄ。，ｄｂ２。０昭０２Ｋ，Ａ２１，（１：？：ｉ３１篡惫。。。，。。
系统．大多数经典的统计方法和最新的技术都可ｄｉｒｅｃｔ ａｎｄ ａｄｊｕｓｔｅｄ ｉｎｄｉｒｅｃｔ ｃ。ｍｐａｒｉｓｏｎ ｔ。ｐｌａｃｅｂｏ： Ａｎ
上述命令是载入“ｎｌｍｅ”扩展包，采用ｌｍｅ（）函 数拟合线性混合效应模型。其中“ｌｏｇｏｒ～０＋ Ｒｉｆａｂｕｔｉｎ＋Ｃｌａｒｉｔｈｒｏｍｙｅｉｎ＋Ｐｌａｃｅｂｏ”为固定效应：
张天嵩，等．网络Ｍｅｔａ分析在Ｒ软件中的实现
１８７
在结果的随机效应估计部分．其截距标准差代 表网络不相关性，其数值越小表明网络越稳定；固定
网络Ｍｅｔａ分析是由Ｌｕｍｌｅｙ建立的一种用于 间接比较的方法。可以获得直接和间接比较的结 果。该方法基于线性混合效应模型．包括两个随机
１８８
循汪医学
２０１２年第１２卷第３期
效应模型：研究水平的异质性、直接比较与间接比
［参考文献］
较结果比较的非一致性（即不相关性），只有当网 络中的比较对象形成一个闭合环（如图１）Ｈｅ，才
Ｋｅｙ ｗｏｒｄｓ：ｎｅｔｗｏｒｋ ｍｅｔａ—ａｎａｌｙｓｉｓ；ｌｉｎｅａｒ ｍｉｘｅｄ—ｅｆｆｅｃｔｓ ｍｏｄｅｌ；ｉｎｄｉｒｅｃｔ ｃｏｍｐａｒｉｓｏｎｓ；Ｒ
在理想情况下，直接比较某一干预措施与对 照或阳性药物的随机对照试验是评价治疗效应的 金标准，但在实际情况下，许多干预措施缺乏或没 有直接比较的证据…。随着医学研究的进展，在没 有直接证据时．越来越需要干预措施间接比较的 证据。采用网络Ｍｅｔａ分析可以同时进行直接与间 接比较．即使是相比较的两种治疗药物从未进行
人类免疫缺陷病毒感染者出现复合鸟分支杆菌事件发生：而对无直接比较证据的克拉霉素与阿奇霉素，间接比
较结果显示两者差异无统计学意义。结论合并多个干预措施效果和没有所有干预措施直接比较证据时，采用
网络Ｍｅｔａ分析方法合并直接和间接比较结果是一种正规的方法。
［关键词］ 网络Ｍｅｔａ分析；线性混合效应模型；间接比较；Ｒ软件
Ａｚｉｔｈｒｏｍｙｃｉｎ”直接比较的证据．但可以通过网络 Ｍｅｔａ分析来计算出间接比较的证据。
表１实例数据
图１药物直接比较来自百度文库间接比较网络图
２数据预处理
从表１中可以看出，有“Ｒｉｆａｂｕｔｉｎ与Ｐｌａｃｅｂｏ”、 “Ｃｌａｒｉｔｈｒｏｍｙｃｉｎ与Ｐｌａｃｅｂｏ”、“Ａｚｉｔｈｒｏｍｙｅｉｎ与 Ｐｌａｃｅｂｏ”、“Ｒｉｆａｂｕｔｉｎ与Ｃｌａｒｉｔｈｒｏｍｙｃｉｎ”直接比较的 证据．图１中以实线表示，但无“Ｃｌａｒｉｔｈｒｏｍｙｅｉｎ与
［１］Ｔｈｉｊｓ Ｖ，Ｌｅｍｍｅｎｓ Ｒ，ＦｉｅｕＷｓ Ｓ·Ｎｅｔｗｏｒｋ ｍｅｔａ－ａｎａｌｙｓｉｓ：
慕１：：＝：。ｒｅｉ。ｃ。ｔ‰ａ－ａｎ血ａｌｙｓ叭ｉｓ。兰：；＝：。８焉节矾繁淼。；ｉ
能估算不相关性１２］，但未提供不相关性大小的切
２００８．２９（９）：１０８６—１０９２．
割点。Ｒ软件属于ＧＮＵ系统的一个自由、免费、开
万方数据
ａｎａｌｙｓｅｓ，ＮＭＡ）”、“混合治疗比较（ｍｉｘｅｄ ｔｒｅａｔｍｅｎｔ ｃｏｍｐａｒｉｓｏｎ，ＭＴＣ）”、“多种干预措施Ｍｅｔａ分析 （ｍｕｌｔｉｐｌｅ ｔｒｅａｔｍｅｎｔｓ ｍｅｔａ—ａｎａｌｙｓｉｓ，ＭＴＣ）”等‘５。，这 些Ｍｅｔａ分析方法可以通过概率统计或贝叶斯方法 来实现［２，６Ｉ。
［摘要］ 目的介绍网络Ｍｅｔａ分析方法及其在Ｒ软件的实现。方法 以实例说明，采用Ｒ软件ｌｍｅ（）函
数，通过对来自５个预防人类免疫缺陷病毒感染者发生复合鸟分支杆菌感染事件随机对照研究的３种药物和安
慰剂进行直接和间接比较，进行网络Ｍｅｔａ分析。结果直接和间接比较结果表明，所有药物较之于安慰剂均可
以减少人类免疫缺陷病毒感染者出现复合鸟分支杆菌事件发生，克拉霉素和阿奇霉素较之于利福平均可以减少
＞ｌｉｂｒａｒｙ（ｆｏｒｅｉｇｎ） ＞ｍｙｄａｔａ＜．ｒｅａｄ．ｔａｂｌｅ（”ｅ：／ｈｉｖ—ｍａ．ｔｘｔ”，ｈｅａｄｅｒ＝ ＴＲＵＥ） ＞ｍｙｄａｔａ＜－ｔｒａｎｓｆｏｒｍ（ｍｙｄａｔａ，ｔｒｔｐａｉｒ ２ ｐａｓｔｅ （，ｄｒｕｇｌ，ｄｎｌ９２））
万方数据
上述命令是将数据读人，建立一个名为 “ｍｙｄａｔａ”的数据集，并经过处理，将两两比较的药 物结成对子，以便于进一步分析。
［作者简介］
张天嵩（１９７０一），男，山东昌邑人，医学博士， 主任医师，副教授，主要从事呼吸系统疾病 的中西医结合治疗及循证医学方法学研究。
万方数据
过直接对比．这种分析方法可以将一系列不同治 疗方法随机临床试验数据汇总，然后就给定的治 疗终点进行点估计及可信区间估计，同时对非相 关性进行评估，Ｌｕｍｌｅｙ教授首先提出了网络Ｍｅｔａ 分析（ｎｅｔｗｏｒｋ ｍｅｔａ．ａｎａｌｙｓｉｓ）的模型和方法１ ２『，可以 在多种统计软件中通过拟合线性混合效应模型而 实现，本文以实例说明以Ｒ１４．１的“ｎｌｍｅ”软件包来 进行网络Ｍｅｔａ分析的具体过程。
较的ＯＲ及９５％ＣＩ为１．０３（０．７６，１．３９）。
面一． ～一
两种药物预防ＭＡＣ事件出现直接和间接比 较结果（ＯＲ及９５％Ｃ１）见表３。表中右上部分是直 接比较的结果（横向排列药物为参照）， “Ｃｌａｒｉｔｈｒｏｍｙｃｉｎ”与“Ａｚｉｔｈｒｏｍｙｃｉｎ”没有直接比较的 结果：左下部分以浅灰色为底标记的数据是间接
结果更为精确些。
４讨论
～酡一 渤
～也一 一～一～一‰～一 一～一一５
在没有直接比较的证据或者在比较两个以上 干预措施有效性时，必然会涉及间接比较．“，间接 比较结果在理论上可能比直接比较结果准确，但 其证据级别低于直接比较的证据ｃ引．如比较Ａ和Ｂ 两种干预措施的有效性，尽管没有Ａ干预措施和 Ｂ干预措施直接比较的随机对照研究，但有Ａ干