生存分析森林图

生存分析森林图

介绍

• 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 <- colonhead(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", &#