plot.libsl: Calibration Plot

Description

A calibration plot of an object of the class libsl (library of survival super learner).

Usage

# S3 method for libsl
plot(x, n.groups=5, pro.time=NULL,
newdata=NULL, ...)

Value

No return value for this S3 method.

Arguments

x: An object returned by a library of survival super learner.
n.groups: A numeric value with the number of groups by their class probabilities. The default is 5.
pro.time: The prognostic time at which the calibration plot of the survival probabilities.
newdata: An optional data frame containing the new sample for validation with covariate values, follow-up times, and event status. The default value is NULL, the calibration plot is performed from the same subjects of the training sample.
...: Additional arguments affecting the plot.

Details

The plot represents the observed survival and the related 95% confidence intervals, which are respectively estimated by the Kaplan and Meier estimator and the Greenwood formula, against the mean of the predictive values for individuals stratified into groups of the same size according to the percentiles. The identity line is usually included for reference.

Examples

Run this code

data("dataDIVAT2")

# The estimation of the model

formula<-Surv(times,failures) ~ age + hla + retransplant + ecd
data=dataDIVAT2[1:150,]
model <- LIB_COXall(formula, data=data)

# The calibration plot from the validation sample of 150 patients
plot(model, n.groups=5, pro.time=12, col=3,
     xlab="Predicted 12-year survival", ylab="Observed 12-year survival",
     newdata=dataDIVAT2[151:300,])