剂量反应关系Meta分析在R软件中的实现

现代预防医学２０１５年第４２卷第５期ＭｏｄｅｒｎＰｒｅｖｅｎｔｉｖｅＭｅｄｉｃｉｎｅ，２０１５，Ｖｏｌ．４２，ＮＯ．５

·７８９·

·流行病与统计方法·

剂量反应关系Ｍｅｔａ分析在Ｒ软件中的实现

周权１，李卉２，罗美玲３，薛静４，黎敏５

１．常德市第一人民医院科教科，湖南常德４１５００３；２．常德市第一人民医院肿瘤科，湖南常德４１５００３；

３．长沙市疾病预防控制中心计划免疫科，湖南长沙４１０００１；４．中南大学公共卫生学院流行病与卫生统计系，湖南长沙４１００７８；

５．南昌大学第二附属医院神经内科，江西南昌３３０００６

摘要：目的介绍如何利用Ｒ软件中实现剂量反应关系的Ｍｅｔａ分析。方法拟通过实例数据的演示，在Ｒ软件中，先利用Ｍｅｔａｆｏｒ程序包进行探索性二分类Ｍｅｔａ分析，在得出意义结果的基础上再利用ｄｏｓｒｅｓｍｅｔａ程序包进行剂量反应的Ｍｅｔａ分析，并绘制非线性（和）或线性剂量反应曲线结果

Ｒ软件分析结果显示：饮酒是结直肠癌的危险因素，

０ｇ／天），每剂量反应关系的非线性检验Ｐ值处于临界水平（０．０４６）。非线性剂量反应模型显示：相对于最低饮酒量（

０．９３～１．１）、１．３（１．１３～１．４）、１．５（１．２６～１．９）；线性剂量日饮酒２０ｇ、４０ｇ、６０ｇ发生结直肠癌的相对危险度为：１（

反应模型显示：在最低饮酒量（０ｇ／天）的基础上，每增加１２ｇ，发生结直肠癌的危险增加８％，非线性或线性剂量反应曲线图形基本相似。结论Ｒ软件可以实现剂量反应关系的Ｍｅｔａ分析，功能强大，实用性强。

关键词：剂量反应关系；Ｍｅｔａ分析；系统评价；Ｒ软件中图分类号：Ｒ１９５

文献标志码：Ａ

文章编号：１００３－８５０７（２０１５）０５－０７８９－０６

Implementationofdose-responseformeta-analysisinRsoftware

ＺＨＯＵＱｕａｎ，ＬＩＨｕｉ，ＬＵＯＭｅｉ－ｌｉｎｇ，ＸＵＥＪｉｎｇ，ＬＩＭｉｎ

Departmentofscienceandeducation,TheFirstPeople'sHospitalofChangde,Changde,Hunan415003,China.

Abstract:ObjectiveＴｏｉｍｐｌｅｍｅｎｔｔｈｅｄｏｓｅ－ｒｅｓｐｏｎｓｅｆｏｒｍｅｔａ－ａｎａｌｙｓｉｓｉｎＲｓｏｆｔｗａｒｅ．MethodsＴｈｅｃｏｄｅｗａｓｉｌｌｕｓｔｒａｔｅｄｗｉｔｈｅｘａｍｐｌｅｓｏｆａｌｃｏｈｏｌｉｎｔａｋｅａｎｄｃｏｌｏｒｅｃｔａｌｃａｎｃｅｒ．Ｆｉｒｓｔｌｙ，ｔｈｅｍｅｔａｆｏｒｐａｃｋａｇｅｗａｓｕｓｅｄｔｏｅｘｐｌｏｒｅｂｉｎａｒｙｃｌａｓｓｉｆｉｃａｔｉｏｎｍｅｔａａｎａｌｙｓｉｓ．ＵｎｄｅｒｔｈｅｃｏｎｄｉｔｉｏｎｏｆａｓｉｇｎｉｆｉｃａｎｔRRｖａｌｕｅ，ｗｅｃｏｎｄｕｃｔｅｄｔｈｅｄｏｓｅ－ｒｅｓｐｏｎｓｅｍｅｔａ－ａｎａｌｙｓｉｓｗｉｔｈｄｏｓｒｅｓｍｅｔａｐｒｏｇｒａｍａｎｄｄｒｅｗｔｈｅｌｉｎｅａｒａｎｄｎｏｎｌｉｎｅａｒｄｏｓｅ－ｒｅｓｐｏｎｓｅｃｕｒｖｅｓ．ResultsＴｈｅｒｅｓｕｌｔｓａｎａｌｙｚｅｄｂｙＲｓｈｏｗｔｈａｔａｌｃｏｈｏｌｉｎｔａｋｅｉｓａｒｉｓｋｆａｃｔｏｒｆｏｒｃｏｌｏｒｅｃｔａｌｃａｎｃｅｒ，ａｎｄｔｈｅPｖａｌｕｅｆｏｒｎｏｎ－ｌｉｎｅａｒｉｔｙｔｅｓｔｏｆｄｏｓｅ－ｒｅｓｐｏｎｓｅｒｅｌａｔｉｏｎｓｈｉｐｉｓ０．０４６．Ｔｈｅｎｏｎ－ｌｉｎｅａｒｉｔｙｍｏｄｅｌｓｈｏｗｓｔｈａｔｃｏｍｐａｒｅｄｗｉｔｈｎｏｎ－ｄｒｉｎｋｅｒｓ，ｔｈｅｐｏｏｌｅｄRRｖａｌｕｅｓｏｆｃｏｌｏｒｅｃｔａｌｃａｎｃｅｒｗｅｒｅ１（０．９３－１．１），１．３（１．１３－１．４），１．５（１．２６－１．９）ｆｏｒ２０，４０，６０ｇｒａｍ／ｄａｙ，ｒｅｓｐｅｃｔｉｖｅｌｙ；Ｔｈｅｌｉｎｅａｒｉｔｙｍｏｄｅｌｓｈｏｗｓｔｈａｔｔｈｅｒｉｓｋｏｆｃｏｌｏｒｅｃｔａｌｃａｎｃｅｒｉｎｃｒｅａｓｅｄ８％ｗｉｔｈａ１２－ｇ／ｄａｙｉｎｃｒｅａｓｅｆｏｒｔｈｅａｌｃｏｈｏｌｉｎｔａｋｅ．Ｔｈｅｒｅｗａｓｎｏｓｉｇｎｉｆｉｃａｎｔｄｉｆｆｅｒｅｎｃｅｂｅｔｗｅｅｎｔｈｅｌｉｎｅａｒａｎｄｎｏｎｌｉｎｅａｒｄｏｓｅ－ｒｅｓｐｏｎｓｅｃｕｒｖｅｓ．ConclusionＤｏｓｅ－ｒｅｓｐｏｎｓｅｍｅｔａ－ａｎａｌｙｓｉｓｃａｎｂｅｉｍｐｌｅｍｅｎｔｅｄｉｎＲｓｏｆｔｗａｒｅ，ａｎｄｉｔｉｓｐｏｗｅｒｆｕｌａｎｄｐｒａｃｔｉｃａｌ．Keywords:Ｄｏｓｅ－ｒｅｓｐｏｎｓｅ；Ｍｅｔａ－ａｎａｌｙｓｉｓ；Ｓｙｓｔｅｍｒｅｖｉｅｗ；Ｒｓｏｆｔｗａｒｅ

剂量反应关系的Ｍｅｔａ分析是一类新型的Ｍｅｔａ分析方法，相比标准的二分类及连续性资料Ｍｅｔａ分析，此类研究可同时处理三个及以上组别的数据，并直接估计暴露因素与疾病的剂量反应关系［1］。但此类研究的统计方法较为复杂，与标准Ｍｅｔａ分析有所不同，也无法用ＲｅｖｉｅｗＭａｎａｇｅｒ软件实现，需要借助Ｓｔａｔａ、ＳＡＳ等进行编程实现［2，3］。国内很少有关介绍剂量反应关系Ｍｅｔａ分析的软件实现方法的文献报道，这在一定程度上给这类新型Ｍｅｔａ分析方法的应用带来了障碍，２０１３年８月瑞典生物统计学家ＡｌｅｓｓｉｏＣｒｉｐｐａ［4］开发了一种基于Ｒ软件的剂量反应Ｍｅｔａ分析的软件包－“ｄｏｓｒｅｓｍｅｔａ”为这一类Ｍｅｔａ分析的实现提供了极大的便捷，本文结合实例在国内首次

作者简介：周权（１９８７－），男，硕士，研究方向：流行病与卫生统

计学

通讯作者：黎敏，Ｅ－ｍａｉｌ：ｍｉｋｅ３０５＠ｑｑ．ｃｏｍ

探讨Ｒ软件实现剂量反应Ｍｅｔａ分析的方法，以期为

此类Ｍｅｔａ分析实现提供方法学指导。1

Ｒ软件介绍及剂量反应Ｍｅｔａ程序包导入

Ｒ软件是一个自由、免费、源代码开放的软件，具有统计分析和统计制图的强大功能［5］。Ｒ软件的最新版本可以从ＣＲＡＮ（ｔｈｅＣｏｍｐｒｅｈｅｎｓｉｖｅＲＡｒｃｈｉｖｅＮｅｔｗｏｒｋ）下载，网址为：ｈｔｔｐ：／／ｃｒａｎ．ｒ－ｐｒｏｊｅｃｔ．ｏｒｇ／，截止到２０１４年６月２５日已最新至Ｒ－３．１．０版本。

Ｒ软件的剂量反应Ｍｅｔａ程序包“ｄｏｓｒｅｓｍｅｔａ”是由瑞典生物统计学家ＡｌｅｓｓｉｏＣｒｉｐｐａ开发，自２０１３年８月推出１．０版本后，之后经不断完善更新，于２０１４年１月１７日更新至１．３．０版。

Ｒ软件安装好后，先在电脑中的任意盘（如Ｄ盘）建一个名为Ｒｗｏｒｋｐｌａｃｅ文件夹，同时在桌面的Ｒ软件图标上单击右键选择“属性”，将“起始位置”修改为：Ｄ：＼Ｒｗｏｒｋｐｌａｃｅ，以后Ｒ软件运行的所

·７９０·

现代预防医学２０１５年第４２卷第５期ＭｏｄｅｒｎＰｒｅｖｅｎｔｉｖｅＭｅｄｉｃｉｎｅ，２０１５，Ｖｏｌ．４２，ＮＯ．５

有数据都会自动储存在这个文件夹中。运行Ｒ软件，（ｉｄ，字符型）、研究类型（ｔｙｐｅ，字符型，分别用点击“Ｐａｃｋａｇｅｓ”下拉菜单后，再选“ＳｅｔＣＲＡＮｃｃ、ｉｒ、ｃｉ表示病例对照研究、人年队列研究和累计

，会出现下载数据源的镜像列表，然后选择ｄｏｓｅ，数值型，一般取原ｍｉｒｒｏｒ”发病数队列研究）、剂量（离自己距离最近的镜像，如Ｃｈｉｎａ（Ｂｅｉｊｉｎｇ１）始暴 …… 此处隐藏：11622字，全部文档内容请下载后查看。喜欢就下载吧 ……