森林图常用于展示对比多个实验结果的效应量和置信区间的差别,特别会在医学的Meta分析中经常用到。forestploter包可以绘制出简洁干净的森林图,并且由于结果是ggplot类,可以通过ggsave()函数输出结果,十分方便。下面我们通过一个假数据来展示如何使用forestploter
首先我们下载并加载forestploter包,并加载tidyverse包,代码如下:
然后我们需要通过下面的代码制作一个假数据。假数据仅作为数据整理结果的参考,并无实际意义,代码如下:
为了简化数据,以便绘图使用,我们需要再稍微整理一下,留下group、treatment、est、se、p.value、or值和相应的置信区间上下限,整理代码如下:
新生成了一列变量名是' '的字符变量,它下面的值全是同一长度的空格,主要目的是为了放置画出来的置信区间。
下面就可以开始绘制森林图了,首先我们不进行太多的细节设置,代码如下:
这是一个最简单的绘图的结果,我们还可以通过forest_theme()函数来设置图的背景,置信区间的细节等,代码如下:
以上就是使用forestploter绘制森林图的方法,只要数据已经被处理好,操作还是十分易于上手的。代码和输出的图片可以通过下面的链接获取:
链接:https://pan.baidu.com/s/1K8aaTICeIFfPUAmNmu1R_Q?pwd=82fe
提取码:82fe