Graphical comparing of cumulative hazards
Comparison of the cumulative hazards functions for a semi-parametric and parametric models.
plotHaz(
sp,
pp,
interval,
main = NULL,
xlab = "Time",
ylab = "Cum. hazards",
leglab,
col = c("blue", "red"),
lty = 1:2,
ylim,
log = "",
printLegend = TRUE
)sp |
An object of type "coxreg" or "phreg", typically output from
|
pp |
An object of type "coxreg" or "phreg", typically output from
|
interval |
Time interval for the plot, if missing, calculated from |
main |
Header for the plot. Default is distribution and "cumulative hazard function" |
xlab |
Label on x axis (default "Time") |
ylab |
Label on y axis (default "Cum. Hazards") |
leglab |
Labels in legend. |
col |
Line colors. should be |
lty |
line types. |
ylim |
Y limits for the plot. |
log |
Argument sent to |
printLegend |
Should a legend be printed? Default is |
For the moment only a graphical comparison. The arguments sp and
pp may be swapped.
No return value.
Göran Broström
check.dist, coxreg and phreg.
data(mort) op <- par(mfrow = c(1, 2)) fit.cr <- coxreg(Surv(enter, exit, event) ~ ses, data = mort) fit.w <- phreg(Surv(enter, exit, event) ~ ses, data = mort) fit.g <- phreg(Surv(enter, exit, event) ~ ses, data = mort, dist = "gompertz") plotHaz(fit.cr, fit.w, interval = c(0, 20), main = "Weibull") plotHaz(fit.cr, fit.g, main = "Gompertz") par(op)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.