data(dataHepatology)
times <- dataHepatology$time
survival <- dataHepatology$survival
study <- dataHepatology$study
plot(times, survival, type="n",
ylim=c(0,1), xlab="Time",ylab="Survival")
for (i in unique(sort(study)))
{
lines(times[study==i], survival[study==i], type="l", col="grey")
points(max(times[study==i]),
survival[study==i & times == max( times[study==i])], pch=15)
}
Run the code above in your browser using DataLab