Nonparametric maximum likelihood estimation for the multisample Wicksell corpuscle problem

Biometrika. 2016 Jun;103(2):273-286. doi: 10.1093/biomet/asw011. Epub 2016 May 6.

Abstract

We study nonparametric maximum likelihood estimation for the distribution of spherical radii using samples containing a mixture of one-dimensional, two-dimensional biased and three-dimensional unbiased observations. Since direct maximization of the likelihood function is intractable, we propose an expectation-maximization algorithm for implementing the estimator, which handles an indirect measurement problem and a sampling bias problem separately in the E- and M-steps, and circumvents the need to solve an Abel-type integral equation, which creates numerical instability in the one-sample problem. Extensions to ellipsoids are studied and connections to multiplicative censoring are discussed.

Keywords: Abel-type integral equation; Expectation-maximization algorithm; Indirect measurement; Particle size.