A Bayesian nonparametric analysis for zero‐inflated multivariate count data with application to microbiome study. Issue 4 (23rd June 2021)

• Record Type: Journal Article
• Authors: Shuler, Kurtis; Verbanic, Samuel; Chen, Irene A.; Lee, Juhee
• Is Part Of: Journal of the Royal Statistical Society. Volume 70:Issue 4(2021)
• Page Start: 961
• Page End: 979
• Publication Date: 2021-06-23
• Subjects: Bayesian nonparametrics -- dependent Dirichlet process -- high‐throughput sequencing -- microbiome -- multivariate count -- normalization -- operational taxonomic unit -- zero inflation; Statistics -- Periodicals; 519.5
• DOI: 10.1111/rssc.12493
• Languages: English
• ISSNs: 0035-9254

Abstract

Abstract: High‐throughput sequencing technology has enabled researchers to profile microbial communities from a variety of environments, but analysis of multivariate taxon count data remains challenging. We develop a Bayesian nonparametric (BNP) regression model with zero inflation to analyse multivariate count data from microbiome studies. A BNP approach flexibly models microbial associations with covariates, such as environmental factors and clinical characteristics. The model produces estimates for probability distributions which relate microbial diversity and differential abundance to covariates, and facilitates community comparisons beyond those provided by simple statistical tests. We compare the model to simpler models and popular alternatives in simulation studies, showing, in addition to these additional community‐level insights, it yields superior parameter estimates and model fit in various settings. The model's utility is demonstrated by applying it to a chronic wound microbiome data set and a Human Microbiome Project data set, where it is used to compare microbial communities present in different environments.