Abstract:Statistical methods for identifying harmful chemicals in a correlated mixture often assume linearity in exposure-response relationships. Non-monotonic relationships are increasingly recognised (e.g., for endocrine-disrupting chemicals); however, the impact of non-monotonicity on exposure selection has not been evaluated. In a simulation study, we assessed the performance of Bayesian kernel machine regression (BKMR), Bayesian additive regression trees (BART), Bayesian structured additive regression with spike-slab priors (BSTARSS), and lasso penalised regression. We used data on exposure to 12 phthalates and phenols in pregnant women from the U.S. National Health and Nutrition Examination Survey to simulate realistic exposure data using a multivariate copula. We simulated datasets of size N = 250 and compared methods across 32 scenarios, varying by model size and sparsity, signal-to-noise ratio, correlation structure, and exposure-response relationship shapes. We compared methods in terms of their sensitivity, specificity, and estimation accuracy. In most scenarios, BKMR and BSTARSS achieved moderate to high specificity (0.56--0.91 and 0.57--0.96, respectively) and sensitivity (0.49--0.98 and 0.25--0.97, respectively). BART achieved high specificity ($\geq$ 0.96), but low to moderate sensitivity (0.13--0.66). Lasso was highly sensitive (0.75--0.99), except for symmetric inverse-U-shaped relationships ($\leq$ 0.2). Performance was affected by the signal-to-noise ratio, but not substantially by the correlation structure. Penalised regression methods that assume linearity, such as lasso, may not be suitable for studies of environmental chemicals hypothesised to have non-monotonic relationships with outcomes. Instead, BKMR and BSTARSS are attractive methods for flexibly estimating the shapes of exposure-response relationships and selecting among correlated exposures.

Bayesian outcome selection modelling

Bayesian variable selection for understanding mixtures in environmental exposures

Bayesian structural equation modeling for data from multiple cohorts

Model Selection of the Effect of Binary Exposures over the Life Course

Bayesian variable selection in linear regression models with instrumental variables

Bayesian analysis for pretest-posttest binary outcomes with adaptive significance levels

A Hierarchical Meta-Analysis for Settings Involving Multiple Outcomes across Multiple Cohorts

Assessing the Longitudinal Impact of Environmental Chemical Mixtures on Children's Neurodevelopment: A Bayesian Approach

Hair-Derived Exposome Exploration of Cardiometabolic Health: Piloting a Bayesian Multitrait Variable Selection Approach

Bayesian continuous-time hidden Markov models with covariate selection for intensive longitudinal data with measurement error.

Bayesian Quantile Regression with Subset Selection: A Posterior Summarization Perspective

Identifying Gene-environment interactions with robust marginal Bayesian variable selection

Estimating Perinatal Critical Windows of Susceptibility to Environmental Mixtures via Structured Bayesian Regression Tree Pairs

An estimating equations approach to fitting latent exposure models with longitudinal health outcomes

Performance of variable and function selection methods for estimating the non-linear health effects of correlated chemical mixtures: a simulation study

Multiple exposure distributed lag models with variable selection

Comparing Stochastic Optimization Methods for Variable Selection in Binary Outcome Prediction, With Application to Health Policy

A Bayesian multivariate factor analysis model for evaluating an intervention by using observational time series data on multiple outcomes

Bayesian Mendelian Randomization Analysis for Latent Exposures Leveraging GWAS Summary Statistics for Traits Co-Regulated by the Exposures

Bayesian adaptive and interpretable functional regression for exposure profiles

Using Bayesian item response theory for multicohort repeated measure design to estimate individual latent change scores.