Print number of observations, number of events, the restricted mean survival and its standard error, and the median survival with confidence limits for the median.
# S3 method for survfit
print(x, scale=1, digits = max(options()$digits - 4,3),
print.rmean=getOption("survfit.print.rmean"),
rmean = getOption('survfit.rmean'),...)x, with the invisible flag set to prevent printing. (The default for all print functions in R is to return the object passed to them; print.survfit complies with this pattern. If you want to capture these printed results for further processing, see the
table component of summary.survfit.)
the result of a call to the survfit function.
a numeric value to rescale the survival time, e.g.,
if the input data to survfit were in days,
scale=365 would scale the printout to years.
Number of digits to print
Options for computation and display of the restricted mean.
for future results
The number of observations, the number of events, the median survival
with its confidence interval, and optionally the
restricted mean survival (rmean) and its standard error, are printed.
If there are multiple curves, there is one line of output for each.
The mean and its variance are based on a truncated estimator. That is, if the
last observation(s) is not a death, then the survival curve estimate does not
go to zero and the mean is undefined.
There are four possible approaches to resolve this, which are selected by the
rmean option.
The first is to set the upper limit to a constant, e.g.,rmean=365.
In this case the reported mean would be the expected number of days, out
of the first 365, that would be experienced by each group. This is
useful if interest focuses on a fixed period.
Other options are "none" (no estimate), "common" and
"individual".
The "common" option uses the maximum time for all curves in the
object as a common upper limit for the auc calculation.
For the "individual"options the mean is computed as the area
under each curve,
over the range from 0 to the maximum observed time for that curve.
Since the end point is random, values for different curves are not
comparable and the printed standard errors are an underestimate as
they do not take into account this random variation. This option is
provided mainly for backwards compatability, as this estimate was the
default (only) one in earlier releases of the code.
Note that SAS (as of version 9.3) uses the integral up to the last
event time of each individual curve; we consider this the worst
of the choices and do not provide an option for that calculation.
The median and its confidence interval are defined by drawing a horizontal line at 0.5 on the plot of the survival curve and its confidence bands. If that line does not intersect the curve, then the median is undefined. The intersection of the line with the lower CI band defines the lower limit for the median's interval, and similarly for the upper band. If any of the intersections is not a point then we use the center of the intersection interval, e.g., if the survival curve were exactly equal to 0.5 over an interval. When data is uncensored this agrees with the usual definition of a median.
Miller, Rupert G., Jr. (1981). Survival Analysis. New York:Wiley, p 71.
summary.survfit, quantile.survfit