A comparison of methods for estimating the benchmark dose based on overdispersed data from developmental toxicity studies.

Developmental anomalies resulting from prenatal toxicity can be manifested in terms of both malformations among surviving offspring and prenatal death. Although these two endpoints have traditionally been analyzed separately in the assessment of risk, multivariate methods of risk characterization have recently been proposed. We examined this and other issues in developmental toxicity risk assessment by evaluating the accuracy and precision of estimates of the effective dose (ED05) and the benchmark dose (BMD05) using computer simulation. Our results indicated that different variance structures (Dirichlet-trinomial and generalized linear model) used to characterize overdispersion yielded comparable results when fitting joint dose response models based on generalized estimating equations. (The choice of variance structure in separate modeling was also not critical.) However, using the Rao-Scott transformation to eliminate overdispersion tended to produce estimates of the ED05 with reduced bias and mean squared error. Because joint modeling ensures that the ED05 for overall toxicity (based on both malformations and prenatal death) is always less than the ED05 for either malformations or prenatal death, joint modeling is preferred to separate modeling for risk assessment purposes.