plot.coxph_mpl function

Plot a coxph_mpl Object

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). utf8

## 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.

See Also

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

Author(s)

Dominique-Laurent Couturier, Maurizio Manuguerra

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)