Modeling Fetal Death and Malformation in Developmental Toxicity Studies

We review approaches to dose-response modeling and risk assessment for binary data from developmental toxicity studies. In particular, we focus on jointly modeling fetal death and malformation and use a continuation ratio formulation of the multinomial distribution to provide a model for risk. Generalized estimating equations are used to account for clustering of animals within litters. The fitted model is then used to calculate doses corresponding to a specified level of excess risk. Two methods of arriving at a lower confidence limit or Benchmark dose are illustrated and compared. We also discuss models based on single binary end points and compare our approach to a binary analysis of whether or not the animal was 'affected' (either dead or malformed). The models are illustrated using data from four developmental toxicity studies in EG, DEHP, TGDM, and DYME conducted through the National Toxicology Program.