Variation in Cancer Risk Estimates for Exposure to Powerline Frequency Electromagnetic Fields: A Meta-Analysis Comparing EMF Measurement Methods

We used meta-analysis to synthesize the findings from eleven case-control studies on cancer risks in humans exposed to 50-60 Hertz powerline electromagnetic fields (EMFs). Pooled estimates of risk are derived for different EMF measurement methods and types of cancer. EMF measurement methods are classified as: wiring configuration codes, distance to power distribution equipment, spot measurements of magnetic fields, and calculated indices based on distance to power distribution equipment and historic load data. Pooled odds ratios depicting the risk of cancer by each measurement type are presented for all cancers combined, leukemia for all age groups and childhood leukemia. The wire code measurement technique was associated with a significantly increased risk for all three cancer types, while spot measures consistently showed non-significant odds ratios. Distance measures and the calculated indices produced risk estimates which were significant only for leukemia.     

Limitations to Empirical Extrapolation Studies: The Case of BMD Ratios

Extrapolation relationships are of keen interest to chemical risk assessment in which they play a prominent role in translating experimentally derived (usually in animals) toxicity estimates into estimates more relevant to human populations. A standard approach for characterizing each extrapolation relies on ratios of pre-existing toxicity estimates. Applications of this "ratio approach" have overlooked several sources of error. This article examines the case of ratios of benchmark doses, trying to better understand their informativeness. The approach involves mathematically modeling the process by which the ratios are generated in practice. Both closed form and simulation-based models of this "data-generating process" (DGP) are developed, paying special attention to the influence of experimental design. The results show the potential for significant limits to informativeness, and revealing dependencies. Future applications of the ratio approach should take imprecision and bias into account. Bootstrap techniques are recommended for gauging imprecision, but more complicated techniques will be required for gauging bias (and capturing dependencies). Strategies for mitigating the errors are suggested.     

Compressed Natural Gas Bus Safety: A Quantitative Risk Assessment

This study assesses the fire safety risks associated with compressed natural gas (CNG) vehicle systems, comprising primarily a typical school bus and supporting fuel infrastructure. The study determines the sensitivity of the results to variations in component failure rates and consequences of fire events. The components and subsystems that contribute most to fire safety risk are determined. Finally, the results are compared to fire risks of the present generation of diesel-fueled school buses. Direct computation of the safety risks associated with diesel-powered vehicles is possible because these are mature technologies for which historical performance data are available. Because of limited experience, fatal accident data for CNG bus fleets are minimal. Therefore, this study uses the probabilistic risk assessment (PRA) approach to model and predict fire safety risk of CNG buses. Generic failure data, engineering judgments, and assumptions are used in this study. This study predicts the mean fire fatality risk for typical CNG buses as approximately 0.23 fatalities per 100-million miles for all people involved, including bus passengers. The study estimates mean values of 0.16 fatalities per 100-million miles for bus passengers only. Based on historical data, diesel school bus mean fire fatality risk is 0.091 and 0.0007 per 100-million miles for all people and bus passengers, respectively. One can therefore conclude that CNG buses are more prone to fire fatality risk by 2.5 times that of diesel buses, with the bus passengers being more at risk by over two orders of magnitude. The study estimates a mean fire risk frequency of 2.2 x 10(-5) fatalities/bus per year. The 5% and 95% uncertainty bounds are 9.1 x 10(-6) and 4.0 x 10(-5), respectively. The risk result was found to be affected most by failure rates of pressure relief valves, CNG cylinders, and fuel piping.     

Model Averaging in Microbial Risk Assessment Using Fractional Polynomials

The alleviation of food-borne diseases caused by microbial pathogen remains a great concern in order to ensure the well-being of the general public. The relation between the ingested dose of organisms and the associated infection risk can be studied using dose-response models. Traditionally, a model selected according to a goodness-of-fit criterion has been used for making inferences. In this article, we propose a modified set of fractional polynomials as competitive dose-response models in risk assessment. The article not only shows instances where it is not obvious to single out one best model but also illustrates that model averaging can best circumvent this dilemma. The set of candidate models is chosen based on biological plausibility and rationale and the risk at a dose common to all these models estimated using the selected models and by averaging over all models using Akaike's weights. In addition to including parameter estimation inaccuracy, like in the case of a single selected model, model averaging accounts for the uncertainty arising from other competitive models. This leads to a better and more honest estimation of standard errors and construction of confidence intervals for risk estimates. The approach is illustrated for risk estimation at low dose levels based on Salmonella typhi and Campylobacter jejuni data sets in humans. Simulation studies indicate that model averaging has reduced bias, better precision, and also attains coverage probabilities that are closer to the 95% nominal level compared to best-fitting models according to Akaike information criterion.     

An Evaluation of Risk Estimation Procedures for Mixtures of Carcinogens

The estimation of health risks from exposure to a mixture of chemical carcinogens is generally based on the combination of information from several available single compound studies. The current practice of directly summing the upper bound risk estimates of individual carcinogenic components as an upper bound on the total risk of a mixture is known to be generally too conservative. Gaylor and Chen (1996, Risk Analysis) proposed a simple procedure to compute an upper bound on the total risk using only the upper confidence limits and central risk estimates of individual carcinogens. The Gaylor-Chen procedure was derived based on an underlying assumption of the normality for the distributions of individual risk estimates. In this paper we evaluated the Gaylor-Chen approach in terms of the coverage probability. The performance of the Gaylor-Chen approach in terms the coverages of the upper confidence limits on the true risks of individual carcinogens. In general, if the coverage probabilities for the individual carcinogens are all approximately equal to the nominal level, then the Gaylor-Chen approach should perform well. However, the Gaylor-Chen approach can be conservative or anti-conservative if some or all individual upper confidence limit estimates are conservative or anti-conservative.     

Potential Airborne Asbestos Exposure and Risk Associated with the Historical Use of Cosmetic Talcum Powder Products

Over time, concerns have been raised regarding the potential for human exposure and risk from asbestos in cosmetic‐talc–containing consumer products. In 1985, the U.S. Food and Drug Administration (FDA) conducted a risk assessment evaluating the potential inhalation asbestos exposure associated with the cosmetic talc consumer use scenario of powdering an infant during diapering, and found that risks were below levels associated with background asbestos exposures and risk. However, given the scope and age of the FDA's assessment, it was unknown whether the agency's conclusions remained relevant to current risk assessment practices, talc application scenarios, and exposure data. This analysis updates the previous FDA assessment by incorporating the current published exposure literature associated with consumer use of talcum powder and using the current U.S. Environmental Protection Agency's (EPA) nonoccupational asbestos risk assessment approach to estimate potential cumulative asbestos exposure and risk for four use scenarios: (1) infant exposure during diapering; (2) adult exposure from infant diapering; (3) adult exposure from face powdering; and (4) adult exposure from body powdering. The estimated range of cumulative asbestos exposure potential for all scenarios (assuming an asbestos content of 0.1%) ranged from 0.0000021 to 0.0096 f/cc‐yr and resulted in risk estimates that were within or below EPA's acceptable target risk levels. Consistent with the original FDA findings, exposure and corresponding health risk in this range were orders of magnitude below upper‐bound estimates of cumulative asbestos exposure and risk at ambient levels, which have not been associated with increased incidence of asbestos‐related disease.     

Model Uncertainty and Risk Estimation for Experimental Studies of Quantal Responses

Experimental animal studies often serve as the basis for predicting risk of adverse responses in humans exposed to occupational hazards. A statistical model is applied to exposure-response data and this fitted model may be used to obtain estimates of the exposure associated with a specified level of adverse response. Unfortunately, a number of different statistical models are candidates for fitting the data and may result in wide ranging estimates of risk. Bayesian model averaging (BMA) offers a strategy for addressing uncertainty in the selection of statistical models when generating risk estimates. This strategy is illustrated with two examples: applying the multistage model to cancer responses and a second example where different quantal models are fit to kidney lesion data. BMA provides excess risk estimates or benchmark dose estimates that reflects model uncertainty.     

Vehicle Emission Unit Risk Factors for Transportation Risk Assessments

When the transportation risk posed by shipments of hazardous chemical and radioactive materials is being assessed, it is necessary to evaluate the risks associated with both vehicle emissions and cargo-related risks. Diesel exhaust and fugitive dust emissions from vehicles transporting hazardous shipments lead to increased air pollution, which increases the risk of latent fatalities in the affected population along the transport route. The estimated risk from these vehicle-related sources can often be as large or larger than the estimated risk associated with the material being transported. In this paper, data from the U.S. Environmental Protection Agency's Motor Vehicle-Related Air Toxics Study are first used to develop latent cancer fatality estimates per kilometer of travel in rural and urban areas for all diesel truck classes. These unit risk factors are based on studies investigating the carcinogenic nature of diesel exhaust. With the same methodology, the current per-kilometer latent fatality risk factor used in transportation risk assessments for heavy diesel trucks in urban areas is revised and the analysis expanded to provide risk factors for rural areas and all diesel truck classes. These latter fatality estimates may include, but are not limited to, cancer fatalities and are based primarily on the most recent epidemiological data available on mortality rates associated with ambient air PM-10 concentrations.     

A Risk Assessment for Occupational Acrylonitrile Exposure Using Epidemiology Data

The extensive data from the Blair et al.((1)) epidemiology study of occupational acrylonitrile exposure among 25460 workers in eight plants in the United States provide an excellent opportunity to update quantitative risk assessments for this widely used commodity chemical. We employ the semiparametric Cox relative risk (RR) regression model with a cumulative exposure metric to model cause-specific mortality from lung cancer and all other causes. The separately estimated cause-specific cumulative hazards are then combined to provide an overall estimate of age-specific mortality risk. Age-specific estimates of the additional risk of lung cancer mortality associated with several plausible occupational exposure scenarios are obtained. For age 70, these estimates are all markedly lower than those generated with the cancer potency estimate provided in the USEPA acrylonitrile risk assessment.((2)) This result is consistent with the failure of recent occupational studies to confirm elevated lung cancer mortality among acrylonitrile-exposed workers as was originally reported by O'Berg,((3)) and it calls attention to the importance of using high-quality epidemiology data in the risk assessment process.     

Determinants of Risk Perceptions of a Hazardous Waste Site

A before-stimulus-after quasi-experimental design is used to assess the factors relating to risk perceptions of a hazardous waste site. First, a pretest obtains measures of attitudes and beliefs about hazardous waste and waste sites. Second, a detailed hypothetical "Superfund" scenario, including a complex cleanup plan, is introduced. Finally, indices of health risk estimates, trust, knowledge, and other pertinent beliefs are obtained. Levels of concern, both before and after cleanup, are the dependent variables. Independent variables include risk management options, health risk estimates, trust, and five sociodemographic characteristics. Concern is extremely high prior to cleanup and moderately high after cleanup. Concern is a clear function of health risk estimates. Toxic chemicals from waste sites are viewed as a major cause of multiple health problems, especially cancers. Accurate health risk estimates moderate fears and are linked to levels of education. Education, however, does not explain concern. Trust is a major factor explaining concern and health risk estimates. The implications of these findings for risk communication is discussed.     

Addressing Nonlinearity in the Exposure-Response Relationship for a Genotoxic Carcinogen: Cancer Potency Estimates for Ethylene Oxide

Ethylene oxide (EO) has been identified as a carcinogen in laboratory animals. Although the precise mechanism of action is not known, tumors in animals exposed to EO are presumed to result from its genotoxicity. The overall weight of evidence for carcinogenicity from a large body of epidemiological data in the published literature remains limited. There is some evidence for an association between EO exposure and lympho/hematopoietic cancer mortality. Of these cancers, the evidence provided by two large cohorts with the longest follow-up is most consistent for leukemia. Together with what is known about human leukemia and EO at the molecular level, there is a body of evidence that supports a plausible mode of action for EO as a potential leukemogen. Based on a consideration of the mode of action, the events leading from EO exposure to the development of leukemia (and therefore risk) are expected to be proportional to the square of the dose. In support of this hypothesis, a quadratic dose-response model provided the best overall fit to the epidemiology data in the range of observation. Cancer dose-response assessments based on human and animal data are presented using three different assumptions for extrapolating to low doses: (1) risk is linearly proportionate to dose; (2) there is no appreciable risk at low doses (margin-of-exposure or reference dose approach); and (3) risk below the point of departure continues to be proportionate to the square of the dose. The weight of evidence for EO supports the use of a nonlinear assessment. Therefore, exposures to concentrations below 37 microg/m3 are not likely to pose an appreciable risk of leukemia in human populations. However, if quantitative estimates of risk at low doses are desired and the mode of action for EO is considered, these risks are best quantified using the quadratic estimates of cancer potency, which are approximately 3.2- to 32-fold lower, using alternative points of departure, than the linear estimates of cancer potency for EO. An approach is described for linking the selection of an appropriate point of departure to the confidence in the proposed mode of action. Despite high confidence in the proposed mode of action, a small linear component for the dose-response relationship at low concentrations cannot be ruled out conclusively. Accordingly, a unit risk value of 4.5 x 10(-8) (microg/m3)(-1) was derived for EO, with a range of unit risk values of 1.4 x 10(-8) to 1.4 x 10(-7) (microg/m3)(-1) reflecting the uncertainty associated with a theoretical linear term at low concentrations.     

Assessing Cancer Risks from Short-Term Exposures in Children

For the vast majority of chemicals that have cancer potency estimates on IRIS, the underlying database is deficient with respect to early-life exposures. This data gap has prevented derivation of cancer potency factors that are relevant to this time period, and so assessments may not fully address children's risks. This article provides a review of juvenile animal bioassay data in comparison to adult animal data for a broad array of carcinogens. This comparison indicates that short-term exposures in early life are likely to yield a greater tumor response than short-term exposures in adults, but similar tumor response when compared to long-term exposures in adults. This evidence is brought into a risk assessment context by proposing an approach that: (1) does not prorate children's exposures over the entire life span or mix them with exposures that occur at other ages; (2) applies the cancer slope factor from adult animal or human epidemiology studies to the children's exposure dose to calculate the cancer risk associated with the early-life period; and (3) adds the cancer risk for young children to that for older children/adults to yield a total lifetime cancer risk. The proposed approach allows for the unique exposure and pharmacokinetic factors associated with young children to be fully weighted in the cancer risk assessment. It is very similar to the approach currently used by U.S. EPA for vinyl chloride. The current analysis finds that the database of early life and adult cancer bioassays supports extension of this approach from vinyl chloride to other carcinogens of diverse mode of action. This approach should be enhanced by early-life data specific to the particular carcinogen under analysis whenever possible.     

Unveiling Variability and Uncertainty for Better Science and Decisions on Cancer Risks from Environmental Chemicals

The National Research Council 2009 “Silver Book” panel report included a recommendation that the U.S. Environmental Protection Agency (EPA) should increase all of its chemical carcinogen (CC) potency estimates by ～7‐fold to adjust for a purported median‐vs.‐mean bias that I recently argued does not exist (Bogen KT. “Does EPA underestimate cancer risks by ignoring susceptibility differences?,” Risk Analysis, 2014; 34(10):1780–1784). In this issue of the journal, my argument is critiqued for having flaws concerning: (1) intent, bias, and conservatism of EPA estimates of CC potency; (2) bias in potency estimates derived from epidemiology; and (3) human‐animal CC‐potency correlation. However, my argument remains valid, for the following reasons. (1) EPA's default approach to estimating CC risks has correctly focused on bounding average (not median) individual risk under a genotoxic mode‐of‐action (MOA) assumption, although pragmatically the approach leaves both inter‐individual variability in CC–susceptibility, and widely varying CC‐specific magnitudes of fundamental MOA uncertainty, unquantified. (2) CC risk estimates based on large epidemiology studies are not systematically biased downward due to limited sampling from broad, lognormal susceptibility distributions. (3) A good, quantitative correlation is exhibited between upper‐bounds on CC‐specific potency estimated from human vs. animal studies (n = 24, r = 0.88, p = 2 × 10^?8). It is concluded that protective upper‐bound estimates of individual CC risk that account for heterogeneity in susceptibility, as well as risk comparisons informed by best predictions of average‐individual and population risk that address CC‐specific MOA uncertainty, should each be used as separate, complimentary tools to improve regulatory decisions concerning low‐level, environmental CC exposures.     

The EPA Health Risk Assessment of Methylcyclopentadienyl Manganese Tricarbonyl (MMT)

This paper describes the U.S. Environmental Protection Agency's assessment of potential health risks associated with the possible widespread use of a manganese (Mn)-based fuel additive, methylcyclopentadienyl manganese tricarbonyl (MMT). This assessment was significant in several respects and may be instructive in identifying certain methodological issues of general relevance to risk assessment. A major feature of the inhalation health risk assessment was the derivation of Mn inhalation reference concentration (RfC) estimates using various statistical approaches, including benchmark dose and Bayesian analyses. The exposure assessment component used data from the Particle Total Exposure Assessment Methodology (PTEAM) study and other sources to estimate personal exposure levels of particulate Mn attributable to the permitted use of MMT in leaded gasoline in Riverside, CA, at the time of the PTEAM study; on this basis it was then possible to predict a distribution of possible future exposure levels associated with the use of MMT in all unleaded gasoline. Qualitative as well as quantitative aspects of the risk characterization are summarized, along with inherent uncertainties due to data limitations.     

The Benefits of Probabilistic Exposure Assessment: Three Case Studies Involving Contaminated Air,Water, and Soil1

Probabilistic risk assessments are enjoying increasing popularity as a tool to characterize the health hazards associated with exposure to chemicals in the environment. Because probabilistic analyses provide much more information to the risk manager than standard “point” risk estimates, this approach has generally been heralded as one which could significantly improve the conduct of health risk assessments. The primary obstacles to replacing point estimates with probabilistic techniques include a general lack of familiarity with the approach and a lack of regulatory policy and guidance. This paper discusses some of the advantages and disadvantages of the point estimate vs. probabilistic approach. Three case studies are presented which contrast and compare the results of each. The first addresses the risks associated with household exposure to volatile chemicals in tapwater. The second evaluates airborne dioxin emissions which can enter the food-chain. The third illustrates how to derive health-based cleanup levels for dioxin in soil. It is shown that, based on the results of Monte Carlo analyses of probability density functions (PDFs), the point estimate approach required by most regulatory agencies will nearly always overpredict the risk for the 95^th percentile person by a factor of up to 5. When the assessment requires consideration of 10 or more exposure variables, the point estimate approach will often predict risks representative of the 99.9^th percentile person rather than the 50^th or 95^th percentile person. This paper recommends a number of data distributions for various exposure variables that we believe are now sufficiently well understood to be used with confidence in most exposure assessments. A list of exposure variables that may require additional research before adequate data distributions can be developed are also discussed.     

Accounting for Both Random Errors and Systematic Errors in Uncertainty Propagation Analysis of Computer Models Involving Experimental Measurements with Monte Carlo Methods

A Monte Carlo method is presented to study the effect of systematic and random errors on computer models mainly dealing with experimental data. It is a common assumption in this type of models (linear and nonlinear regression, and nonregression computer models) involving experimental measurements that the error sources are mainly random and independent with no constant background errors (systematic errors). However, from comparisons of different experimental data sources evidence is often found of significant bias or calibration errors. The uncertainty analysis approach presented in this work is based on the analysis of cumulative probability distributions for output variables of the models involved taking into account the effect of both types of errors. The probability distributions are obtained by performing Monte Carlo simulation coupled with appropriate definitions for the random and systematic errors. The main objectives are to detect the error source with stochastic dominance on the uncertainty propagation and the combined effect on output variables of the models. The results from the case studies analyzed show that the approach is able to distinguish which error type has a more significant effect on the performance of the model. Also, it was found that systematic or calibration errors, if present, cannot be neglected in uncertainty analysis of models dependent on experimental measurements such as chemical and physical properties. The approach can be used to facilitate decision making in fields related to safety factors selection, modeling, experimental data measurement, and experimental design.     

Effect of detection methods on risk estimates of exposure to biosolids-associated human enteric viruses

This study illustrates the effect of virus detection methods on estimates of risks of infection of biosolids-associated viruses for occupational workers and residential population during a hypothetical exposure of biosolids. Five gastroenteritis-associated human enteric viruses--enteroviruses (echovirus-12, enteroviruse types 68-71), adenoviruses, rotaviruses, and noroviruses genotype--I-were considered to represent human enteric viruses for risk estimation purposes. Ingested viral doses were calculated using literature-reported total infectious virus concentrations (based on BGM and A549 cell lines) and genome copies (GCs) in Michigan dewatered and class B biosolids. Cell-line-based infectivity parameters (i.e., ratio of total infectious virus concentration to GCs) were developed for different viruses in biosolids to use GCs for calculating ingested viral dose, addressing the issue of integration of molecular methods with biosolids-based virus risk assessment. Use of virus concentrations from molecular methods (with and without using cell-line-based infectivity parameter) resulted in higher risk estimates than culture methods, indicating the effect of the virus detection method on risk estimates. Further, use of virus concentrations from A549 cell lines resulted in higher risk estimates compared to those from BGM cell lines, suggesting the need for a proper choice of cell lines in determining infectious viral dose. The Monte Carlo uncertainty analyses of estimates for risk of infection due to enteroviruses showed that enteroviruses concentration was the most important parameter influencing risk estimates, indicating the need for reducing associated uncertainty. More work is required to develop cell-line-based infectivity parameters for different virus concentration levels and sample matrix types using a cut-off-based approach.     

Probabilistic Health Risk Assessment for Exposures to Estuary Sediments and Biota Contaminated with Polychlorinated Biphenyls, Polychlorinated Terphenyls and Other Toxic Substances

A Latin Hypercube probabilistic risk assessment methodology was employed in the assessment of health risks associated with exposures to contaminated sediment and biota in an estuary in the Tidewater region of Virginia. The primary contaminants were polychlorinated biphenyls (PCBs), polychlorinated terphenyls (PCTs), polynuclear aromatic hydrocarbons (PAHs), and metals released into the estuary from a storm sewer system. The exposure pathways associated with the highest contaminant intake and risks were dermal contact with contaminated sediment and ingestion of contaminated aquatic and terrestrial biota from the contaminated area. As expected, all of the output probability distributions of risk were highly skewed, and the ratios of the expected value (mean) to median risk estimates ranged from 1.4 to 14.8 for the various exposed populations. The 99th percentile risk estimates were as much as two orders of magnitude above the mean risk estimates. For the sediment exposure pathways, the stability of the median risk estimates was found to be much greater than the stability of the expected value risk estimates. The interrun variability in the median risk estimate was found to be +/-1.9% at 3000 iterations. The interrun stability of the mean risk estimates was found to be approximately equal to that of the 95th percentile estimates at any number of iterations. The variation in neither contaminant concentrations nor any other single input variable contributed disproportionately to the overall simulation variance. The inclusion or exclusion of spatial correlations among contaminant concentrations in the simulation model did not significantly effect either the magnitude or the variance of the simulation risk estimates for sediment exposures.     

Integrated Uncertainty Analysis for Ambient Pollutant Health Risk Assessment: A Case Study of Ozone Mortality Risk

The U.S. Environmental Protection Agency (EPA) uses health risk assessment to help inform its decisions in setting national ambient air quality standards (NAAQS). EPA's standard approach is to make epidemiologically‐based risk estimates based on a single statistical model selected from the scientific literature, called the “core” model. The uncertainty presented for “core” risk estimates reflects only the statistical uncertainty associated with that one model's concentration‐response function parameter estimate(s). However, epidemiologically‐based risk estimates are also subject to “model uncertainty,” which is a lack of knowledge about which of many plausible model specifications and data sets best reflects the true relationship between health and ambient pollutant concentrations. In 2002, a National Academies of Sciences (NAS) committee recommended that model uncertainty be integrated into EPA's standard risk analysis approach. This article discusses how model uncertainty can be taken into account with an integrated uncertainty analysis (IUA) of health risk estimates. It provides an illustrative numerical example based on risk of premature death from respiratory mortality due to long‐term exposures to ambient ozone, which is a health risk considered in the 2015 ozone NAAQS decision. This example demonstrates that use of IUA to quantitatively incorporate key model uncertainties into risk estimates produces a substantially altered understanding of the potential public health gain of a NAAQS policy decision, and that IUA can also produce more helpful insights to guide that decision, such as evidence of decreasing incremental health gains from progressive tightening of a NAAQS.     

Quantitative Assessment of the Risk of Lung Cancer Associated with Occupational Exposure to Refractory Ceramic Fibers

We present the results of a quantitative assessment of the lung cancer risk associated with occupational exposure to refractory ceramic fibers (RCF). The primary sources of data for our risk assessment were two long-term oncogenicity studies in male Fischer rats conducted to assess the potential pathogenic effects associated with prolonged inhalation of RCF. An interesting feature of the data was the availability of the temporal profile of fiber burden in the lungs of experimental animals. Because of this information, we were able to conduct both exposure–response and dose–response analyses. Our risk assessment was conducted within the framework of a biologically based model for carcinogenesis, the two-stage clonal expansion model, which allows for the explicit incorporation of the concepts of initiation and promotion in the analyses. We found that a model positing that RCF was an initiator had the highest likelihood. We proposed an approach based on biological considerations for the extrapolation of risk to humans. This approach requires estimation of human lung burdens for specific exposure scenarios, which we did by using an extension of a model due to Yu. Our approach acknowledges that the risk associated with exposure to RCF depends on exposure to other lung carcinogens. We present estimates of risk in two populations: (1) a population of nonsmokers and (2) an occupational cohort of steelworkers not exposed to coke oven emissions, a mixed population that includes both smokers and nonsmokers.