Bayesian benchmark dose analysis

An important objective in environmental risk assessment is the estimation of minimum exposure levels called benchmark doses (BMDs), which induce a pre-specified benchmark response in a target population. Established inferential approaches for BMD analysis typically involve one-sided, frequentist confidence limits, leading in practice to what are called benchmark dose lower Limits (BMDLs). Appeal to Bayesian modeling and credible limits for building BMDLs is far less developed, however. Indeed, for the few existing forms of Bayesian BMDs, informative prior information is seldom incorporated. We develop reparameterized quantal-response models that explicitly describe the BMD as a target parameter. Our goal is to obtain an improved estimation and calculation archetype for the BMD and for the BMDL by quantifying prior beliefs representing parameter uncertainty in the statistical model. Implementation is facilitated via a Monte Carlo-based adaptive Metropolis algorithm to approximate the posterior distribution. An example from environmental carcinogenicity testing illustrates the methodology. Copyright © 2015 John Wiley & Sons, Ltd.

[1]  George Casella,et al.  A Short History of Markov Chain Monte Carlo: Subjective Recollections from Incomplete Data , 2008, 0808.2902.

[2]  Kenny S. Crump,et al.  Calculation of Benchmark Doses from Continuous Data , 1995 .

[3]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[4]  Michael A. West,et al.  BAYESIAN MODEL ASSESSMENT IN FACTOR ANALYSIS , 2004 .

[5]  Louise Ryan,et al.  A Bayesian hierarchical model for risk assessment of methylmercury , 2003 .

[6]  A. Gelman.  Prior distributions for variance parameters in hierarchical models (comment on article by Browne and Draper) , 2004 .

[7]  Christophe Andrieu,et al.  A tutorial on adaptive MCMC , 2008, Stat. Comput..

[8]  Toxicology and carcinogenesis studies of cumene (CAS No. 98-82-8) in F344/N rats and B6C3F1 mice (inhalation studies). , 2009, National Toxicology Program technical report series.

[9]  Lingling An,et al.  The impact of model uncertainty on benchmark dose estimation , 2012, Environmetrics.

[10]  Matthew W Wheeler,et al.  Benchmark Dose Estimation Incorporating Multiple Data Sources , 2009, Risk analysis : an official publication of the Society for Risk Analysis.

[11]  Jeffrey S Gift,et al.  Introduction to benchmark dose methods and U.S. EPA's benchmark dose software (BMDS) version 2.1.1. , 2011, Toxicology and applied pharmacology.

[12]  Mitchell J. Small,et al.  Statistical Evaluation of Toxicological Experimental Design for Bayesian Model Averaged Benchmark Dose Estimation with Dichotomous Data , 2012 .

[13]  Kan Shao,et al.  Potential Uncertainty Reduction in Model‐Averaged Benchmark Dose Estimates Informed by an Additional Dose Study , 2011, Risk analysis : an official publication of the Society for Risk Analysis.

[14]  Malay Ghosh,et al.  Objective Priors: An Introduction for Frequentists , 2011, 1108.2120.

[15]  Lingling An,et al.  Information‐theoretic model‐averaged benchmark dose analysis in environmental risk assessment , 2013, Environmetrics.

[16]  Rutger van Haasteren,et al.  Gibbs Sampling , 2010, Encyclopedia of Machine Learning.

[17]  Xiao-Li Meng,et al.  SIMULATING RATIOS OF NORMALIZING CONSTANTS VIA A SIMPLE IDENTITY: A THEORETICAL EXPLORATION , 1996 .

[18]  J. Borwein,et al.  Two-Point Step Size Gradient Methods , 1988 .

[19]  A. John Bailer,et al.  Comparing model averaging with other model selection strategies for benchmark dose estimation , 2009, Environmental and Ecological Statistics.

[20]  A. Brix.  Bayesian Data Analysis, 2nd edn , 2005 .

[21]  S Zacks,et al.  Cancer phase I clinical trials: efficient dose escalation with overdose control. , 1998, Statistics in medicine.

[22]  A John Bailer,et al.  Monotonic Bayesian Semiparametric Benchmark
 Dose Analysis , 2012, Risk analysis : an official publication of the Society for Risk Analysis.

[23]  K S Crump,et al.  A new method for determining allowable daily intakes. , 1984, Fundamental and applied toxicology : official journal of the Society of Toxicology.

[24]  Ralph L. Kodell.  Managing uncertainty in health risk assessment , 2005 .

[25]  Bradley P. Carlin,et al.  Markov Chain Monte Carlo conver-gence diagnostics: a comparative review , 1996 .

[26]  Cody L. Wilson,et al.  Bayesian derivation of an oral cancer slope factor distribution for 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK). , 2009, Regulatory toxicology and pharmacology : RTP.

[27]  G. Casella,et al.  A History of Markov Chain Monte Carlo–Subjective Recollections from Incomplete Data– , 2008 .

[28]  John Geweke,et al.  Evaluating the accuracy of sampling-based approaches to the calculation of posterior moments , 1991 .

[29]  Qijun Fang,et al.  Hierarchical Bayesian Benchmark Dose Analysis , 2014 .

[30]  Jefferson Fowles,et al.  A benchmark dose analysis for sodium monofluoroacetate (1080) using dichotomous toxicity data. , 2007, Regulatory toxicology and pharmacology : RTP.

[31]  Ravi Varadhan,et al.  BB: An R Package for Solving a Large System of Nonlinear Equations and for Optimizing a High-Dimensional Nonlinear Objective Function , 2009 .

[32]  Anindya Roy,et al.  Nonparametric Bayesian Methods for Benchmark Dose Estimation , 2013, Risk analysis : an official publication of the Society for Risk Analysis.

[33]  Leonhard Held.  Simultaneous Inference in Risk Assessment; a Bayesian Perspective , 2004 .

[34]  Walter W. Piegorsch,et al.  Confidence limits on one-stage model parameters in benchmark risk assessment , 2009, Environmental and Ecological Statistics.

[35]  Kan Shao,et al.  A comparison of three methods for integrating historical information for Bayesian model averaged benchmark dose estimation. , 2012, Environmental toxicology and pharmacology.

[36]  Joseph G. Ibrahim,et al.  Bayesian Model Averaging With Applications to Benchmark Dose Estimation for Arsenic in Drinking Water , 2006 .

[37]  Herbert A. Sturges,et al.  The Choice of a Class Interval , 1926 .