自定义网络元分析的森林图

时间:2017-05-01 13:48:10

标签: r forestplot

我是荟萃分析的新手,我想知道是否可以自定义生成的森林图

data(Senn2013)
net1 <- netmeta(TE, seTE, treat1, treat2, studlab,
        data=Senn2013, sm="MD", reference="plac")

forest(net1, ref="plac", digits=1, just="right")

enter image description here

例如:i)去除安慰剂参考治疗; ii)更改结果数字; iii)添加一栏,列出每项治疗的研究数量; iv)按估计的效应大小对治疗进行分类。

我真的很感激这方面的帮助...

这个情节基于@Guido Schwarzer回复(正如我所需要的)完成:

enter image description here

1 个答案:

答案 0 :(得分:1)

自定义i) - iv)不能在CRAN上使用当前版本的netmeta。实际上,ii)数字应该有效(但没有),而i),iii)和iv)没有实现。

我在GitHub上的netmeta开发版中更改了forest.netmeta(),以便所有自定义工作。

## install.packages("devtools")
devtools::install_github("guido-s/meta")
devtools::install_github("guido-s/netmeta", ref = "develop")

library(netmeta)
data(Senn2013)
net1 <- netmeta(TE, seTE, treat1, treat2, studlab,
                data=Senn2013, sm="MD", reference="plac")

forest(net1, digits = 1, sortvar = TE,
       drop.reference.group = TRUE,
       leftcols = c("studlab", "k"),
       leftlabs = c("Contrast\nto Placebo", "Direct\nComparisons"),
       just.studlab = "right", just.addcols = "left")

安装GitHub版本的meta并不是绝对必要的,但修复了带有研究编号的列的列标签中的forest.meta()中的一个小错误(这是一个额外的专栏,以forest.meta为标准) ))。

请注意列&#34; k&#34;打印提供直接证据的研究数量,用于比较治疗与参考(此处:安慰剂),而不是治疗的总体研究数量。这似乎更适合我作为森林情节&#34; only&#34;显示治疗与参考比较的结果。原则上你可以使用新的参数add.data来添加一个列,其中包含每个治疗的总研究数量,但是,这些数字的解释对我来说似乎很难。