用R进行meta分析（metafor包）

这里以我的一篇meta分析为例，详细描述meta分析的一般步骤，该例子实现的是效应量β的合并

R包:metafor或meta包，第一个例子以metafor包为例

1.准备数据集

2.异质性检验

install.packages("metafor")
library(metafor)
metamod<-rma(yi=β,data=metabirth3,sei=se,method="DL")
#β是回归系数的值，se是标准误，DL表示选择的是随机效应模型，固定效应模型是method="FE"
summary(metamod)
#输出结果

可以看到I^2为28.34%，属于低度异质性，可采用固定效应模型。异质性低的时候可以采用固定效应模型和随机效应模型，结果差别不大，但高异质性只能选择随机效应模型，否则会使结果外推性受到约束。此处选择随机效应模型是出于保守情况考虑。

fixed-effect model是基于跨研究间不存在异质性的假设，也就是其合并时，默认运算是认为所纳入合并的研究均为同质；
random-effect model是基于跨研究间存在异质性的假设，该合并模型承认研究间异质性的存在，但是不对异质性加以处理；
如果纳入合并的研究间存在异质性，尽管未达到我们常规设定的I^2>50%，但是在用fixed-effect model合并时，默认运算直接忽略这一部分异质性的存在，这样合并的结果会造成假阳性误差，而选用random-effect model合并时，尽管不处理异质性，但是其默认运算承认异质性的存在，合并结果更可信！

3.森林图绘制

forestplot<-forest(metamod,refline = 1,mlab="Random-effect Model for All Studies",slab=paste(metabirth3$author,metabirth3$year,sep=","),xlab="β",showweights = T)
text(-500,6:1,pos=2,metabirth3$country)
text(c(-1600,-500,300,800),8,pos=c(4,2,4,4),c("Author(s) and Year", "Location","Weight","β[95%CI]"),cex=1,font=2)#添加标签

结果如下图

黑色方块表示的是权重大小，对应着右方的weight，水平横线是95%置信区间。菱形是合并效应值，95%置信区间小于0，意味着meta分析结果又统计学意义。

4.文献发表偏倚检验

可以用漏斗图及Begger's及Egger's检验

funnel(metamod)
ranktest(metamod)#Begg's检验#
regtest(metamod)#Egger's检验#

可能是由于文献数量的限制，漏斗图呈现明显不对称模式，无法判断是由于发表偏倚所致还是文献数量太少的缘故

可以看出Begg's检验及Egger's 检验的结果，P值都是大于0.05的，也就意味着没有发表偏倚。

5.敏感性检验

leavelout(metamod,digits = 3)

用R进行meta分析（metafor包）相关推荐

熟练掌握R语言的Meta分析全流程和不确定性分析，并结合机器学习等方法讲解Meta分析在文献大数据的延伸应用
Meta分析是针对某一科研问题,根据明确的搜索策略.选择筛选文献标准.采用严格的评价方法,对来源不同的研究成果进行收集.合并及定量统计分析的方法,最早出现于"循证医学",现已广泛应 ...
基于R语言的Meta分析【全流程、不确定性分析】方法与Meta机器学习技术应用
Meta分析是针对某一科研问题,根据明确的搜索策略.选择筛选文献标准.采用严格的评价方法,对来源不同的研究成果进行收集.合并及定量统计分析的方法,最早出现于"循证医学",现已广泛应 ...
基于R语言的Meta分析【全流程、不确定性分析】方法与Meta机器学习应用
Meta分析是针对某一科研问题,根据明确的搜索策略.选择筛选文献标准.采用严格的评价方法,对来源不同的研究成果进行收集.合并及定量统计分析的方法,最早出现于"循证医学",现已广泛应 ...
r语言找不到cochrane函数_网状meta必备技能之6-利用R中的meta包实现meta分析
我们仅仅是代码的编辑者.整合者.搬运工,仅免费传授方法,下文数据和代码取自于网络和免费软件"R语言说明书",如果您觉得我们侵犯了您的版权,请通知我们撤稿.请大家谅解,谢谢! 相信大 ...
用R进行gwas meta分析,原来如此简单
欢迎关注"生信修炼手册"! 在生物信息数据分析中,R语言是必备技能,简洁的语法,丰富的生态,美观的可视化,种种优势使得其成为该领域中使用最广泛的编程语言之一.用R进行meta分析当 ...
统计语言类教程：R语言贝叶斯统计学、Copula、SEM、极值统计学、混合效应模型、R\Python\matlab机器学习、科研数据可视化、线性回归、分位数回归、GAMS、meta分析、近红外光谱等.
查看原文>>>统计语言类教程:贝叶斯统计学.Copula.SEM.极值统计学.混合效应模型.PyTorch深度学习.科研数据可视化以下给大家整理了一些常用的统计学内容和python ...
基于R语言 meta/metafor 包进行荟萃分析（一）
目录目录 1.引言 2.安装,数据的读取和保存不讲.这个有别的大佬讲了,没有创新点. 3.标准方法 3.1 固定效应模型与随机效应模型 3.1.1 连续型结果的数据测量 3.1.2 固定效应模型 ...
贝叶斯软件genle教程_手把手教你用R的gemtc包对生存数据进行贝叶斯网状Meta分析...
大家好,本教程将介绍如何使用R的gemtc包对生存数据(HR为效应量)进行贝叶斯网状Meta分析. 前提条件: 需要下载R软件(推荐使用的R版本为3.5.3),以及RStudio(一个R的友好交互界面 ...
r语言和metawin_手把手教你用R的gemtc包对连续变量进行贝叶斯网状Meta分析
来源:微信公众号小侃数据大家好,本教程将介绍如何使用R的gemtc包对连续变量进行贝叶斯网状Meta分析. 前提条件: 准备工作完成后,打开RStudio后,在console里输入install.p ...
细思极恐——R语言forestplot包画meta分析群体药动学常用森林图
细思极恐--R语言forestplot包画meta分析群体药动学常用森林图今天,笔者想分享一下最近科研作图的经历,最主要的就是用于群体药动学模型建立的森林图,其百度百科定义为: 森林图是以统计指标和 ...

用R进行meta分析（metafor包）

用R进行meta分析（metafor包）相关推荐

最新文章

热门文章