我是荟萃分析的新手,我想知道是否可以自定义生成的森林图
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回复(正如我所需要的)完成:
答案 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来添加一个列,其中包含每个治疗的总研究数量,但是,这些数字的解释对我来说似乎很难。