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