plot.cp=function(ab,reverse=F,xlab="Time to event",
ylab="Conditional probability")
{
#######################################
## plots the conditional probability  #
## aa is an object created with Cpvar #
#######################################
aa=ab
if (is.data.frame(ab)) aa=list(ab)

ng=length(aa)
m=0

for (i in 1:ng)
{
m=max(m,aa[[i]]$time)
}

for (i in 1:ng)
{
doubletime=sort(c(aa[[i]]$time,aa[[i]]$time))
doublecp=sort(c(aa[[i]]$cp,aa[[i]]$cp))
doubletime=c(0,doubletime)
doublecp=c(0,0,doublecp[1:(length(doublecp)-1)])
if (reverse) doublecp=1-doublecp
if (i==1)
{plot(doubletime,doublecp,xlim=c(0,m),ylim=c(0,1),
xlab=xlab,ylab=ylab,type="l",lty=1)}
else {lines(doubletime,doublecp,lty=i)}
}
}