Cytometry inference through adaptive atomic deconvolution

Manon Costa (1) , Sébastien Gadat (2) , Pauline Gonnord (3) , Laurent Risser (1)

1 IMT - Institut de Mathématiques de Toulouse UMR5219
2 TSE-R - Toulouse School of Economics
3 CPTP - Centre de Physiopathologie Toulouse Purpan

In this paper we consider a statistical estimation problem known as atomic deconvolution. Introduced in reliability, this model has a direct application when considering biological data produced by flow cytometers. In these experiments, biologists measure the fluorescence emission of treated cells and compare them with their natural emission to study the presence of specific molecules on the cells' surface. They observe a signal which is composed of a noise (the natural fluorescence) plus some additional signal related to the quantity of molecule present on the surface if any. From a statistical point of view, we aim at inferring the percentage of cells expressing the selected molecule and the probability distribution function associated with its fluorescence emission. We propose here an adap-tive estimation procedure based on a previous deconvolution procedure introduced by [?, ?]. For both estimating the mixing parameter and the mixing density automatically, we use the Lepskii method based on the optimal choice of a bandwidth using a bias-variance decomposition. We then derive some concentration inequalities for our estimators and obtain the convergence rates, that are shown to be minimax optimal (up to some log terms) in Sobolev classes. Finally, we apply our algorithm on simulated and real biological data.

Mots clés: 62P10; Adaptive kernel estimators; Inverse problems MSC2010: Primary: 62G07 Secondary: 62G20; Atomic deconvolution; Mixture models; Inverse problems

Domaines: Statistiques [math.ST]

Manon Costa, Sébastien Gadat, Pauline Gonnord, Laurent Risser. Cytometry inference through adaptive atomic deconvolution. 2017. ⟨hal-01613723⟩