手机下载视频网站模板,七牛云wordpress缓存配置,济南企业免费建站,做网站坚持原创文章有什么好处近来年#xff0c;网状 Meta 分析相关研究不断涌现#xff0c;此类研究不但能发表在国内各大核心期刊上#xff0c;还能在SCI期刊甚至医学4大刊上看到其身影。随手在pubmed上面一搜索#xff0c;就能得到一万多篇相关文献。俨然成为医学文献研究的“大杀器”#xff01; P…
近来年网状 Meta 分析相关研究不断涌现此类研究不但能发表在国内各大核心期刊上还能在SCI期刊甚至医学4大刊上看到其身影。随手在pubmed上面一搜索就能得到一万多篇相关文献。俨然成为医学文献研究的“大杀器” PubMed网状 Meta 分析检索截图
本文就网状 Meta 分析的概述、用途及操作步骤R语言方面进行了介绍。篇幅会比较长内容较多请多点耐心噢。 1、Meta分析简介
Meta分析(Meta-analysis)常译为汇总分析或荟萃分析是系统评价(systematic review,SR)中经常用到的统计方法是运用恰当的统计方法对若干个研究的结果数据进行合并在医学上常常用于比较干预措施的有效性和安全性。
2、常规Meta分析的不足
常规的meta分析通常只能处理直接比较比如下图之中的文献评价了虚拟现实VR技术对乳腺癌患者术后上肢功能康复的影响其纳入的5项RCT随机对照试验均为两组设计的研究都是VR干预与常规干预的比较。 知网截图
3、网状 Meta 分析
但是如果说想比较干预A与干预B的效果但仅有比较A与CB与C的研究那么怎么办呢此时就需要网状 Meta 分析了
NMA是什么 普通Meta分析是合并直接比较的结果比如你想研究A药和B药对C疾病的治疗效果如果有A药和B药直接对比的临床研究那么你就可以做普通的Meta分析直接合并这些临床研究得到A药和B药的对比结果。 但是很多时候会缺少A药和B药直接对比的临床研究但是有A药和C药直接对比以及B药和C药直接对比的临床研究这样C药就是共同对照以C药为中介我们可以得到A药和B药间接比较的结果。 知网截图
上面的文献运用了网状 Meta 分析评价了10 种非药物干预方法在减少孕产妇分娩恐惧方面的效果并得出了干预效果排序。这就是网状Meta分析最大的特色。
RCT研究1A药可让患者血压降低10mmHgC药可让患者血压降低2mmHg所以A药比C药能多降低8mmHgdA-C
RCT研究2B药可让患者血压降低15mmHgC药可让患者血压降低5mmHg所以B药比C药能多降低10mmHgdB-C
C是A药和B药的共同对照dA-C-dB-C-2mmHg所以B药相对于A药患者平均血压能多降低2mmHg。这里举的是连续性变量如果是二分类变量OR值就是取对数相减。 上面说的是间接比较的概念有的人以为NMA就是间接比较那就错了实际上有的时候既有直接比较的文献又有间接比较的文献这时候NMA是综合直接比较和间接比较的结果。就如下面的图。 NMA的定义就是基于多个研究分析两个以上干预措施之间间接比较结果或者直接比较结果与间接比较结果的合并结果混合比较的Meta分析。
那么什么情况下适合做NMA呢
主要就是要满足三个基本假设同质性假设、相似性假设和一致性假设
同质性假设其实就是普通Meta分析中的异质性原始研究之间临床特征、方法学、统计学没有异质性。
相似性假设就是所有原始研究间及不同对照组间影响效应量的因素相似。大概意思就是A和C比较的原始研究要和B和C比较的原始研究之间是相似的
一致性假设假如A和B之间的比较既有直接比较也有间接比较那么直接比较的结果和间接比较的结果要一致。
那么NMA到底怎么做呢
其实跟普通的Meta分析步骤区别并不是很大只不过展示的结果会多一些除了直接比较的内容你还要展示间接比较的内容。
基本步骤
1.选题
2.文献检索
3.数据提取
4.质量评价
5.统计分析比普通Meta分析复杂
6.报告结果间接比较或混合比较
4、贝叶斯meta分析
关于贝叶斯内容实在太多个人仅了解了一点皮毛仅作简单描述有兴趣的童鞋可以去看看相关统计教材。
meta分析方法主要包括频率学法和贝叶斯法这是因为当前统计学可以分为两个主要学派经典统计学与贝叶斯统计学派。
频率学方法基于概率的频率解释将概率解释为大量重复试验后频率的稳定值贝叶斯法认为某一事件的发生概率不是固定的而是基于观察数据及一些条件来对先验概率进行修正并对所有分析的干预措施进行排序。
5、适用场景 对多种干预措施进行比较 无直接比较的原始研究 有直接比较的原始研究但数量或质量不理想
6、结果报告与解读
报告结果包括网状关系图一致性检测构建模型排序图异质性分析判断发表偏倚等。这里以李腾等学者发表的论文《非药物干预减少分娩恐惧的网状 Meta 分析》为例了解下网状Meta分析通常应输出的结果。
①网状关系图与一致性检验
网状关系图图中的节点表示不同的干预措施连线表示两种干预方法有直接比较节点越大表示该干预措施直接比较出现的次数越多样本最大形成的闭环表示既有直接比较又有间接比较。
一致性检验是比较直接比较与混合比较两个模型的一致性P0.05表示一致性良好。 图片来源《非药物干预减少分娩恐惧的网状 Meta 分析》
②排序图/表
排序图/表的形式多样采用不同软件也会有差异。这里附上2份排序表的截图。 图片来源《非药物干预减少分娩恐惧的网状 Meta 分析》不同颜色代表不同效果 图片来源《药物治疗不宁腿综合征有效性和安全性的网状 Meta 分析》最下方的一行疗效最好 图片来源《药物治疗不宁腿综合征有效性和安全性的网状 Meta 分析》看整体排序排序为1表示最好2次之。
其余的结果报告与普通的meta分析类似限于篇幅这里就不一一描述了。 网状Meta分析实操——基于R语言的实现
Stata、R语言、WinBUGS是贝叶斯网状Meta分析应用最多的工具。Stata收费这里以开源的R语言进行实操分享。以下的步骤参照了向虹等学者发表的论文《应用R语言BUGSnet程序包实现贝叶斯网状Meta分析》个人把此论文的代码跑了一遍然后加了一些便于理解的注释顺便分享了一些可能会碰到的问题可能踩的坑。
文章后面附上了我自己录入的数据和代码文件的获取方式。代码经个人测试与论文的结果保持一致。下面是详细的步骤。
一 、准备数据
数据来源于一项关于类固醇佐剂治疗天疱疮的网状Meta分析研究该研究纳入 10 个试验纳入了592 例患者评估了 7 种方法治疗天疱疮的效果见下表。 图片来源于论文《应用R语言BUGSnet程序包实现贝叶斯网状Meta分析》
这里有个坑按照这个表格录入数据会有变量名不匹配的问题要求是英文且为了与参照的论文中的代码保持一致还是要注意下。PS其实自己随便取个英文名也是可以的只是说这样的话后续的代码不能直接用了要做针对性地调整。
以下是我手动录入的表格数据建议严格按照下图的变量名即表格中的列名进行设置最后一行因截屏问题没有显示出来大家录入的时候补上就好。 数据文件截图
变量名我给附上来了练习时建议保持一致可以直接copy。
StudyTreatmentNumber.of.patientsNumber.of.remissionsNumber.of.relapsesNumber.of.adverse.event.related.withdrawalsCumulative.dose.of.steroid_MeanCumulative.dose.of.steroid_SD
变量名
二、安装软件和包
如果是第一次接触R语言请看下面的安装过程依次安装3个软件如果已经有基础或已安装建议跳过。 R语言我用的R语言版本为当前的最新版R4.2.2直接到官网下载https://www.r-project.org。下载后默认设置进行安装即可也可以手动改下安装的路径不熟悉的可以去百度一下具体方法。 R StudioR Studio是R语言的IDE很好用我们基本就是用它来写R代码直接到官网下载https://www.rstudio.com/products/rstudio/。下载完安装。 JAGS-4.3.0 R语言BUGSnet包实际就是来调用JAGS软件来实现贝叶斯网状meta分析的为何不直接用JAGS那是因为R语言可以提供很漂亮的绘图这个软件下载很慢文末有软件获取方法这里特别感谢小师妹分享的JAGS软件包。 成功安装后打开R Studio软件萌新可以先百度学一下基本的使用方法。这里建议大家做一些必要的设置
1、将镜像源更改为国内
步骤Tools - GlobalOptions-Options - Packages 单击change可修改成chinakaitou镜像源这个操作可以极大地提升安装包的速度。
2、设置code编码
Tools - Global Options-Default- code - Saving text - encoding框单击change将相关设置修改为UTF-8。有编程基础的童鞋应该都清楚utf-8是国际通用的编码可以避免换一台电脑打开就乱码的困扰。 在R Studio中安装一些包
注意#是注释的意思是给我们看的说明文字不会被执行。代码如下。
#安装 Rtoolinstall.packages(pkgbuild)pkgbuild::has_build_tools()#安装BUGSnet包install.packages(c(remotes, knitr)) #为了实现下一步在github中安装包remotes::install_github(audrey-b/BUGSnetv1.1.0, upgrade TRUE, build_vignettes TRUE, dependencies TRUE) #从github中安装BUGSnet 1.1.0#没有安装下面这些包的请执行下面代码进行安装install.packages(c(xlsx, dplyr, tidyr)) #以上是原文的代码目的是为了导入excel格式的文档#由于我这里安装时报错为了简单最后换用了基础包下的read.csv无须安装 三 、网状meta分析具体实现过程
1、录入、整理数据
library(dplyr)#导入dplyr包library(tidyr)#导入tidyr包目的是实现管道操作及一系列tidy风格的函数
data_total-read.csv(meta.csv, 1) #导入编辑好的数据文件,csv格式meta.csv是文件名称str(data_total) #查看导入的数据信息data_remission- data_total[, c(1:4)] #拿到第1~4列并命名为data_remissiondata_relapse-data_total[, c(1:3, 5)] #拿到第1~3列和第5列并重命名data_withdrawal- data_total[, c(1:3, 6)]#拿到第1~3列和第6列并重命名data_dose-data_total[, c(1:3, 7:8)]#拿到第1~3列和第7~8列并重命名
#下面代码重命名结局指标#data_remission代表缓解病例数据为主要疗效指标data_remission-data_remission %% rename(sampleSizeNumber.of.patients, events Number.of.remissions) #将data_remission内的变量重命名然后将整个数据依然命名为data_remission
# data_relapse代表复发病例数据为次要疗效指标data_relapse-data_relapse %% drop_na() %% rename(sampleSizeNumber.of.patients, events Number.of.relapses) #重命名
# data_withdrawal代表不良事件相关退出病例数据为次要安全性指标data_withdrawal-data_withdrawal %% drop_na() %% rename(sampleSize Number.of.patients, eventsNumber.of.adverse.event.related.withdrawals)#重命名
#data_dose代表固醇的累积使用剂量为主要安全性指标data_dose-data_dose %% drop_na() %% rename(sampleSizeNumber.of.patients, mean Cumulative.dose.of.steroid_Mean, SDCumulative.dose.of.steroid_SD)#重命名 2、绘制网络关系图
library(BUGSnet) #导入BUGSnet包#预处理指定目标数据及变量名data_remission - data.prep(arm.data data_remission, varname.t Treatment, varname.s Study) data_dose - data.prep(arm.data data_dose, varname.t Treatment, varname.s Study)#绘制描述性统计图网络关系图par(mfrow c(1,2)) #在同一张图上显示1x2个图net.plot(data_remission, node.colour darkgrey) #绘制散点图连线颜色为darkgreynet.plot(data_dose, node.colour darkgrey) 输出的网络关系图
3、设置模型类型
#使用nma.model函数设置拟合哪种模型(固定效应模型或随机效应模型)#设置固定效应模型fixed_effects_model_remission - nma.model(data_remission, outcomeevents, NsampleSize, referenceSteroid alone, familybinomial,linklogit, effectsfixed)
fixed_effects_model_dose - nma.model(datadata_dose, outcomemean, NsampleSize,sdSD, referenceSteroid alone, familynormal,linkidentity, effectsfixed) #nma.model参数解释#outcome指定结局变量# reference设置对照组本例为Steroid alone单独使用类固醇治疗# family用于指定结果的分布类型binomial二分类资料normal计量资料poisson计数资料# link用于指定NMA模型使用的函数logit用于ORlog用于RR或计数数据的比率Rate Ratio# cloglog用于二分类数据危险比HRidentity用于计量资料# effects用于设置效应模型的类型。
#设置随机效应模型random_effects_model_remission - nma.model(data_remission, outcomeevents, NsampleSize, referenceSteroid alone, familybinomial, linklogit, effectsrandom) random_effects_model_dose - nma.model(datadata_dose, outcomemean, NsampleSize, sdSD, referenceSteroid alone, familynormal, linkidentity, effectsrandom) 4、执行贝叶斯网状Meta分析
# 使用nma.run函数执行贝叶斯网状Meta分析set.seed(20222022) #随机数字生成器fixed_effects_results_remission - nma.run(fixed_effects_model_remission, n.adapt1000, n.burnin1000, n.iter10000)
fixed_effects_results_dose - nma.run(fixed_effects_model_dose, n.adapt1000, n.burnin1000, n.iter10000) random_effects_results_remission - nma.run(random_effects_model_remission, n.adapt1000, n.burnin1000, n.iter10000)
random_effects_results_dose - nma.run(random_effects_model_dose, n.adapt1000, n.burnin1000, n.iter10000)
#评估效应模型并输出杠杆图par(mfrow c(2,2)) nma.fit(fixed_effects_results_remission)nma.fit(fixed_effects_results_dose)nma.fit(random_effects_results_remission)nma.fit(random_effects_results_dose) 代码输出的图
①绘制排名图SUCRA 和 rankogram
#使用nma.rank函数 绘制治疗排名结果图largerbetter表示更大更好sucra_out_remission- nma.rank(random_effects_results_remission, largerbetterTRUE)sucra_out_dose- nma.rank(fixed_effects_results_dose, largerbetterFALSE)
# 使用ggplot2绘图神奇进行绘图install.packages(ggplot2) #安装包如果已有请忽略library(ggplot2) #导入包
#设置颜色z-scale_colour_manual(valuesc(AZA#FF6600, MMF#FF0088, Steroid alone#00FF00, CP_P #007700, Cyclosporine#6633CC, DCP_C(6M) #00CCFF, DCP_C(12M)#FFFF00, Rituximab #AA0000)) #绘制累积排序概率图下面积SUCRAs1- sucra_out_remission$sucraplotzs4- sucra_out_dose$sucraplotz
#绘制排序概率图rankogramr1- sucra_out_remission$rankogramtheme(axis.text.x element_text(angle 30, hjust 1, vjust 1))# theme参数是设置主题包括文字等r4- sucra_out_dose$rankogramtheme(axis.text.x element_text(angle 30, hjust 1, vjust 1))
#输出图片install.packages(gridExtra)#安装包library(gridExtra) #导包grid.arrange(s1, s4, r1, r4, ncol 2) #2x2SUCRA 图和 rankogram 图解读前两幅图是越高越好后两幅颜色越淡越好总结Rituximab疗效最好 ②绘制排名表热图
# nma.league() 函数用于获得排名表热图 league.out_remission- nma.league(random_effects_results_remission, central.tdcymedian, order sucra_out_remission$order, log.scale TRUE, low.colour springgreen4, mid.colour white,high.colour red, digits 2)
#参数解释#central.tdcy用于设置统计数据可以设置“mean”或“median”#order用来排序#log.scaleTRUE 表示用对数形式显示数据#colour用于设置颜色本例低、中、高设置了3种不同颜色#digits设定小数点后的显示位数league.out_dose- nma.league(fixed_effects_results_dose, central.tdcymedian, order sucra_out_dose$order, log.scale TRUE, low.colour springgreen4, mid.colour white, high.colour red, digits 2)
#绘制排名表热图l1-league.out_remission$heatplotl4-league.out_dose$heatplotgrid.arrange(l1, l4, ncol 1) #图片放在一列上输出
#可以将排名表热图作为表格输出league.out_remission$tableleague.out_remission$longtable 排名表热图解读不同颜色代表不同的效果
③ 绘制森林图
# 使用nma.forest() 函数绘制 森林图#参数解释与排名表热图一样f1- nma.forest(random_effects_results_remission, central.tdcymedian, log.scale TRUE, comparator Steroid alone)f4- nma.forest(random_effects_results_remission, central.tdcymedian, log.scale TRUE,comparator Steroid alone)grid.arrange(f1, f4, ncol 2) #输出森林图输出的森林图解读Rituximab是最有效的疾病缓解干预措施 关于实操部分的代码是可以直接复制下来的大家也可以按照推文中参考文献手动输入代码。如果觉得手动敲代码太累可以关注本公众号并回复关键词网状meta就可以直接获取数据文件及代码文件的链接。后续会继续更新有关内容
[1]向虹,岳磊,赵红,李华,谢天宏,杨婷,王杰,张宇昊,谢忠平.应用R语言BUGSnet程序包实现贝叶斯网状Meta分析[J].中国循证医学杂志,2022,22(05):600-608.
[2]张茜,陈项婷,周长青.药物治疗不宁腿综合征有效性和安全性的网状Meta分析[J].中国循证医学杂志,2022,22(08):908-916.
[3]李腾,张永爱,张海苗,曾娟.非药物干预减少分娩恐惧的网状Meta分析[J].中国循证医学杂志,2022,22(12):1410-1418.
