pehaz {muhaz}

Estimates piecewise exponential hazard function from right-censored data.

Description

Divides the time domain into bins of equal width, and then estimates the hazard in each bin as the number of events in that bin divided by the total follow-up time in that bin.

Usage

pehaz(times, delta=NA, width=NA, min.time=0, max.time=NA)

Arguments

times

Vector of survival times. Does not need to be sorted.

delta

Vector indicating censoring 0 - censored (alive) 1 - uncensored (dead) If status is missing, all observations are assumed uncensored.

width

Bin width. Default value is that recommended by Mueller, width = (max.time-min.time) / (8*(nu)^0.2) where nu is the number of uncensored observations.

min.time

Left bound of the time domain used in the analysis. If missing, min.time is considered 0.

max.time

Right bound of the time domain used in the analysis. If missing, max.time is considered max(times).

Value

Returns an object of class 'pehaz', containing input and output values. Methods working on such an object are: plot, lines, print. For a detailed description of its components, see object.pehaz.

See Also

pehaz.object

Examples

data(cancer, package="survival")
attach(ovarian)
fit <- pehaz(futime, fustat)
plot(fit)