Revisiting Cohen et al. 2015, Cohen et al. 2014 and Waalkes et al. 2014: a bayesian re-analysis of tumor incidences
Revisiting Cohen et al. 2015, Cohen et al. 2014 and Waalkes et al. 2014: a bayesian re-analysis of tumor incidences
复制标题
回顾科恩等人。
DOI:
10.1007/s00204-016-1749-0
复制
发表时间:
2016
影响因子:
6.1
通讯作者:
L. Burgoon
中科院分区:
文献类型:
--
作者:
Ingrid L. Druwe;L. Burgoon
The control animal tumor incidences reported by Tokar et al. in 2011 are not different from the study reported by the same group in the Waalkes et al. 2014 publication. Our full analysis can be found in Burgoon and Druwe (2015); however, we briefly discuss the approach and results here. To test the null hypothesis, we used a Bayesian approach. Initially, we examined whether or not a difference existed without using any prior knowledge of what the tumor incidence in Waalkes’ laboratory was; therefore, we used a flat prior distribution. We modeled the control tumor incidences from the Tokar 2011 and Waalkes studies as Bernoulli distributions. The posterior distributions were calculated and are shown in Burgoon and Druwe (2015). There is a difference of about 13% between the means of these two distributions. However, that fact alone does not mean that the incidences are from different distributions. In fact, what we could be observing is a case where the incidences for each study were taken from different sides of the same distribution. Thus, to test the hypothesis, we took samples from the posterior distributions and calculated the difference to obtain a distribution of the differences. If the incidences in both studies were from the same distribution, we would expect a difference of 0, or close to 0, to be a credible value. In order to accomplish this, we used an approach that sets a region of practical equivalence (ROPE) around the zero difference, and the 95% highest density interval (HDI) of the difference distribution. The ROPE demarcates a region around zero difference that is functionally equivalent to no difference. In general, if any part of the 95% HDI is within the ROPE, then we accept the null hypothesis that the control tumor incidences from the studies are the same. Else, we reject the null hypothesis. A complete explanation of our decision rules can be found in Burgoon and Druwe (2015).Over the past year, this journal has published a lively discussion in the form of letters to the editor between Waalkes et al. and Cohen et al. regarding an article published in this journal titled “Lung Tumors in Mice induced by ‘whole Life’inorganic arsenic exposure at human relevant doses” by Waalkes et al. in 2014. Cohen et al. raised a series of thoughtful questions with respect to the reproducibility of the control animal tumor incidences in the arsenic exposure studies published by Tokar et al.(2011) and Waalkes et al.(2014). In addition, Cohen et al. brought into question whether the development of lung tumors in the mice was related to the genetic background of the mice used in the study rather than arsenic exposure. Many of the questions raised by Cohen and colleagues centered around what they deemed to be uncertainty in the tumor incidences in the control animals used in both studies, and by extension, the quality of the studies performed by Tokar and Waalkes. If the assertions made by Cohen et al. proved true, risk assessors would be unable to use Waalkes et al.’s (2014) study in hazard and dose–response assessments of inorganic arsenic. Thus, we performed the analysis requested