R/plot.ipcwsurvivalROC.R

In timeROC: Time-Dependent ROC Curve and AUC for Censored Survival Data

plot.ipcwsurvivalROC <- function(x,time,col="red",add=FALSE,title=TRUE,...){
  if (missing(time)){
    stop("time is missing. \n") }
  if (length(time)!=1){
    stop("time must be an unique value. \n It indicates the time point at which the ROC curve is plotted. \n") }
  if ( ! time %in% x$times){
    stop("time must be one of the value included in the vector times that you used  for estimation.\n") }  
  if (add==FALSE) {
    plot(0,0,type="l",xlim=c(0,1),ylim=c(0,1),xlab="1-Specificity",ylab="Sensitivity")
    if(title==TRUE){
      title(paste("ROC at time t=",time,", AUC=",round(x$AUC[which(x$times==time)]*100,1),sep=""))
    }
  }
  lines(x$FP[,which(x$times==time)],x$TP[,which(x$times==time)],col=col,type="l",...)
  abline(0,1,lty=2)
}