蒙特卡罗方法及其在辐射剂量计算中的应用

概述了蒙特卡罗方法的产生与发展,阐述了蒙特卡罗方法的基本特点,最后就蒙特卡罗方法在辐射剂量计算上的应用进行了讨论。     

Low dose risk estimation via simultaneous statistical inferences

Summary.  The paper develops and studies simultaneous confidence bounds that are useful for making low dose inferences in quantitative risk analysis. Application is intended for risk assessment studies where human, animal or ecological data are used to set safe low dose levels of a toxic agent, but where study information is limited to high dose levels of the agent. Methods are derived for estimating simultaneous, one-sided, upper confidence limits on risk for end points measured on a continuous scale. From the simultaneous confidence bounds, lower confidence limits on the dose that is associated with a particular risk (often referred to as a bench-mark dose ) are calculated. An important feature of the simultaneous construction is that any inferences that are based on inverting the simultaneous confidence bounds apply automatically to inverse bounds on the bench-mark dose.     

Benchmark Analysis: Shopping with Proper Confidence

We discuss the issue of using benchmark doses for quantifying (excess) risk associated with exposure to environmental hazards. The paradigm of low-dose risk estimation in dose-response modeling is used as the primary application scenario. Emphasis is placed on making simultaneous inferences on benchmark doses when data are in the form of proportions, although the concepts translate easily to other forms of outcome data.     

Predictive Bayesian Microbial Dose-Response Assessment Based on Suggested Self-Organization in Primary Illness Response: Cryptosporidium parvum

The probability of illness caused by very low doses of pathogens cannot generally be tested due to the numbers of subjects that would be needed, though such assessments of illness dose response are needed to evaluate drinking water standards. A predictive Bayesian dose-response assessment method was proposed previously to assess the unconditional probability of illness from available information and avoid the inconsistencies of confidence-based approaches. However, the method uses knowledge of the conditional dose-response form, and this form is not well established for the illness endpoint. A conditional parametric dose-response function for gastroenteric illness is proposed here based on simple numerical models of self-organized host-pathogen systems and probabilistic arguments. In the models, illnesses terminate when the host evolves by processes of natural selection to a self-organized critical value of wellness. A generalized beta-Poisson illness dose-response form emerges for the population as a whole. Use of this form is demonstrated in a predictive Bayesian dose-response assessment for cryptosporidiosis. Results suggest that a maximum allowable dose of 5.0 x 10(-7) oocysts/exposure (e.g., 2.5 x 10(-7) oocysts/L water) would correspond with the original goals of the U.S. Environmental Protection Agency Surface Water Treatment Rule, considering only primary illnesses resulting from Poisson-distributed pathogen counts. This estimate should be revised to account for non-Poisson distributions of Cryptosporidium parvum in drinking water and total response, considering secondary illness propagation in the population.     

Bootstrap Estimation of Benchmark Doses and Confidence Limits with Clustered Quantal Data

The benchmark dose (BMD) is an exposure level that would induce a small risk increase (BMR level) above the background. The BMD approach to deriving a reference dose for risk assessment of noncancer effects is advantageous in that the estimate of BMD is not restricted to experimental doses and utilizes most available dose-response information. To quantify statistical uncertainty of a BMD estimate, we often calculate and report its lower confidence limit (i.e., BMDL), and may even consider it as a more conservative alternative to BMD itself. Computation of BMDL may involve normal confidence limits to BMD in conjunction with the delta method. Therefore, factors, such as small sample size and nonlinearity in model parameters, can affect the performance of the delta method BMDL, and alternative methods are useful. In this article, we propose a bootstrap method to estimate BMDL utilizing a scheme that consists of a resampling of residuals after model fitting and a one-step formula for parameter estimation. We illustrate the method with clustered binary data from developmental toxicity experiments. Our analysis shows that with moderately elevated dose-response data, the distribution of BMD estimator tends to be left-skewed and bootstrap BMDL s are smaller than the delta method BMDL s on average, hence quantifying risk more conservatively. Statistically, the bootstrap BMDL quantifies the uncertainty of the true BMD more honestly than the delta method BMDL as its coverage probability is closer to the nominal level than that of delta method BMDL. We find that BMD and BMDL estimates are generally insensitive to model choices provided that the models fit the data comparably well near the region of BMD. Our analysis also suggests that, in the presence of a significant and moderately strong dose-response relationship, the developmental toxicity experiments under the standard protocol support dose-response assessment at 5% BMR for BMD and 95% confidence level for BMDL.     

Modeling for Risk Assessment of Neurotoxic Effects

The regulation of noncancer toxicants, including neurotoxicants, has usually been based upon a reference dose (allowable daily intake). A reference dose is obtained by dividing a no-observed-effect level by uncertainty (safety) factors to account for intraspecies and interspecies sensitivities to a chemical. It is assumed that the risk at the reference dose is negligible, but no attempt generally is made to estimate the risk at the reference dose. A procedure is outlined that provides estimates of risk as a function of dose. The first step is to establish a mathematical relationship between a biological effect and the dose of a chemical. Knowledge of biological mechanisms and/or pharmacokinetics can assist in the choice of plausible mathematical models. The mathematical model provides estimates of average responses as a function of dose. Secondly, estimates of risk require selection of a distribution of individual responses about the average response given by the mathematical model. In the case of a normal or lognormal distribution, only an estimate of the standard deviation is needed. The third step is to define an adverse level for a response so that the probability (risk) of exceeding that level can be estimated as a function of dose. Because a firm response level often cannot be established at which adverse biological effects occur, it may be necessary to at least establish an abnormal response level that only a small proportion of individuals would exceed in an unexposed group. That is, if a normal range of responses can be established, then the probability (risk) of abnormal responses can be estimated. In order to illustrate this process, measures of the neurotransmitter serotonin and its metabolite 5-hydroxyindoleacetic acid in specific areas of the brain of rats and monkeys are analyzed after exposure to the neurotoxicant methylene-dioxymethamphetamine. These risk estimates are compared with risk estimates from the quantal approach in which animals are classified as either abnormal or not depending upon abnormal serotonin levels.     

Summary of the Workshop on Issues in Risk Assessment: Quantitative Methods for Developmental Toxicology

This report summarizes the proceedings of a conference on quantitative methods for assessing the risks of developmental toxicants. The conference was planned by a subcommittee of the National Research Council's Committee on Risk Assessment Methodology ⁴ in conjunction with staff from several federal agencies, including the U.S. Environmental Protection Agency, U.S. Food and Drug Administration, U.S. Consumer Products Safety Commission, and Health and Welfare Canada. Issues discussed at the workshop included computerized techniques for hazard identification, use of human and animal data for defining risks in a clinical setting, relationships between end points in developmental toxicity testing, reference dose calculations for developmental toxicology, analysis of quantitative dose-response data, mechanisms of developmental toxicity, physiologically based pharmacokinetic models, and structure-activity relationships. Although a formal consensus was not sought, many participants favored the evolution of quantitative techniques for developmental toxicology risk assessment, including the replacement of lowest observed adverse effect levels (LOAELs) and no observed adverse effect levels (NOAELs) with the benchmark dose methodology.     

Incorporating Time‐Dose‐Response into Legionella Outbreak Models

A novel method was used to incorporate in vivo host–pathogen dynamics into a new robust outbreak model for legionellosis. Dose‐response and time‐dose‐response (TDR) models were generated for Legionella longbeachae exposure to mice via the intratracheal route using a maximum likelihood estimation approach. The best‐fit TDR model was then incorporated into two L. pneumophila outbreak models: an outbreak that occurred at a spa in Japan, and one that occurred in a Melbourne aquarium. The best‐fit TDR from the murine dosing study was the beta‐Poisson with exponential‐reciprocal dependency model, which had a minimized deviance of 32.9. This model was tested against other incubation distributions in the Japan outbreak, and performed consistently well, with reported deviances ranging from 32 to 35. In the case of the Melbourne outbreak, the exponential model with exponential dependency was tested against non‐time‐dependent distributions to explore the performance of the time‐dependent model with the lowest number of parameters. This model reported low minimized deviances around 8 for the Weibull, gamma, and lognormal exposure distribution cases. This work shows that the incorporation of a time factor into outbreak distributions provides models with acceptable fits that can provide insight into the in vivo dynamics of the host‐pathogen system.     

Estimation of Benchmark Dose of Lifetime Cadmium Intake for Adverse Renal Effects Using Hybrid Approach in Inhabitants of an Environmentally Exposed River Basin in Japan

The aim of this study is to estimate the reference level of lifetime cadmium intake (LCd) as the benchmark doses (BMDs) and their 95% lower confidence limits (BMDLs) for various renal effects by applying a hybrid approach. The participants comprised 3,013 (1,362 men and 1,651 women) and 278 (129 men and 149 women) inhabitants of the Cd‐polluted and nonpolluted areas, respectively, in the environmentally exposed Kakehashi River basin. Glucose, protein, aminonitrogen, metallothionein, and β₂‐microglobulin in urine were measured as indicators of renal dysfunction. The BMD and BMDL that corresponded to an additional risk of 5% were calculated with background risk at zero exposure set at 5%. The obtained BMDLs of LCd were 3.7 g (glucose), 3.2 g (protein), 3.7 g (aminonitrogen), 1.7 g (metallothionein), and 1.8 g (β₂‐microglobulin) in men and 2.9 g (glucose), 2.5 g (protein), 2.0 g (aminonitrogen), 1.6 g (metallothionein), and 1.3 g (β₂‐microglobulin) in women. The lowest BMDL was 1.7 g (metallothionein) and 1.3 g (β₂‐microglobulin) in men and women, respectively. The lowest BMDL of LCd (1.3 g) was somewhat lower than the representative threshold LCd (2.0 g) calculated in the previous studies. The obtained BMDLs may contribute to further discussion on the health risk assessment of cadmium exposure.     

Quantal Risk Assessment Database: A Database for Exploring Patterns in Quantal Dose‐Response Data in Risk Assessment and its Application to Develop Priors for Bayesian Dose‐Response Analysis

Quantitative risk assessments for physical, chemical, biological, occupational, or environmental agents rely on scientific studies to support their conclusions. These studies often include relatively few observations, and, as a result, models used to characterize the risk may include large amounts of uncertainty. The motivation, development, and assessment of new methods for risk assessment is facilitated by the availability of a set of experimental studies that span a range of dose‐response patterns that are observed in practice. We describe construction of such a historical database focusing on quantal data in chemical risk assessment, and we employ this database to develop priors in Bayesian analyses. The database is assembled from a variety of existing toxicological data sources and contains 733 separate quantal dose‐response data sets. As an illustration of the database's use, prior distributions for individual model parameters in Bayesian dose‐response analysis are constructed. Results indicate that including prior information based on curated historical data in quantitative risk assessments may help stabilize eventual point estimates, producing dose‐response functions that are more stable and precisely estimated. These in turn produce potency estimates that share the same benefit. We are confident that quantitative risk analysts will find many other applications and issues to explore using this database.