Last chance! 50% off unlimited learning
Sale ends in
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.
pehaz(times, delta=NA, width=NA, min.time=0, max.time=NA)
Vector of survival times. Does not need to be sorted.
Vector indicating censoring 0 - censored (alive) 1 - uncensored (dead) If status is missing, all observations are assumed uncensored.
Bin width. Default value is that recommended by
Mueller,
Left bound of the time domain used in the analysis. If missing, min.time is considered 0.
Right bound of the time domain used in the analysis. If missing, max.time is considered max(times).
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
.
# NOT RUN {
data(cancer, package="survival")
attach(ovarian)
fit <- pehaz(futime, fustat)
plot(fit)
# }
Run the code above in your browser using DataLab