[R] survfit
Terry Therneau
therneau at mayo.edu
Wed Dec 15 14:21:49 CET 2010
> testpred
Call:
coxph(formula = Surv(time, status) ~ x, method = object$method)
coef exp(coef) se(coef) z p
karno -0.00756 0.992 0.0364 -0.208 0.84
Likelihood ratio test=0.89 on 1 df, p=0.344 n= 7
> newd1 <- list(testall$x[1,])
> survtest <- survfit(testpred,newdata=newd1,individual=T)
--------------
Do not use the individual=T argument. It only applies in the case of
time-dependent covariates. (In my latest version of the code this would
have been detected and gives an error message.)
To limit the printout to 30 and 80 days use summary(survest,
times=c(30,80))
Terry Therneau
More information about the R-help
mailing list