Usage
survplot(fit, ...)
## S3 method for class 'rms':
survplot(fit, \dots, xlim,
ylim=if(loglog) c(-5, 1.5) else if
(what == "survival" & missing(fun)) c(0, 1),
xlab, ylab, time.inc,
what=c("survival","hazard"),
type=c("tsiatis","kaplan-meier"),
conf.type=c("log","log-log","plain","none"),
conf.int=FALSE, conf=c("bands","bars"),
add=FALSE, label.curves=TRUE,
abbrev.label=FALSE, levels.only=FALSE,
lty, lwd=par("lwd"),
col=1, col.fill=gray(seq(.95, .75, length=5)),
adj.subtitle=TRUE, loglog=FALSE, fun,
n.risk=FALSE, logt=FALSE, dots=FALSE, dotsize=.003,
grid=NULL, srt.n.risk=0, sep.n.risk=0.056, adj.n.risk=1,
y.n.risk, cex.n.risk=.6, pr=FALSE)
## S3 method for class 'survfit':
survplot(fit, xlim,
ylim, xlab, ylab, time.inc,
conf=c("bands","bars","diffbands","none"), add=FALSE,
label.curves=TRUE, abbrev.label=FALSE,
levels.only=FALSE, lty,lwd=par('lwd'),
col=1, col.fill=gray(seq(.95, .75, length=5)),
loglog=FALSE,fun,n.risk=FALSE,logt=FALSE,
dots=FALSE,dotsize=.003, grid=NULL,
srt.n.risk=0,sep.n.risk=.056,adj.n.risk=1,
y.n.risk,cex.n.risk=.6, pr=FALSE, ...)
survdiffplot(fit, order=1:2, fun=function(y) y,
xlim, ylim, xlab, ylab="Difference in Survival Probability",
time.inc, conf.int, conf=c("shaded", "bands","diffbands","none"),
add=FALSE, lty=1, lwd=par('lwd'), col=1,
n.risk=FALSE, grid=NULL,
srt.n.risk=0, adj.n.risk=1,
y.n.risk, cex.n.risk=.6)Arguments
fit
result of fit (cph, psm, survfit,
survest.psm). For survdiffplot, fit must be the
result of survfit.
...
list of factors with names used in model. For fits from survfit,
these arguments do not appear - all strata are plotted. Otherwise the
first factor listed is the factor used to determine different survival
curves. Any other factors are used
xlim
a vector of two numbers specifiying the x-axis range for follow-up time.
Default is (0,maxtime) where maxtime was the pretty()d version
of the maximum follow-up time
in any stratum, stored in fit$maxtime
ylim
y-axis limits. Default is c(0,1) for survival, and
c(-5,1.5) if loglog=TRUE. If fun or
loglog=TRUE are given and ylim is not, the limits will be
computed from the data. For
xlab
x-axis label. Default is units attribute of failure time
variable given to Surv.
ylab
y-axis label. Default is "Survival Probability" or
"log(-log Survival Probability)". If fun is given, the default
is "". For what="hazard", the default is "Hazard Function".
time.inc
time increment for labeling the x-axis and printing numbers at risk.
If not specified, the value
of time.inc stored with the model fit will be used.
type
specifies type of estimates, "tsiatis" (the default) or "kaplan-meier".
"tsiatis" here corresponds to the Breslow
estimator. This is ignored if survival estimates stored with surv=TRUE are
being used. Fo
conf.type
specifies the basis for confidence limits. This argument is
ignored for fits from survfit.
conf.int
Default is FALSE. Specify e.g. .95 to plot 0.95 confidence bands.
For fits from parametric survival models, or Cox models with
x=TRUE and y=TRUE specified to the fit, the exact
asymptotic formulas will
conf
"bars" for confidence bars at each time.inc time point. If the fit
was from cph(..., surv=TRUE), the time.inc used will be
that stored with the fit. Use conf="bands" (the default) for
bands
what
defaults to "survival" to plot survival estimates. Set to
"hazard" or an abbreviation to plot the hazard function (for
psm fits only). Confidence intervals are not available for
what="hazard".
add
set to TRUE to add curves to an existing plot.
label.curves
default is TRUE to use labcurve to label curves where they are farthest
apart. Set label.curves to a list to specify options to
labcurve, e.g., label.curves=list(method="arrow", cex=.
abbrev.label
set to TRUE to abbreviate() curve labels that are plotted
levels.only
set to TRUE to remove variablename= from the start of
curve labels.
lty
vector of line types to use for different factor levels. Default is
c(1,3,4,5,6,7,...).
lwd
vector of line widths to use for different factor levels. Default is
current par setting for lwd.
col
color for curve, default is 1. Specify a vector to assign different
colors to different curves.
col.fill
a vector of colors to used in filling confidence bands
adj.subtitle
set to FALSE to suppress plotting subtitle with levels of adjustment factors
not plotted. Defaults to TRUE.
This argument is ignored for survfit.
loglog
set to TRUE to plot log(-log Survival) instead of Survival
fun
specifies any function to translate estimates and confidence limits
before plotting.
logt
set to TRUE to plot log(t) instead of t on the x-axis
n.risk
set to TRUE to add number of subjects at risk for each curve, using the
surv.summary created by cph or using the failure times used in
fitting the model if y=TRUE was specified to the fit or if the fit
w
srt.n.risk
angle of rotation for leftmost number of subjects at risk (since this number
may run into the second or into the y-axis). Default is 0.
adj.n.risk
justification for leftmost number at risk. Default is 1 for right
justification.
Use 0 for left justification, .5 for centered.
sep.n.risk
multiple of upper y limit - lower y limit for separating lines of text
containing number of subjects at risk. Default is .056*(ylim[2]-ylim[1]).
y.n.risk
When n.risk=TRUE, the default is to place numbers of patients at risk above
the x-axis. You can specify a y-coordinate for the bottom line of the
numbers using y.n.risk.
cex.n.risk
character size for number of subjects at risk (when n.risk is TRUE)
dots
set to TRUE to plot a grid of dots. Will be plotted at every
time.inc (see cph) and at survival increments of .1 (if
d>.4), .05 (if .2 < d <= .4<="" code="">), or .025 (if d <= .2<="" code="">),
wher
dotsize
size of dots in inches
grid
defaults to NULL (not drawing grid lines). Set to TRUE to
plot gray(.8) grid lines, or specify any color.
pr
set to TRUE to print survival curve coordinates used in the plots
order
an integer vector of length two specifying the order of groups when
computing survival differences. The default of 1:2 indicates
that the second group is subtracted from the first. Specify
order=2:1 to instead subtract th