Risk Assessment for Toxicity Experiments with Discrete and Continuous Outcomes: A Bayesian Nonparametric Approach

We present a Bayesian nonparametric modeling approach to inference and risk assessment for developmental toxicity studies. The primary objective of these studies is to determine the relationship between the level of exposure to a toxic chemical and the probability of a physiological or biochemical response. We consider a general data setting involving clustered categorical responses on the number of prenatal deaths, the number of live pups, and the number of live malformed pups from each laboratory animal, as well as continuous outcomes (e.g., body weight) on each of the live pups. We utilize mixture modeling to provide flexibility in the functional form of both the multivariate response distribution and the various dose–response curves of interest. The nonparametric model is built from a structured mixture kernel and a dose-dependent Dirichlet process prior for the mixing distribution. The modeling framework enables general inference for the implied dose–response relationships and for dose-dependent correlations between the different endpoints, features which provide practical advances relative to traditional parametric models for developmental toxicology. We use data from a toxicity experiment that investigated the toxic effects of an organic solvent (diethylene glycol dimethyl ether) to demonstrate the range of inferences obtained from the nonparametric mixture model, including comparison with a parametric hierarchical model. Supplementary materials accompanying this paper appear on-line.