Is it possible to get survfit to produce the survival line for a single strata like preddow <- survfit(modall,newdata=newdat,se.fit=F,strata=2) # the strata argument is being ignored in the call above Or even get a more economical/faster calculation of the hazard directly from the coxph object?? Stephen [[alternative HTML version deleted]]