plot_km: Plot Survival Curve

In PWEXP: Piecewise Exponential Distribution Prediction Model

Plot Survival Curve

Description

Plot KM curve with right censoring data or the survival curve of a fitted piecewise exponential model.

Usage

## Default S3 method:
plot_survival(time, event, add=FALSE, conf.int=FALSE, mark.time=TRUE,
       lwd=2, xlab='Follow-up time', ylab='Survival function', ...)
## S3 method for class 'pwexp.fit'
plot_survival(time, add=TRUE, show_breakpoint=TRUE,
       breakpoint_par=NULL, ...)
## S3 method for class 'boot.pwexp.fit'
plot_survival(time, add=TRUE, alpha=0.1, show_breakpoint=TRUE,
       breakpoint_par=NULL, show_CI=TRUE, CI_par=NULL, ...)

Arguments

time	observed time from randomization or a pwexp.fit/ boot.pwexp.fit object.
event	the status indicator, normally 0=censor, 1=event. Other choices are TRUE/FALSE (TRUE = event).
add	logical; if TRUE add lines to current plot.
show_breakpoint	logical; if TRUE add vertial dashed lines to indicate breakpoints.
breakpoint_par	a list of parameters to control the apperance of vertical lines of breakpoionts. The values pass to abline.
alpha	the significance level of the confidence interval.
show_CI	logical; if TRUE add confidence interval of the estimated curve. For KM esitmator, use conf.int=TRUE to show CI band.
CI_par	a list of parameters to control the apperance of lines of confidence intervals. The values pass to lines.
conf.int	determines whether pointwise confidence intervals will be plotted. Passed over to plot.survfit.
mark.time	controls the labeling of the curves. Passed over to plot.survfit.
lwd	line width of the KM curve.
xlab	x label.
ylab	y label.
...	other arguments are passed over to plot.survfit (default method) or plot (for class pwexp.fit).

Details

For the default method, this a wrapper of plot.survfit function to plot right censoring data.

For class pwexp.fit, parameters in ... are passed to plot function to control the appearance of the survival curve; parameters in breakpoint_par are passed to abline function to control the appearance of vertical lines of breakpoints. See examples for usage.

For class boot.pwexp.fit, parameters in ... are passed to plot function to control the appearance of the survival curve; parameters in breakpoint_par are passed to abline function to control the appearance of vertical lines of breakpoints; parameters in CI_par are passed to lines function to control the appearance of confidence intervals. See examples for usage.

Value

No return value.

Author(s)

Tianchen Xu [e-mail]

See Also

plot_event

Examples

event_dist <- function(n)rpwexp(n, rate = c(0.1, 0.01, 0.2), breakpoint =  c(5,14))
dat <- simdata(rand_rate = 20, drop_rate = 0.03,  total_sample = 1000,
               advanced_dist = list(event_dist=event_dist),
               add_column = c('censor_reason','event','followT','followT_abs'))

plot_survival(dat$followT, dat$event, xlim=c(0,40))

fit_res <- pwexp.fit(dat$followT, dat$event, nbreak = 2)
plot_survival(fit_res, col='red', lwd=3, breakpoint_par = list(col='grey', lwd=2.5))

PWEXP documentation built on May 29, 2024, 2:19 a.m.