coxrobust (version 1.0)

plot.coxr: Plot Diagnostics for a coxr Object

Description

Graphical tool which in a series of 5 graphs let us compare how well data are explained by the estimated proportional hazards model with non-robust (black color) and robust method (green color). The first graph gives standardized difference of two estimated survival functions; one via the Cox model and the other via Kaplan Meier estimator. The following four graphs show the same differences for four strata, defined by the quartiles of the estimated linear predictor. Comparison of estimation results along with analysis of the graphs leads frequently to a very detailed information about the model fit (see examples).

Usage

# S3 method for coxr
plot(x, caption = c("Full data set", "First quartile",
    "Second quartile", "Third quartile", "Fourth quartile"), main = NULL,
    xlab = "log time", ylab = "standardized survival differences", …,
    color = TRUE)

Arguments

x

coxr object, typically result of coxr.

caption

captions to appear above the plots.

xlab

title for the x axis.

ylab

title for the y axis.

main

overall title for the plot.

…

other parameters to be passed through to plotting functions.

color

if FALSE grayscale mode is used.

See Also

coxr

Examples

# NOT RUN {
#use the lung cancer data at Mayo Clinic 
#to compare results of non-robust and robust estimation
result <- coxr(Surv(time, status) ~ age + sex + ph.karno + meal.cal + wt.loss, data = lung)
plot(result, main = "Mayo Clinic Lung Cancer Data")
# }