On Tue, 2 Nov 2004, Min-Han Tan wrote:
Hi,
Sorry to trouble the list. I have a problem which I'm not sure how to resolve.
I have a Cox model with 1 independent variable with 2 categories (and thus 2 survival curves on plotting survfit)
How can I get an estimate of survival for each category at a particular time point, with standard error?
Here's an example
fit <- coxph( Surv(futime,fustat)~rx,data=ovarian)
summary( survfit( fit, newdata=list(rx=c(1,2))))
Call: survfit.coxph(object = fit, newdata = list(rx = c(1, 2)))
time n.risk n.event survival1 survival2 59 26 1 0.952 0.973 115 25 1 0.903 0.945 156 24 1 0.855 0.917 268 23 1 0.806 0.888 329 22 1 0.758 0.858 353 21 1 0.710 0.828 365 20 1 0.663 0.797 431 17 1 0.613 0.764 464 15 1 0.560 0.727 475 14 1 0.509 0.689 563 12 1 0.454 0.648 638 11 1 0.403 0.606
-thomas
______________________________________________ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
