我正在使用metafor包进行元分析。由于我的数据集包含了许多研究,因此森林情节变得非常长,我几乎无法阅读。
我有什么方法可以分割森林情节吗?
谢谢。
答案 0 :(得分:1)
您可以调整输出并绘制子集图,如下所示:
library(metafor)
set.seed(5132)
k <- 50
vi <- rchisq(k, df=1) * .03
yi <- rnorm(k, rnorm(k, 0.5, 0.4), sqrt(vi))
res <- rma(yi, vi)
# plot of studies 1:50
forest(res)
res2 <- res
res2$vi.f <- res2$vi.f[1:25]
# plot of studies 1:25
forest(res2, xlim = c(-3, 4))
res3 <- res
res3$vi.f <- res3$vi.f[26:50]
res3$yi.f <- res3$yi.f[26:50]
res3$slab <- res3$slab[26:50]
# plot of studies 26:50
forest(res3, xlim = c(-3, 4))