plot.coxph_mpl {survivalMPL}

Plot a coxph_mpl Object

Description

Plot the bases used to estimate the baseline hazard parameters, as well as the estimate and confidence interval of the baseline hazard, cumulative baseline hazard and baseline survival functions (plots are selectable by which).

Usage

## S3 method for class 'coxph_mpl'
plot(x, se="M2QM2", ask=TRUE, which=1:4, upper.quantile=.95,...)

Arguments

x

an object inheriting from class coxph_mpl, representing a fitted Cox proportional hazard model.

se

an inference method (to build confidence intevals for the baseline hazard, cumulative baseline hazard and baseline survival functions). Possibilites are "H", "M2QM2" and "M2HM2". Refer to the Details Section of coxph_mpl. Default is se="M2QM2".

ask

logical. If TRUE, the user is asked to hit the enter keyboard before each plot. See par(ask=.). Default is ask=TRUE.

which

integer vector indicating the list of wished plots. If a subset of the plots is required, specify a subset of the numbers 1:4. By default, all plots are provided.

upper.quantile

quantile of the model response defining the upper limit of the x-axis of the plots of the baseline hazard, cumulative baseline hazard and baseline survival functions. Default is upper.quantile=.95.

...

other parameters to be passed through to plotting functions.

Details

In the first plot, the bases corresponding to zero (or close to zero) estimates appear in dashed line. An estimate is considered as a zero if it is smaller than min.Theta (See coxph_mpl.control).

Confidence intervals for the baseline hazard, cumulative baseline hazard and baseline survival functions are obtained using the delta method.

Author(s)

Dominique-Laurent Couturier, Maurizio Manuguerra

See Also

coxph_mpl, coxph_mpl.control, coxph_mpl.object and summary.coxph_mpl.

Examples

## Not run: 
data(lung)
fit_mpl <- coxph_mpl(Surv(time, status == 2) ~ age + sex + ph.karno + wt.loss, data = lung)

par(mfrow=c(2,2))
plot(fit_mpl, ask=FALSE, cex.main=.75)

## End(Not run)