Model averaging quantiles from data censored by a limit of detection

Abstract

In chemical risk assessment, it is important to determine the quantiles of the distribution of concentration data. The selection of an appropriate distribution and the estimation of particular quantiles of interest are largely hindered by the omnipresence of observations below the limit of detection, leading to left-censored data. The log-normal distribution is a common choice, but this distribution is not the only possibility and alternatives should be considered as well. Here, we focus on several distributions that are related to the log-normal distribution or that are seminonparametric extensions of the log-normal distribution. Whereas previous work focused on the estimation of the cumulative distribution function, our interest here goes to the estimation of quantiles, particularly in the left tail of the distribution where most of the left-censored data are located. Two different model averaged quantile estimators are defined and compared for different families of candidate models. The models and methods of selection and averaging are further investigated through simulations and illustrated on data of cadmium concentration in food products. The approach is extended to include covariates and to deal with uncertainty about the values of the limit of detection. These extensions are illustrated with (134) cesium measurements from Fukushima Prefecture, Japan. We can conclude that averaged models do achieve good performance characteristics in case no useful prior knowledge about the true distribution is available; that there is no structural difference in the performance of the direct and indirect method; and that, not surprisingly, only the true or closely approximating model can deal with extremely high percentages of censoring.

Editorial: Special issue on models for continuous data with a spike at zero

“…A simulation study shows that the proposed test has good properties in terms of both Type I error rate and power. Nysen et al (2016) discuss how to derive estimates for quantiles of interest in a chemical risk assessment study, when the observed (say continuous) distribution suffers from the presence of a spike at zero. These authors discuss cases where this spike is associated to left-censored data, that is to data below the limit of detection.…”