R语言 森林图

案例

今天介绍的用于画森林图的R包为forestmodel，用到其中的函数forest_model，这个函数可以接收lm、glm、coxph等函数拟合的模型。我们先来画一个cox模型的森林图

library(forestmodel)
library(survival)


lung$ph.ecog <- factor(lung$ph.ecog)
lung$sex <- factor(lung$sex)
model <- coxph(Surv(time, status) ~ ., lung) #拟合模型
pic <- forest_model(model)
pic

效果图如下:

cox效果图展示

自定义展示内容

案例部分我们建立了一个生存分析模型，使用函数的默认设置得出了漂亮的森林图。forest_model这个函数的参数非常多，我认为比较重要的参数为panels。panels决定我们森林图展示的大概样子，有默认值default_forest_panels，我们也可以自己定义panels这个参数来定制自己的森林图。

下面重点介绍panels这个参数的设置（list with details of the panels that make up the plot ）。

可以将森林图想象为是由一系列的"列"拼起来的。我们需要通过panels这个参数定义我们的每一列是什么。如下面这张图，1是一个空列，2是我们的变量名，3是分类变量的level&#x