R可视化：生存分析森林图

生存分析森林图

生存分析森林图

介绍

• ggplot2绘制生存分析森林图

加载R包

knitr::opts_chunk$set(message = FALSE, warning = FALSE)

library(tidyverse)
library(survival)
library(survminer)
library(tableone)
library(forestplot)

# rm(list = ls())
options(stringsAsFactors = F)
options(future.globals.maxSize = 10000 * 1024^2)

Importing data

data("colon")

inpudata <- colon

head(inpudata)

数据预处理

predata <- within(inpudata %>%
                     dplyr::filter(!is.na(differ)), {
  sex <- factor(sex, levels = c(0, 1), labels = c("female", "male"))
  differ <- factor(differ, levels = c(1, 2, 3), labels = c("well", "moderate", &#