Bayesian methods for a three-state model for rodent carcinogenicity studies

The objective of a chronic rodent bioassay is to assess the impact of a chemical compound on the development of tumors. However, most tumor types are not observable prior to necropsy, making direct estimation of the tumor incidence rate problematic. In such cases, estimation can proceed only if the study incorporates multiple interim sacrifices or we make use of simplified parametric or nonparametric models. In addition, it is widely accepted that other factors, such as weight, can be related to both dose level and tumor onset, confounding the association of interest. However, there is not typically enough information in the current study to assess such effects. The addition of historical data can help alleviate this problem. In this article, we propose a novel Bayesian semiparametric model for the analysis of data from rodent carcinogenicity studies. We develop informative prior distributions for covariate effects through the use of historical control data and outline a Gibbs sampling scheme. We implement the model by analyzing data from a National Toxicology Program chronic rodent bioassay.