The Two-Dimensional Monte Carlo: A New Methodologic Paradigm for Dose Reconstruction for Epidemiological Studies

Abstract

Retrospective dose estimation, particularly dose reconstruction that supports epidemiological investigations of health risk, relies on various strategies that include models of physical processes and exposure conditions with detail ranging from simple to complex. Quantification of dose uncertainty is an essential component of assessments for health risk studies since, as is well understood, it is impossible to retrospectively determine the true dose for each person. To address uncertainty in dose estimation, numerical simulation tools have become commonplace and there is now an increased understanding about the needs and what is required for models used to estimate cohort doses (in the absence of direct measurement) to evaluate dose response. It now appears that for dose-response algorithms to derive the best, unbiased estimate of health risk, we need to understand the type, magnitude and interrelationships of the uncertainties of model assumptions, parameters and input data used in the associated dose estimation models. Heretofore, uncertainty analysis of dose estimates did not always properly distinguish between categories of errors, e.g., uncertainty that is specific to each subject (i.e., unshared error), and uncertainty of doses from a lack of understanding and knowledge about parameter values that are shared to varying degrees by numbers of subsets of the cohort. While mathematical propagation of errors by Monte Carlo simulation methods has been used for years to estimate the uncertainty of an individual subject’s dose, it was almost always conducted without consideration of dependencies between subjects. In retrospect, these types of simple analyses are not suitable for studies with complex dose models, particularly when important input data are missing or otherwise not available. The dose estimation strategy presented here is a simulation method that corrects the previous deficiencies of analytical or simple Monte Carlo error propagation methods and is termed, due to its capability to maintain separation between shared and unshared errors, the two-dimensional Monte Carlo (2DMC) procedure. Simply put, the 2DMC method simulates alternative, possibly true, sets (or vectors) of doses for an entire cohort rather than a single set that emerges when each individual’s dose is estimated independently from other subjects. Moreover, estimated doses within each simulated vector maintain proper inter-relationships such that the estimated doses for members of a cohort subgroup that share common lifestyle attributes and sources of uncertainty are properly correlated. The 2DMC procedure simulates inter-individual variability of possibly true doses within each dose vector and captures the influence of uncertainty in the values of dosimetric parameters across multiple realizations of possibly true vectors of cohort doses. The primary characteristic of the 2DMC approach, as well as its strength, are defined by the proper separation between uncertainties shared by members of the entire cohort or members of defined cohort su...

Thyroid Dose Estimates for a Cohort of Belarusian Children Exposed to131I from the Chernobyl Accident: Assessment of Uncertainties

“…This procedure is similar to and generally consistent with the 2-dimensional Monte Carlo method (14). For a specific dose realization some of the model parameter values were in common among members of subgroups, i.e., shared among subjects of those groups, implying that any error made on this parameter was shared by all subjects to whom it applied.…”

Concerted Uranium Research in Europe (CURE): toward a collaborative project integrating dosimetry, epidemiology and radiobiology to study the effects of occupational uranium exposure

“…The assessment and treatment of uncertainties is a complex area and approaches to this are still developing as knowledge improves. For instance, new approaches have been proposed since the end of CURE to take into account both shared and unshared sources of uncertainties on dose estimates in epidemiological analyses (Land et al 2015, Simon et al 2015, Stram et al 2015, Kwon et al 2016. It was recognised from the start of CURE that, due to the limited time and resources available within the project it would not be possible to develop the definitive final methodology for the treatment of uncertainties as part of future analyses, but also that this should be a priority for future research.…”

Accounting for Shared and Unshared Dosimetric Uncertainties in the Dose Response for Ultrasound-Detected Thyroid Nodules after Exposure to Radioactive Fallout

“…The current analysis utilizes a 2DMC approach that explicitly takes into account substantial and complex uncertainties in the individual parameters of the dose estimation model, including explicit separation of shared and unshared components (30, 31). …”

“…Following this characterization, the state of knowledge about each parameter used for dose estimation is quantitatively described by a probability density function (PDF) that characterizes random inter-individual variability and uncertainty in the center and spread of the PDF, as well as how the PDF is shared among relevant subgroups of the study population. The calculations to produce multiple realizations of simulated individual doses follow a detailed Monte Carlo sampling plan that maintains separation of shared from unshared sources of errors (31). In this study, we have used the 2DMC method to generate 5,000 realizations of dose-vector pairs, each pair representing a realization of possibly true cumulative external and internal doses for each of the 2,376 individuals.…”