Implements product-limit (Kaplan-Meier) method for time-to-event data with censoring.
numeric vector with individual observation times
integer vector indicating status at the end of the observation time: 1 = event, 0 = censored
This function implements the product-limit method for estimating survival probability for time-to-event data with censoring:
tj are event times (i.e., times at which one or more events
nj are the number at risk at time
tj (by convention,
subjects censored at time
tj are considered at-risk and included in
dj are the number of events at time
A primary purpose of this function was to demonstrate the use of
available R functions to implement a simple statistical method. For example,
dimnames. Studying this function carefully helps one understand
and appreciate the utility of R functions to implement simple methods.
For serious survival analysis load the
survival package. The
survfit function in this package implements the product-limit method
and much more. See examples.
Returns an individual-level data frame
Visit http://medepi.com for the latest.
Selvin S. Statistical Analysis of Epidemiologic Data (Monographs in Epidemiology and Biostatistics, V. 35). Oxford University Press; 3rd edition (May 1, 2004)
1 2 3 4