Open Access
March 2019 Sequential Dirichlet process mixtures of multivariate skew $t$-distributions for model-based clustering of flow cytometry data
Boris P. Hejblum, Chariff Alkhassim, Raphael Gottardo, François Caron, Rodolphe Thiébaut
Ann. Appl. Stat. 13(1): 638-660 (March 2019). DOI: 10.1214/18-AOAS1209


Flow cytometry is a high-throughput technology used to quantify multiple surface and intracellular markers at the level of a single cell. This enables us to identify cell subtypes and to determine their relative proportions. Improvements of this technology allow us to describe millions of individual cells from a blood sample using multiple markers. This results in high-dimensional datasets, whose manual analysis is highly time-consuming and poorly reproducible. While several methods have been developed to perform automatic recognition of cell populations most of them treat and analyze each sample independently. However, in practice individual samples are rarely independent, especially in longitudinal studies. Here we analyze new longitudinal flow-cytometry data from the DALIA-1 trial, which evaluates a therapeutic vaccine against HIV, by proposing a new Bayesian nonparametric approach with Dirichlet process mixture (DPM) of multivariate skew $t$-distributions to perform model based clustering of flow-cytometry data. DPM models directly estimate the number of cell populations from the data, avoiding model selection issues, and skew $t$-distributions provides robustness to outliers and nonelliptical shape of cell populations. To accommodate repeated measurements, we propose a sequential strategy relying on a parametric approximation of the posterior. We illustrate the good performance of our method on simulated data and on an experimental benchmark dataset. This sequential strategy outperforms all other methods evaluated on the benchmark dataset and leads to improved performance on the DALIA-1 data.


Download Citation

Boris P. Hejblum. Chariff Alkhassim. Raphael Gottardo. François Caron. Rodolphe Thiébaut. "Sequential Dirichlet process mixtures of multivariate skew $t$-distributions for model-based clustering of flow cytometry data." Ann. Appl. Stat. 13 (1) 638 - 660, March 2019.


Received: 1 July 2017; Revised: 1 July 2018; Published: March 2019
First available in Project Euclid: 10 April 2019

zbMATH: 07057442
MathSciNet: MR3937443
Digital Object Identifier: 10.1214/18-AOAS1209

Keywords: Automatic gating , Bayesian nonparametrics , Dirichlet process , flow cytometry , HIV , mixture model , skew $t$-distribution

Rights: Copyright © 2019 Institute of Mathematical Statistics

Vol.13 • No. 1 • March 2019
Back to Top