修改 ggsurvplot 子图

[ShixiangWang]

备份一下问题和答案。

#setwd("C:\\Users\\89793\\Desktop\\数据")

library("survminer")
PA=read.table("数据.txt",sep="\t",header=T,check.names=F)
require("survival")
fit <- survfit(Surv(days, status) ~ group, data = PA)

pdf(file="OS.pdf",width=10,height=10, onefile = FALSE)

ggsurv <- ggsurvplot( 
fit, # survfit object with calculated statistics. 
data = PA, # data used to fit survival curves. 
risk.table = TRUE, # show risk 
pval = TRUE, # show p-value of log-rank test. 
conf.int = TRUE, # show confidence intervals for  
                           # point estimates of survi
palette = c("#E7B800", "#2E9FDF"), 
xlim = c(0,150), # present narrower X axis, but not affect
                         # survival estimates. 
xlab = "Time in months", # customize X axis label. 
break.time.by = 15, # break X axis in time intervals by 500. 
ggtheme = theme_light(), # customize plot and risk table with a theme. 
risk.table.y.text.col = T,# colour risk table text annotations. 
risk.table.height = 0.20, # the height of the risk table 
risk.table.y.text = FALSE,# show bars instead of names in text annotations 
                                        # in legend of risk table. 
ncensor.plot = TRUE, # plot the number of censored subjects at time t 

ncensor.plot.height = 0.20, 
conf.int.style = "step", # customize style of confidence intervals 
surv.median.line = "hv", # add the median survival pointer. 
legend.labs = c("male", "female") # change legend labels. 
)


# Changing Labels
# %%%%%%%%%%%%%%%%%%%%%%%%%%
ggsurv$plot <- ggsurv$plot + labs(
  title    = "survival curves",                     
  subtitle = "Based on Kaplan-Meier estimates"  
  )

# Labels for Risk Table 
ggsurv$table <- ggsurv$table + labs(
  title    = "Note the risk set sizes")

# Labels for ncensor plot 
ggsurv$ncensor.plot <- ggsurv$ncensor.plot + labs( 
  title    = "Number of censorings")

# 这些图每部分都与ggplot相关，所以可以通过ggplot方式修改
# class(ggsurv$ncensor.plot)
# ggsurv$ncensor.plot + scale_y_continuous(breaks = seq(0,30,5))
ggsurv$ncensor.plot <- ggsurv$ncensor.plot + scale_y_continuous(breaks = seq(0,30,5))

ggsurv
dev.off()