R/TrioWPCCurve.R

In WPC: Weighted Predictiveness Curve

TrioWPCCurve <-
function(wpc1,wpc2,wpc3,xlab="Marker",ylab="Survival Rate",main="Weighted Predictiveness Curve",ylim=c(0,1),xlim=c(0,100),type="l",col1="red",
				col2="blue",col3="black",lwd=2,legendloc="bottomright",legendtxt=c("Method1","Method2","Method3"),confi="N",ptsest="N")
	{	   plot(wpc1$x, wpc1$s, xlab=xlab, ylab=ylab, main=main, ylim=ylim, xlim=xlim, type=type,col=col1,lwd=lwd)
		   points(wpc2$x, wpc2$s, type=type, col=col2, lwd=lwd)
		   points(wpc3$x, wpc3$s, type=type, col=col3, lwd=lwd)
		   if(confi=="N"){legend(legendloc,legendtxt,col = c(col1,col2,col3),bty="n",lty=c(1,1,1),lwd=c(lwd,lwd,lwd))}	
		   if(confi=="Y"){
 		   	points(wpc1$x,wpc1$lb,lwd=2,lty=2,col=col1,type="l")
		   	points(wpc1$x,wpc1$ub,lwd=2,lty=2,col=col1,type="l")	 
			color_transparent <- adjustcolor(col1, alpha.f = 0.1) 	
			polygon(c(wpc1$x, rev(wpc1$x)), c(wpc1$lb, rev(wpc1$ub)),col = color_transparent, border = NA)   
 		   	points(wpc2$x,wpc2$lb,lwd=2,lty=2,col=col2,type="l")
		   	points(wpc2$x,wpc2$ub,lwd=2,lty=2,col=col2,type="l") 	
			color_transparent <- adjustcolor(col2, alpha.f = 0.1) 	
			polygon(c(wpc2$x, rev(wpc2$x)), c(wpc2$lb, rev(wpc2$ub)),col = color_transparent, border = NA)  	   
		   	points(wpc3$x,wpc3$lb,lwd=2,lty=2,col=col3,type="l")
		   	points(wpc3$x,wpc3$ub,lwd=2,lty=2,col=col3,type="l")
			color_transparent <- adjustcolor(col3, alpha.f = 0.1) 	
			polygon(c(wpc3$x, rev(wpc3$x)), c(wpc3$lb, rev(wpc3$ub)),col = color_transparent, border = NA)  
		   	legend(legendloc,legendtxt,col = c(col1,col1,col2,col2,col3,col3),bty="n",lty=c(1,2,1,2,1,2),lwd=c(lwd,2,lwd,2,lwd,2))	}
		    if(ptsest=="Y"){
			points(wpc1$x,wpc1$y,lwd=3,col=col1)
			points(wpc2$x,wpc2$y,lwd=3,col=col2)
			points(wpc3$x,wpc3$y,lwd=3,col=col3)}}