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

Question

我想使用metafor包从模型的统计摘要中保存的对象来绘制森林图。在我的对象中，将有5列: Group、sub groups、estimate、上限和下限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作为误差栏。

谢谢!

Answer 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"))
            )

​

​