为网络元分析定制森林图

Question

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

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

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

例如：i)去除安慰剂参考治疗； ii) 改变结果位数； iii) 添加一列，其中包含每种治疗的研究数量； iv) 按估计的效果大小对处理进行排序。

非常感谢在这方面的帮助...

此图根据@Guido Schwarzer 回复完成（正是我所需要的）：

Answer 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()).

行话专栏

请注意，"k" 列打印的是为治疗与参考（此处：安慰剂）的比较提供直接证据的研究数量，而不是治疗研究的总数。这对我来说似乎更合适，因为森林图 "only" 显示了与参考进行治疗比较的结果。原则上，您可以使用新参数 add.data 添加一列，其中包含每种治疗的研究总数，但是，对这些数字的解释对我来说似乎很困难。

为网络元分析定制森林图

Customizing forest plot for a network meta-analysis

r

forestplot