Survival analysis: up from Kaplan–Meier–Greenwood
- 9 September 2008
- journal article
- Published by Springer Nature in European Journal of Epidemiology
- Vol. 23 (9) , 585-592
- https://doi.org/10.1007/s10654-008-9278-7
Abstract
In the type of survival analysis that now is routine, only the points of follow-up at which deaths from the cause at issue occurred make contributions to the Greenwood standard error (SE) of the survival rate’s Kaplan–Meier (KM) point estimate. An equivalent of this ‘KMG’ analysis draws from defined subintervals of the survival period being addressed. The data on each subinterval consist of the number of deaths from the cause at issue and the amount of population–time of follow-up, d j and T j , together with the duration of the interval, t j . The KM point estimate is replicated by \( \exp [-{\sum \nolimits_{j} ({{{d_{j} }/ \mathord{{\vphantom {{d_{j}}{T_{j}}}} \kern-0em} {T_{j} }}})t_{j}}], \) and the KMG interval estimate is replicated by treating the {d j } as a set of point estimates of Poisson parameters {λ j }, thus taking the SE of \( \sum \nolimits_{j}({{{d_{j} } \mathord{{\vphantom {{d_{j} } {T_{j} }}} \kern-0em} {T_{j}}}})t_{j}\) to be \([{\sum \nolimits _{j} d_{j} ({{{t_{j} } / \mathord{{\vphantom {{t_{j} } {T_{j} }}} \kern-0em} {T_{j} }}})^{2} }]^{{{1/ \mathord{{\vphantom {1 2}} \kern-0em} 2}}}.\) In both the KMG analysis and this equivalent of it, the SE used to derive the survival rate’s lower confidence limit needs to be augmented by a factor that accounts for the loss of information due to censorings subsequent to the last ‘failure’ in the survival period at issue. But, SE-based interval estimation of survival rate actually needs to be replaced by a first-principles counterpart of it. A suitable point of departure in this is first-principles asymptotic interval estimation of the Poisson parameter \( \lambda =\sum \nolimits _{j}{\lambda_{j}}, \) if not the exact counterpart of this. A confidence limit for the survival rate can then be based on suitable augmentation or contraction of the {d j } set to \( \{ d_{j}^{*} \} \) consistent with a given limit for λ, the corresponding survival-rate limit being \( \exp [-{\sum \nolimits _{j} ({{{d_{j}^{ *} } \mathord{/ {\vphantom {{d_{j}^{ *} } {T_{j} }}} \kern-0em} {T_{j}}}})t_{j}}]. \) Suitable augmentation is constituted by an identical addition to each \(d_{j}^{1/2},\) suitable contraction by an identical subtraction from each \(d_{j}^{1/2} \ge 1.\)
Keywords
This publication has 5 references indexed in Scilit:
- Survival of Patients with Stage I Lung Cancer Detected on CT ScreeningNew England Journal of Medicine, 2006
- A simple hybrid variance estimator for the Kaplan–Meier survival functionStatistics in Medicine, 2004
- Design and analysis of randomized clinical trials requiring prolonged observation of each patient. II. Analysis and examplesBritish Journal of Cancer, 1977
- CIGARETTE-SMOKING AND NONFATAL MYOCARDIAL INFARCTION: RATE RATIO IN RELATION TO AGE, SEX AND PREDISPOSING CONDITIONS1American Journal of Epidemiology, 1976
- Nonparametric Estimation from Incomplete ObservationsJournal of the American Statistical Association, 1958