Email-запись: Bayesian hierarchical negative binomial models for multivariable analyses with applications to human microbiome count data