[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