plot.coxr {coxrobust} | R Documentation |
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 class '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
#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")
[Package
coxrobust version 1.0
Index]