lifetab {biostat3} | R Documentation |
Create cohort life table
Description
Create cohort life table.
Usage
lifetab(tis, ninit, nlost, nevent)
Arguments
tis |
a vector of end points of time intervals, whose length is 1 greater than nlost and nevent. |
ninit |
the number of subjects initially entering the study. |
nlost |
a vector of the number of individuals lost follow or withdrawn alive for whatever reason. |
nevent |
a vector of the number of individuals who experienced the event |
Details
This is a minor update of the lifetab
function from the
KMsurv package, where the start and stop times of the intervals
are now included in the return value.
Value
A data.frame with the following columns:
tstart |
interval start time. |
tstop |
interval end time. |
nsubs |
the number of subject entering the intervals who have not experienced the event. |
nlost |
the number of individuals lost follow or withdrawn alive for whatever reason. |
nrisk |
the estimated number of individuals at risk of experiencing the event. |
nevent |
the number of individuals who experienced the event. |
surv |
the estimated survival function at the start of the intervals. |
pdf |
the estimated probability density function at the midpoint of the intervals. |
hazard |
the estimated hazard rate at the midpoint of the intervals. |
se.surv |
the estimated standard deviation of survival at the beginning of the intervals. |
se.pdf |
the estimated standard deviation of the prbability density function at the midpoint of the intervals. |
se.hazard |
the estimated standard deviation of the hazard function at the midpoint of the intervals |
The row.names are the intervals.
Author(s)
Jun Yan jyan@stat.uconn.edu
Examples
tis <- c(0, 2, 3, 5, 7, 11, 17, 25, 37, 53, NA)
nsubs <- c(927, 848, 774, 649, 565, 449, 296, 186, 112, 27)
nlost <- c(2, 3, 6, 9, 7, 5, 3, rep(0, 3))
nevent <- c(77, 71, 119, 75, 109, 148, 107, 74, 85, 27)
lifetab(tis, nsubs[1], nlost, nevent)