R语言快速绘制1、2年生存曲线并计算AUC

最新推荐文章于 2023-09-07 08:05:42 发布

weixin_49320263

最新推荐文章于 2023-09-07 08:05:42 发布

阅读量333

点赞数

文章标签： r语言

本文链接：https://blog.csdn.net/weixin_49320263/article/details/130469884

版权

library(timeROC)
library(survival)
library(ggplot2)
train12<-na.omit(lung)
table(train12$status)
moxing<-coxph(Surv(time,status)~age+sex,train12)
train12$marker<-predict(moxing)
ROC_1<-timeROC(T=train12$time,delta=train12$status,
               marker=train12$marker,cause=2,
               weighting="cox",
               times=c(365,730),ROC=TRUE)
ROC_1$FP_1
ggplot() +
  geom_line(aes(x=ROC_1[["FP_1"]][,1],y=ROC_1[["TP"]][,1],color="1 year"),size=1.2)+
  geom_line(aes(x=ROC_1[["FP_1"]][,2],y=ROC_1[["TP"]][,2],color="2 years"),size=1.2)+
  geom_abline(intercept = 0, slope = 1, linetype="dotted")+
  labs(title="ROC Curve", x="1 - Specificity", y="Sensitivity", color="Timepoint")+
  scale_color_manual(values=c("#009E73","#D55E00","#0072B2"))+
  annotate("text", x=0.8, y=0.2, label=paste0("AUC (1 year) = ",round(ROC_1$AUC_1,2)), color="black", size=3.5) +
  annotate("text", x=0.8, y=0.15, label=paste0("AUC (2 years) = ",round(ROC_1$AUC_2,2)), color="red", size=3.5)