生存分析森林圖

生存分析森林圖

介紹

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