eha (version 2.8.5)

plotHaz: Graphical comparing of cumulative hazards

Description

Comparison of the cumulative hazards functions for a semi-parametric and parametric models.

Usage

plotHaz(
  sp,
  pp,
  interval,
  main = NULL,
  xlab = "Time",
  ylab = "Cum. hazards",
  col = c("blue", "red"),
  lty = 1:2,
  ylim,
  printLegend = TRUE
)

Arguments

sp

An object of type "coxreg" or "phreg", typically output from coxreg or phreg.

pp

An object of type "coxreg" or "phreg", typically output from coxreg or phreg.

interval

Time interval for the plot, if missing, calculated from sp.

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")

col

Line colors. should be NULL (black lines) or of length 2

lty

line types.

ylim

Y limits for the plot.

printLegend

Should a legend be printed? Default is TRUE.

Value

No return value.

Details

For the moment only a graphical comparison. The arguments sp and pp may be swapped.

See Also

check.dist, coxreg and phreg.

Examples

# NOT RUN {
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)

# }

Run the code above in your browser using DataCamp Workspace