\HeaderA{print.survfit}{Print a Short Summary of a Survival Curve}{print.survfit}
\keyword{survival}{print.survfit}
\begin{Description}\relax
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.
\end{Description}
\begin{Usage}
\begin{verbatim}
## S3 method for class 'survfit':
print(x, scale=1,digits = max(options()$digits - 4,
3), print.n=getOption("survfit.print.n"), show.rmean=getOption("survfit.print.mean"),...)
\end{verbatim}
\end{Usage}
\begin{Arguments}
\begin{ldescription}
\item[\code{x}] the result of a call to the \code{survfit} function.

\item[\code{print.n}] What to use for number of subjects (see below)
\item[\code{digits}] Number of digits to print
\item[\code{scale}] a numeric value to rescale the survival time, e.g.,
if the input data to survfit were in days,
\code{scale=365} would scale the printout to years.

\item[\code{show.rmean}] Show the restricted mean survival?
\item[\code{...}] other unused arguments
\end{ldescription}
\end{Arguments}
\begin{Details}\relax
The restricted mean (\code{rmean}) and its standard error \code{se(rmean)} are based on a truncated estimator.  If the
last observation(s) is not a death, then the survival curve estimate does not
go to zero and the mean survival time cannot be estimated. Instead, the
quantity reported is the mean of survival restricted to the time before
the last censoring. When the last censoring time is not random this
quantity is occasionally of interest.

Any randomness in the last censoring time is not taken into account in
computing the standard error of the restricted mean.  The restricted
mean is available mainly for compatibility with S, and is not shown by default.


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.
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 smallest point of intersection,
e.g., if the survival curve were exactly equal to 0.5 over an
interval. 

The "number of observations" is not well-defined for counting process
data. Previous versions of this code used the number at risk at the
first time point. This is misleading if many individuals enter late or
change strata. The original S code for the current version uses the
number of records, which is misleading when the counting process data
actually represent a fixed cohort with time-dependent covariates.

Four possibilities are provided, controlled by \code{print.n} or by
\code{options(survfit.print.n)}: \code{"none"} prints \code{NA},
\code{"records"} prints the number of records, \code{"start"} prints the
number at the first time point and \code{"max"} prints the maximum
number at risk. The initial default is \code{"start"}.
\end{Details}
\begin{Value}
\code{x}, with the invisible flag set.
\end{Value}
\begin{Section}{Side Effects}
The number of observations (see Details), the number of events,  the median survival
with its confidence interval, and optionally the
restricted mean survival (\code{rmean}) and its standard error, are printed.  If there are multiple curves, there is one line of output for each.
\end{Section}
\begin{SeeAlso}\relax
\code{\LinkA{summary.survfit}{summary.survfit}}, \code{\LinkA{survfit.object}{survfit.object}}, \code{\LinkA{survfit}{survfit}}
\end{SeeAlso}
\begin{Examples}
\begin{ExampleCode}
##effect of print.n and show.rmean

a<-coxph(Surv(start,stop,event)~age+strata(transplant),data=heart)
b<-survfit(a)
print(b,print.n="none")
print(b,print.n="records")
print(b,print.n="start")
print(b,print.n="max")
print(b,show.rmean=TRUE)
\end{ExampleCode}
\end{Examples}

