Data with some values below a limit of detection (LOD) can be analyzed using methods of survival analysis for left-censored data. The reverse Kaplan-Meier (KM) estimator provides an effective method for estimating the distribution function and thus population percentiles for such data. Although developed in the 1970s and strongly advocated since then, it remains rarely used, partly due to limited software availability.
In this paper, the reverse KM estimator is described and is illustrated using serum dioxin data from the University of Michigan Dioxin Exposure Study (UMDES) and the National Health and Nutrition Examination Survey (NHANES). Percentile estimates for left-censored data using the reverse KM estimator are compared with replacing values below the LOD with the LOD/2 or LOD/√2.
When some LODs are in the upper range of the complete values, and/or the percent censored is high, the different methods can yield quite different percentile estimates. The reverse KM estimator, which is the nonparametric maximum likelihood estimator, is the preferred method. Software options are discussed: The reverse KM can be calculated using software for the KM estimator. The JMP and SAS (SAS Institute, Cary, NC) and Minitab (Minitab, Inc, State College, PA), software packages calculate the reverse KM directly using their Turnbull estimator routines.
The reverse KM estimator is recommended for estimation of the distribution function and population percentiles in preference to commonly used methods such as substituting LOD/2 or LOD/√2 for values below the LOD, assuming a known parametric distribution, or using imputation to replace the left-censored values.