A Quantitative Risk Assessment of Waterborne Cryptosporidiosis in France Using Second-Order Monte Carlo Simulation

A pragmatic quantitative risk assessment (QRA) of the risks of waterborne Cryptosporidium parvum infection and cryptosporidiosis in immunocompetent and immunodeficient French populations is proposed. The model takes into account French specificities such as the French technique for oocyst enumeration performance and tap water consumption. The proportion of infective oocysts is based on literature review and expert knowledge. The probability of infection for a given number of ingested viable oocysts is modeled using the exponential dose-response model applied on published data from experimental infections in immunocompetent human volunteers challenged with the IOWA strain. Second-order Monte Carlo simulations are used to characterize the uncertainty and variability of the risk estimates. Daily risk of infection and illness for the immunocompetent and the immunodeficient populations are estimated according to the number of oocysts observed in a single storage reservoir water sample. As an example, the mean daily risk of infection in the immunocompetent population is estimated to be 1.08 x 10(-4) (95% confidence interval: [0.20 x 10(-4); 6.83 x 10(-4)]) when five oocysts are observed in a 100 L storage reservoir water sample. Annual risks of infection and disease are estimated from a set of oocyst enumeration results from distributed water samples, assuming a negative binomial distribution of day-to-day contamination variation. The model and various assumptions used in the model are fully explained and discussed. While caveats of this model are well recognized, this pragmatic QRA could represent a useful tool for the French Food Safety Agency (AFSSA) to define recommendations in case of water resource contamination by C. parvum whose infectivity is comparable to the IOWA strain.     

Daily Soil Ingestion Estimates for Children at a Superfund Site

Ingestion of contaminated soil by children may result in significant exposure to toxic substances at contaminated sites. Estimates of such exposure are based on extrapolation of short-term-exposure estimates to longer time periods. This article provides daily estimates of soil ingestion on 64 children between the ages of 1 and 4 residing at a Superfund site; these values are employed to estimate the distribution of 7-day average soil ingestion exposures (mean, 31 mg/day; median, 17 mg/day) at a contaminated site over different time periods. Best linear unbiased predictors of the 95th-percentile of soil ingestion over 7 days, 30 days, 90 days, and 365 days are 133 mg/day, 112 mg/day, 108 mg/day and 106 mg/day, respectively. Variance components estimates (excluding titanium and outliers, based on Tukey's far-out criteria) are given for soil ingestion between subjects (59 mg/day)2, between days on a subject (95 mg/day)2, and for uncertainty on a subject-day (132 mg/day)2. These results expand knowledge of potential exposure to contaminants among young children from soil ingestion at contaminated sites. They also provide basic distributions that serve as a starting point for use in Monte Carlo risk assessments.     

Principles of Good Practice for the Use of Monte Carlo Techniques in Human Health and Ecological Risk Assessments

We propose 14 principles of good practice to assist people in performing and reviewing probabilistic or Monte Carlo risk assessments, especially in the context of the federal and state statutes concerning chemicals in the environment. Monte Carlo risk assessments for hazardous waste sites that follow these principles will be easier to understand, will explicitly distinguish assumptions from data, and will consider and quantify effects that could otherwise lead to misinterpretation of the results. The proposed principles are neither mutually exclusive nor collectively exhaustive. We think and hope that these principles will evolve as new ideas arise and come into practice.     

Development of a Standard Soil-to-Skin Adherence Probability Density Function for Use in Monte Carlo Analyses of Dermal Exposure

It has recently been suggested that "standard" data distributions for key exposure variables should be developed wherever appropriate for use in probabilistic or "Monte Carlo" exposure analyses. Soil-on-skin adherence estimates represent an ideal candidate for development of a standard data distribution: There are several readily available studies which offer a consistent pattern of reported results, and more importantly, soil adherence to skin is likely to vary little from site-to-site. In this paper, we thoroughly review each of the published soil adherence studies with respect to study design, sampling, and analytical methods, and level of confidence in the reported results. Based on these studies, probability density functions (PDF) of soil adherence values were examined for different age groups and different sampling techniques. The soil adherence PDF developed from adult data was found to resemble closely the soil adherence PDF based on child data in terms of both central tendency (mean = 0.49 and 0.63 mg-soil/cm²-skin, respectively) and 95th percentile values (1.6 and 2.4 mg-soil/cm²-skin, respectively). Accordingly, a single, "standard" PDF is presented based on all data collected for all age groups. This standard PDF is lognormally distributed; the arithmetic mean and standard deviation are 0.52 ± 0.9 mg-soil/cm²-skin. Since our review of the literature indicates that soil adherence under environmental conditions will be minimally influenced by age, sex, soil type, or particle size, this PDF should be considered applicable to all settings. The 50th and 95th percentile values of the standard PDF (0.25 and 1.7 mg-soil/cm²-skin, respectively) are very similar to recent U.S. EPA estimates of "average" and "upper-bound" soil adherence (0.2 and 1.0 mg-soil/cm²-skin, respectively).     

Monte Carlo Simulation of Rodent Carcinogenicity Bioassays

In this paper we describe a simulation, by Monte Carlo methods, of the results of rodent carcinogenicity bioassays. Our aim is to study how the observed correlation between carcinogenic potency (beta or 1n2/TD50) and maximum tolerated dose (MTD) arises, and whether the existence of this correlation leads to an artificial correlation between carcinogenic potencies in rats and mice. The validity of the bioassay results depends upon, among other things, certain biases in the experimental design of the bioassays. These include selection of chemicals for bioassay and details of the experimental protocol, including dose levels. We use as variables in our simulation the following factors: (1) dose group size, (2) number of dose groups, (3) tumor rate in the control (zero-dose) group, (4) distribution of the MTD values of the group of chemicals as specified by the mean and standard deviation, (5) the degree of correlation between beta and the MTD, as given by the standard deviation of the random error term in the linear regression of log beta on log (1/MTD), and (6) an upper limit on the number of animals with tumors. Monte Carlo simulation can show whether the information present in the existing rodent bioassay database is sufficient to reject the validity of the proposed interspecies correlations at a given level of stringency. We hope that such analysis will be useful for future bioassay design, and more importantly, for discussion of the whole NCI/NTP program.     

Recommendations on the Testing and Use of Pseudo-Random Number Generators Used in Monte Carlo Analysis for Risk Assessment

Monte Carlo simulation requires a pseudo-random number generator with good statistical properties. Linear congruential generators (LCGs) are the most popular and well-studied computer method for generating pseudo-random numbers used in Monte Carlo studies. High quality LCGs are available with sufficient statistical quality to satisfy all but the most demanding needs of risk assessors. However, because of the discrete, deterministic nature of LCGs, it is important to evaluate the randomness and uniformity of the specific pseudo-random number subsequences used in important risk assessments. Recommended statistical tests for uniformity and randomness include the Kolmogorov-Smirnov test, extreme values test, and the runs test, including runs above and runs below the mean tests. Risk assessors should evaluate the stability of their risk model's output statistics, paying particular attention to instabilities in the mean and variance. When instabilities in the mean and variance are observed, more stable statistics, e.g., percentiles, should be reported. Analyses should be repeated using several non-overlapping pseudo-random number subsequences. More simulations than those traditionally used are also recommended for each analysis.     

Monte Carlo Sampling for Generalized Knowledge Dependence with Application to Human Reliability

A general discussion of knowledge dependence in risk calculations shows that the assumption of independence underlying standard Monte Carlo simulation in uncertainty analysis is frequently violated. A model is presented for performing Monte Carlo simulation when the variabilities of the component failure probabilities are either negatively or positively coupled. The model is applied to examples in human reliability analysis and the results are compared to the results of Sandia Laboratories as published in the Peer Review Study and to recalculations using more recent methods of uncertainty analysis.     

存在退保时分红寿险定价的最小二乘蒙特卡罗模拟

分红型人寿保险保单可以视作由三部分构成:固定收益债券、分红权和退保权.退保权的存在使保单具有美式期权的性质,给定价带来困难.本文用最小二乘蒙特卡罗模拟,建立了计算保单价值的模型,给出了模拟计算结果.     

Effects of Uncertainties on Exposure Estimates to Methylmercury: A Monte Carlo Analysis of Exposure Biomarkers versus Dietary Recall Estimation

This article presents a general model for estimating population heterogeneity and "lack of knowledge" uncertainty in methylmercury (MeHg) exposure assessments using two-dimensional Monte Carlo analysis. Using data from fish-consuming populations in Bangladesh, Brazil, Sweden, and the United Kingdom, predictive model estimates of dietary MeHg exposures were compared against those derived from biomarkers (i.e., [Hg]hair and [Hg]blood). By disaggregating parameter uncertainty into components (i.e., population heterogeneity, measurement error, recall error, and sampling error) estimates were obtained of the contribution of each component to the overall uncertainty. Steady-state diet:hair and diet:blood MeHg exposure ratios were estimated for each population and were used to develop distributions useful for conducting biomarker-based probabilistic assessments of MeHg exposure. The 5th and 95th percentile modeled MeHg exposure estimates around mean population exposure from each of the four study populations are presented to demonstrate lack of knowledge uncertainty about a best estimate for a true mean. Results from a U.K. study population showed that a predictive dietary model resulted in a 74% lower lack of knowledge uncertainty around a central mean estimate relative to a hair biomarker model, and also in a 31% lower lack of knowledge uncertainty around central mean estimate relative to a blood biomarker model. Similar results were obtained for the Brazil and Bangladesh populations. Such analyses, used here to evaluate alternative models of dietary MeHg exposure, can be used to refine exposure instruments, improve information used in site management and remediation decision making, and identify sources of uncertainty in risk estimates.     

半参数空间ZISF的估计、分类与蒙特卡罗模拟

零无效率随机前沿模型（ZISF）包含随机前沿模型和回归模型,两模型各有一定的发生概率,适用于技术无效生产单元和技术有效生产单元同时存在的情形。本文在ZISF的生产函数中引入空间效应和非参函数,并假设回归模型的发生概率为非参函数,构建了半参数空间ZISF。该模型可有效避免忽略空间效应导致的有偏且不一致估计量,也避免了线性模型的拟合不足。本文对非参函数采用B样条逼近,使用极大似然方法和JLMS法分别估计参数和技术效率。蒙特卡罗结果表明：①本文方法的估计精度和分类精度均较高。随着样本容量的增大,精度增加。②忽略空间效应或者非参数效应,估计精度和分类精度降低,文中模型有存在必要性。③忽略发生概率的非参数效应会严重降低估计和分类精度,远大于忽略生产函数的非参数效应的影响。     

Sensitivity Analysis, Monte Carlo Risk Analysis, and Bayesian Uncertainty Assessment

Standard statistical methods understate the uncertainty one should attach to effect estimates obtained from observational data. Among the methods used to address this problem are sensitivity analysis, Monte Carlo risk analysis (MCRA), and Bayesian uncertainty assessment. Estimates from MCRAs have been presented as if they were valid frequentist or Bayesian results, but examples show that they need not be either in actual applications. It is concluded that both sensitivity analyses and MCRA should begin with the same type of prior specification effort as Bayesian analysis.     

A Trichloroethylene Risk Assessment Using a Monte Carlo Analysis of Parameter Uncertainty in Conjunction with Physiologically-Based Pharmacokinetic Modeling

A Monte Carlo simulation is incorporated into a risk assessment for trichloroethylene (TCE) using physiologically-based pharmacokinetic (PBPK) modeling coupled with the linearized multistage model to derive human carcinogenic risk extrapolations. The Monte Carlo technique incorporates physiological parameter variability to produce a statistically derived range of risk estimates which quantifies specific uncertainties associated with PBPK risk assessment approaches. Both inhalation and ingestion exposure routes are addressed. Simulated exposure scenarios were consistent with those used by the Environmental Protection Agency (EPA) in their TCE risk assessment. Mean values of physiological parameters were gathered from the literature for both mice (carcinogenic bioassay subjects) and for humans. Realistic physiological value distributions were assumed using existing data on variability. Mouse cancer bioassay data were correlated to total TCE metabolized and area-under-the-curve (blood concentration) trichloroacetic acid (TCA) as determined by a mouse PBPK model. These internal dose metrics were used in a linearized multistage model analysis to determine dose metric values corresponding to 10^-6 lifetime excess cancer risk. Using a human PBPK model, these metabolized doses were then extrapolated to equivalent human exposures (inhalation and ingestion). The Monte Carlo iterations with varying mouse and human physiological parameters produced a range of human exposure concentrations producing a 10^-6 risk.     

Uncertainty and Variability in Human Exposures to Soil Contaminants Through Home-Grown Food: A Monte Carlo Assessment

This paper presents a general model for exposure to homegrown foods that is used with a Monte Carlo analysis to determine the relative contributions of variability (Type A uncertainty) and true uncertainty (Type B uncertainty) to the overall variance in prediction of the dose-to-concentration ratio. Although classification of exposure inputs as uncertain or variable is somewhat subjective, food consumption rates and exposure duration are judged to have a predicted variance that is dominated by variability among individuals by age, income, culture, and geographical region. Whereas, biotransfer factors and partition factors are inputs that, to a large extent, involve uncertainty. Using ingestion of fruits, vegetables, grains, dairy products, and meat and soils assumed to be contaminated by hexachlorbenzene (HCB) and benzo(a)pyrene (BaP) as cases studies, a Monte Carlo analysis is used to explore the relative contribution of uncertainty and variability to overall variance in the estimated distribution of potential dose within the population that consumes homegrown foods. It is found that, when soil concentrations are specified, variances in ratios of dose-to-concentration for HCB are equally attributable to uncertainty and variability, whereas for BaP, variance in these ratios is dominated by true uncertainty.     

基于幂风险谱和蒙特卡洛模拟的贷款优化配置模型

极端风险对于银行资产配置至关重要,尤其是次贷危机之后尾部风险以惨重的代价引起了金融机构的重视,传统条件风险价值CVaR、风险价值VaR不能有效度量尾部极端风险,因此本文基于幂风险谱和蒙特卡洛模拟构建了贷款组合优化配置模型,同时控制尾部极端风险和信用风险。本文一是通过损失-X_i越大、其风险权重φ_i也就越大的思路,构建幂风险谱PSR （Power Spectral Risk）最小的目标函数对极端风险进行控制,即弥补了CVaR同等看待尾部风险、忽略风险较大的损失应予以更大权重的不足,也同时弥补了VaR仅提供某一置信水平下资产损失的最大值、无法反映一旦超过这一数值的可能损失的弊端。二是通过蒙特卡洛模拟信用等级迁移引起贷款收益的变化情景,并以信用等级迁移后贷款组合损失越大、则风险厌恶权重越大的思路构建幂风险谱PSR最小为目标函数,以贷款组合的收益大于目标收益为约束,构建贷款优化配置模型,改变了现有研究贷款配置时没有同时控制信用风险和尾部风险的不足。对比分析结果表明：本文模型能够实现更高的收益风险比,即在单位幂风险谱PSR下能够实现较高的收益。     

Exposure Assessment: Methods of Analysis for Environmental Carcinogens1

Human populations are exposed to environmental carcinogens in both indoor and outdoor atmospheres. Recent studies indicate that pollutant concentrations are generally higher in indoor atmospheres than in outdoor. Environmental pollutants that occur in indoor air from a variety of sources include radon, asbestos, organic and inorganic compounds, and certain particles (e.g., tobacco smoke). Some of the gases or vapors are adsorbed on suspended particulate matter, whereas others exist entirely in the gas phase or are distributed between the latter and a particle-bound state. Because of differences in chemical and physical properties, each class of carcinogens generally requires different sampling and analytical methods. In addition, a single indoor environment may contain a wide variety of air pollutants from different sources. Unfortunately, no single best approach currently exists for the quantitative determination of such complex mixtures and, for practical reasons, only the more toxic or the more abundant pollutants are usually measured. This paper summarizes the currently available monitoring methods for selected environmental pollutants found in indoor atmospheres. In addition, some possible sources for those pollutants are identified.     

Mlonte Carlo Techniques for Quantitative Uncertainty Analysis in Public Health Risk Assessments

Most public health risk assessments assume and combine a series of average, conservative, and worst-case values to derive a conservative point estimate of risk. This procedure has major limitations. This paper demonstrates a new methodology for extended uncertainty analyses in public health risk assessments using Monte Carlo techniques. The extended method begins as do some conventional methods--with the preparation of a spreadsheet to estimate exposure and risk. This method, however, continues by modeling key inputs as random variables described by probability density functions (PDFs). Overall, the technique provides a quantitative way to estimate the probability distributions for exposure and health risks within the validity of the model used. As an example, this paper presents a simplified case study for children playing in soils contaminated with benzene and benzo(a)pyrene (BaP).     

Estimates of Lifetime-Absorbed Daily Doses from the Use of Personal-Care Products Containing Polyacrylamide: A Monte Carlo Analysis

Estimates of the lifetime-absorbed daily dose (LADD) of acrylamide resulting from use of representative personal-care products containing polyacrylamides have been developed. All of the parameters that determine the amount of acrylamide absorbed by an individual vary from one individual to another. Moreover, for some parameters there is uncertainty as to which is the correct or representative value from a range of values. Consequently, the parameters used in the estimation of the LADD of acrylamide from usage of a particular product type (e.g., deodorant, makeup, etc.) were represented by distributions evaluated using Monte Carlo analyses.((1-4)) From these data, distributions of values for key parameters, such as the amount of acrylamide in polyacrylamide, absorption fraction, etc., were defined and used to provide a distribution of LADDs for each personal-care product. The estimated total acrylamide LADD (across all products) for males and females at the median, mean, and 95th percentile of the distribution of individual LADD values were 4.7 x 10(-8), 2.3 x 10(-7), and 7.3 x 10(-7) mg/kg/day for females and 3.6 x 10(-8), 1.7 x 10(-7), and 5.4 x 10(-7) mg/kg/day for males. The ratio of the LADDs to risk-specific dose corresponding to a target risk level of 1 x 10(-5), the acceptable risk level for this investigation, derived using approaches typically used by the FDA, the USEPA, and proposed for use by the European Union (EU) were also calculated. All ratios were well below 1, indicating that all the extra lifetime cancer risk from the use of polyacrylamide-containing personal-care products, in the manner assumed in this assessment, are well below acceptable levels. Even if it were assumed that an individual used all of the products together, the estimated LADD would still provide a dose that was well below the acceptable risk levels.     

Quantifying the Distribution of Inhalation Exposure in Human Populations: 2. Distributions of Time Spent by Adults, Adolescents, and Children at Home, at Work, and at School

Using distributions of time spent at various ventilation levels, ranges of inhalation exposure in the population can be established. Distributions of exposure time were determined using results of a study by the California Air Resources Board (CARB) which focused on time spent by humans participating in various activities and the locations where the activities occurred. The daily at-home activities from the CARB study were assigned to one of three ventilation levels, generating aggregate time periods. Distinct age and gender populations were identified, and distributions for aggregate time were established for these populations at each of the ventilation levels. In addition to aggregate time spent at home, distributions for various ages and genders were established for aggregate time spent at school and work. By combining distributions of aggregate time with corresponding ventilation rates, the distribution of inhalation rates can be established for at home, at work, and at school exposures.     

Interpretation of Monte Carlo Simulations Using Parameter Space Plots

Multimodal distribution functions that result from Monte Carlo simulations can be interpreted by superimposing joint probability density functions onto the contour space of the simulated calculations. The method is demonstrated by analysis of the pathway of a radioactive groundwater contaminant using an analytical solution to the transport equation. Simulated concentrations at a fixed time and distance produce multimodal histograms, which are understood with reference to the parameter space for the two random variables-velocity and dispersivity. Numerical integration under the joint density function up to the contour of the analytical solution gives the probability of contaminant exceeding a target concentration. This technique is potentially more efficient than Monte Carlo simulation for low probability events. Visualization of parameter space is restricted to two random variables. Nevertheless, analyzing the two most pertinent random variables in a simulation might still offer insights into the multimodal nature of output histograms.