check.dist: Graphical goodness-of-fit test

In goranbrostrom/eha: Event History Analysis

Graphical goodness-of-fit test

Description

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

Usage

check.dist(sp, pp, main = NULL, col = 1:2, lty = 1:2, printLegend = TRUE)

Arguments

sp	An object of type "coxreg", typically output from coxreg
pp	An object of type "phreg", typically output from phreg
main	Header for the plot. Default is distribution and "cumulative hazard function"
col	Line colors. should be NULL (black lines) or of length 2
lty	line types.
printLegend	Should a legend be printed? Default is TRUE.

Details

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

Value

No return value.

Author(s)

Göran Broström

See Also

coxreg and phreg.

Examples

data(mort)
oldpar <- par(mfrow = c(2, 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")
fit.ev <- phreg(Surv(enter, exit, event) ~ ses, data = mort,
dist = "ev")
check.dist(fit.cr, fit.w)
check.dist(fit.cr, fit.g)
check.dist(fit.cr, fit.ev)
par(oldpar)

goranbrostrom/eha documentation built on June 19, 2024, 1:39 a.m.