利用survfit函数拟合生存曲线,ggsurvplot函数对生存分析结果可视化
代码直接放在这里啦↓↓↓
####1_加载包####
library(survival) #提供生存分析数据集、拟合生存曲线
library(survminer) #绘图可视化
####2_生存模型
#time生存时间,status状态,sex性别分组
fit <- survfit(Surv(time,status) ~ sex,data = lung)
ggsurvplot(fit)
ggsurvplot(fit,
pval = T, #p值
conf.int = T, #增加置信区间
conf.int.alpha = 0.1, #置信区间透明度
surv.median.line = "hv", #增加中位生存时间
palette = c("black","red"), #线的颜色
risk.table = T, #不同时间点风险人数表
xlab = "Survival time(days)", #x轴标题
legend.title = "Sex") #图例标题