Abstract:Background: Estimating the health effects of multi-pollutant mixtures is of increasing interest in environmental epidemiology. Recently, a new approach for estimating the health effects of mixtures, Bayesian kernel machine regression (BKMR), has been developed. This method estimates the multivariable exposure-response function in a flexible and parsimonious way, conducts variable selection on the (potentially high-dimensional) vector of exposures, and allows for a grouped variable selection approach that can accommodate highly correlated exposures. However, the application of this novel method has been limited by a lack of available software, the need to derive interpretable output in a computationally efficient manner, and the inability to apply the method to non-continuous outcome variables. Methods: This paper addresses these limitations by (i) introducing an open-source software package in the R programming language, the bkmr R package, (ii) demonstrating methods for visualizing high-dimensional exposure-response functions, and for estimating scientifically relevant summaries, (iii) illustrating a probit regression implementation of BKMR for binary outcomes, and (iv) describing a fast version of BKMR that utilizes a Gaussian predictive process approach. All of the methods are illustrated using fully reproducible examples with the provided R code. Results: Applying the methods to a continuous outcome example illustrated the ability of the BKMR implementation to estimate the health effects of multi-pollutant mixtures in the context of a highly nonlinear, biologically-based dose-response function, and to estimate overall, single-exposure, and interactive health effects. The Gaussian predictive process method led to a substantial reduction in the runtime, without a major decrease in accuracy. In the setting of a larger number of exposures and a dichotomous outcome, the probit BKMR implementation was able to correctly identify the variables included in the exposure-response function and yielded interpretable quantities on the scale of a latent continuous outcome or on the scale of the outcome probability. Conclusions: This newly developed software, integrated suite of tools, and extended methodology makes BKMR accessible for use across a broad range of epidemiological applications in which multiple risk factors have complex effects on health.

Bayesian Method for Improving Logistic Regression Estimates under Group-Based Exposure Assessment with Additive Measurement Errors

Bootstrap method-based estimation on the confidence interval for additive interaction in cohort studies

Robust Bayesian Inference for Berkson and Classical Measurement Error Models

Bayesian estimation of the measurement of interactions in epidemiological studies

A scalable two-stage Bayesian approach accounting for exposure measurement error in environmental epidemiology

A Bayesian Adjustment for Multiplicative Measurement Errors for a Calibration Problem with Application to a Stem Cell Study

A generalisation of the method of regression calibration and comparison with the Bayesian 2-dimensional Monte Carlo method

Adjusting for misclassification of an exposure in an individual participant data meta-analysis

A Bayesian approach to the g-formula

Bayesian multiple index models for environmental mixtures

Estimation and Inference for Case-Control Studies with Multiple Non-Gold Standard Exposure Assessments: with an Occupational Health Application.

Learning Models from Data with Measurement Error: Tackling Underreporting

Statistical software for analyzing the health effects of multiple concurrent exposures via Bayesian kernel machine regression

A Comparison of Statistical Methods for Studying Interactions of Chemical Mixtures

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

Empirical Bayes Estimates for a 2-Way Cross-Classified Additive Model

Bayesian outcome selection modelling

A latent functional approach for modeling the effects of multidimensional exposures on disease risk

Estimation of odds ratio from group testing data with misclassified exposure

Bayesian kernel machine regression for estimating the health effects of multi-pollutant mixtures

Estimating the health effects of environmental mixtures using Bayesian semiparametric regression and sparsity inducing priors