用户头像 |
2021年5月3日 20:44
求助生存资料诊断试验的R语言代码

我需要用于生存资料的多因素诊断试验的R语言代码——时间依赖的ROC曲线

如下是单因素生存资料的ROC曲线的R语言代码,其中mayoscore4是一个单变量,研究其和生存结局的关系。而我需要的是加入更多的变量比如年龄或性别等,进行多因素ROC曲线分析的R语言代码

library(survivalROC)

data(mayo)

cutoff <- 365

Mayo4.2= survivalROC(Stime=mayo$time,

           status=mayo$censor,

           marker = mayo$mayoscore4,

           predict.time = cutoff,

           method="KM")

Mayo4.2= survivalROC(Stime=mayo$time, 

status=mayo$censor, 

 marker = mayo$mayoscore4, 

predict.time = cutoff, method="KM")

plot(Mayo4.2$FP, Mayo4.2$TP, 

type="l",col="red",xlim=c(0,1), ylim=c(0,1), 

 xlab=paste( "FP", "\n", "AUC = ",round(Mayo4.2$AUC,3)), 

ylab="TP",

 main="Mayoscore 4, Method = KM \n Year = 1")

abline(0,1,col="gray",lty=2)

cut.op2=Mayo4.2$cut.values[which.max(Mayo4.2$TP-Mayo4.2$FP)]