案例
今天介绍的用于画森林图的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