r语言 随机森林_宇宙最全R语言科研统计/作图代码森林图

小编又来了~~！

R语言科研统计/作图第二期    森林图

PS：上一期推出之后，有不少小伙伴建议小编讲讲使用场景。那我们就从这一期开始，小编会在开始的部分以对话的形式向大家简单介绍使用场景。

因为太多忠实funs，小编觉得后面每一期都直接上中阶版，回馈大家的厚爱。希望大家帮忙分享/转发给你们身边有需要的人，谢谢！

Let's Go!

---

读者：什么是森林图？该什么时候用呢？

小编： 森林图常见于元分析，但其使用绝不仅如此，比如你现在想要研究的对象有诸多HR结果，表格太复杂，想要汇总为一张图，又要一目了然，森林图就是个非常好的选择。

---

今日(2019-9-1)----森林图绘制

??????

项目	基础版 (2个分组变量)	中阶版 (2个分组变量)	高级版 (多个分组变量)	懒人版
生存曲线样图
获取方式	原代码	操作文档+原代码	视频+操作文档+原代码	点击详解

上干货

??????

原代码

library(forestplot)

#数据读入及检查-----------------------------------------------------

mydata

#矩阵建立------------------------------------------------------------

aas.matrix(mydata[,1:3])

#绘制森林图------------------------------------------------------------

forestplot(a,mydata$HR,mydata$LowerCI, mydata$UpperCI, 

           graph.pos = 2, zero=1, 

           graphwidth = "auto", lineheight = "auto", 

           boxsize = 0.2, 

           col = fpColors(all.elements = "black"))

---

教程详解

??????

第一部分，数据准备，森林图用的数据需要我们前期处理

数据包括以下几个部分：6列

变量，显示组合(HR [95%CI区间])，P值，HR值，95%CI区间下界，95%CI区间上界

详细的表格见下图所示

第二部分，数据读取和准备

#数据读入及检查-----------------------------------------------------

mydata

#矩阵建立------------------------------------------------------------

a

PS:小编特别提醒下大家，矩阵的主要目标把数据变成文本，否则部分数据会因R语言中对数据的特定格式要求而受限，比“ ”， “-”，“/”都部分出现或者出现报错！请各位特别留心和注意，谢谢

第二部分，深林图表的绘制

#绘制森林图------------------------------------------------------------

#a是矩阵文本,mydata$HR/LowerCI/UpperCI，都是列表中对应的列

forestplot(a,mydata$HR,mydata$LowerCI, mydata$UpperCI, 

#graph.pos表示图位于哪一个部分，Zero表示切割线在哪，一般默认为1。

           graph.pos = 2, zero=1, 

#graphwidth表的宽度，可以自己设定，也可以选择“自动”

           graphwidth = "auto", lineheight = "auto", 

#boxsize图标的大小，建议0.1-0.4之间，主要看各位的变量数量，大家可以根据实际调整。

           boxsize = 0.2, 

#col图,图标的颜色，小编比较偏爱黑，简单明了，大家也可以更具自己喜欢修改。

           col = fpColors(all.elements = "black"))

这样图表就完成了，看看效果图吧

读者：小编，小编，我想要多样的，变量的，和你一样技能满满?

小编：转发朋友圈，截图+回复”SLT" , 小编立马为你送上高阶版教程?

想学习R语言的，不要错过我们的公众号哦~9月份起，每3天一期持续更新

点击原文，关注公众号，查看所有R语言教程