Frank E Harrell Jr schrieb: > Matthias Gondan wrote: >>> data(colon) >>> s = survfit(Surv(time, status) ~ rx, data=colon) >>> plot(s) >>> plot(s, col=1:3) >> >> By the way: Does anyone know a neat way to indicate the number of >> patients under risk in >> the curve? > That's one of many options in the survplot function in the Design > package, which uses the survival package.
Ah, thanks! Nice. data(colon) s = survfit(Surv(time, status) ~ rx, data=colon) survplot(s, n.risk=TRUE, conf='none', time.inc=500) Is there an option for drawing small vertical lines if a patient has been censored? Best wishes, Matthias ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.