R - 如何使用 metafor 包绘制森林图？

【问题标题】：R - How to plot forest plot using metafor package?R - 如何使用 metafor 包绘制森林图？
【发布时间】：2020-05-16 03:37:22
【问题描述】：

我想使用 metafor 包从模型的统计摘要中保存的对象绘制森林图。在我的对象中，将有 5 列：组、子组、估计、上限和下限 CI，如下所示：

Group Subgroup  est   lCI hCI
A         1     0.2    ~   ~
A         2     0.3    ~   ~
A         3     0.8    ~   ~
B         1      ~     ~   ~
B         2      ~     ~   ~
B         3      ~     ~   ~

如何绘制将子组组合在一起并进行分类的森林帖子？以 CI 作为误差线。

谢谢！

【问题讨论】：

标签： r plot model summary metafor

【解决方案1】：

我之前只使用过forestplot 包，您需要先将其旋转，然后将其绘制如下：

library(forestplot)
library(tidyr)
library(dplyr)

results = data.frame(Group=rep(c("A","B"),each=3),
Subgroup = rep(1:3,2),est = runif(6,min=-2,max=2))

results$lCI = results$est - 0.1
results$hCI = results$est + 0.1

df_wide = pivot_wider(results,id_cols=Group,names_from=Subgroup,values_from=c("est","lCI","hCI"))

forestplot(list(Var=df_wide[["Group"]]), 
           legend = 1:3,
           fn.ci_norm = c(fpDrawNormalCI, fpDrawCircleCI),
           mean = select(df_wide,contains("est")),
           lower = select(df_wide,contains("lCI")),
           upper = select(df_wide,contains("hCI")),
           col=fpColors(box=c("blue", "darkred"))
            )

【讨论】：