Survival analysis: up from Kaplan–Meier–Greenwood

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[-Sigma(j)(d(j)/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 {lambda(j)}, thus taking the SE of Sigma(j)(d(j)/T(j))t(j) to be [Sigma(j)(d(j)/T(j))2](1/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=Sigma(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 lambda, the corresponding survival-rate limit being exp[-Sigma(j)(d(j)*/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) >or= 1.