Next Article in Journal
Application of Long-Short-Term-Memory Recurrent Neural Networks to Forecast Wind Speed
Previous Article in Journal
A Disc-Type High Speed Rotary Ultrasonic Motor with Internal Contact Teeth
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Risk Assessment for Ozone Regulation Based on Statistical Rollback

1
Department of Statistics, Kyungpook National University, Daegu 41566, Korea
2
Department of Statistics, Colorado State University, Fort Collins, CO 80523, USA
*
Author to whom correspondence should be addressed.
Submission received: 13 February 2021 / Revised: 27 February 2021 / Accepted: 4 March 2021 / Published: 8 March 2021
(This article belongs to the Section Environmental Sciences)

Abstract

:
In environmental studies, it is important to assess how regulatory standards for air pollutants affect public health. High ozone levels contribute to harmful air pollutants. The EPA regulates ozone levels by setting ozone standards to protect public health. It is thus crucial to assess how various regulatory ozone standards affect non-accidental mortality related to respiratory deaths during the ozone season. The original rollback approach provides an adjusted ozone process under a new regulation scenario in a deterministic fashion. Herein, we consider a statistical rollback approach to allow for uncertainty in the rollback procedure by adopting the quantile matching method so that it provides flexible rollback sets. Hierarchical Bayesian models are used to predict the potential effects of different ozone standards on human health. We apply the method to epidemiologic data.

1. Introduction

Regulating high ozone levels is essential as exposure to ozone can affect the risk of respiratory diseases or related deaths. Tropospheric ozone, also known as ground-level ozone, is one of the main harmful air pollutants that can cause adverse health effects. Many areas in the United States have been observed to exceed the current ozone National Ambient Air Quality Standards (NAAQS). Elevated ozone concentrations have also been a growing concern for rapidly developing nations where emissions of ozone precursors have been risen from expanding transportation networks [1]. Recently, U.S. EPA promulgated a new ozone NAAQS as 0.070 ppm (parts per million) and this change has prompted studies on the effects of the new ozone standard [2,3,4]. To investigate how different ozone NAAQS affect public health and welfare can be of great interest in epidemiologic studies, toxicological studies, or controlled human exposure studies. Based on the reviews of the air quality criteria for ozone ( O 3 ) and related photochemical oxidants and the NAAQS for O 3 , a modification in the current ozone regulatory standards provides the required protection for public health and welfare (see also [5,6]). Bell et al. [7,8] disapproved of the EPA’s scientific reviews reporting that designations based on air quality data from 2006 to 2008 would be effective in air quality data to take effect in 2010 for the 2007 8-h ozone standard. Therefore, it would be interesting to investigate whether the current regulations are sufficiently stringent to prevent respiratory-related mortality or not.
To assess how changes in the ozone regulations affect mortality, ground-level ozone must be adjusted by strengthening the air quality standards. Although rollback functions, namely, air quality adjustment procedures proposed by the EPA [9], can be useful for adjusting the ozone process, it cannot introduce sufficient variability in the rollback adjustment as the adjustment is deterministic and the EPA regulatory standards are based on the average of three consecutive years’ AQI values. Thus, we consider a parametric rollback approach with quantile matching method to allow for uncertainty in the rollback procedure
To conduct the risk assessment, we consider hierarchical Bayesian models that provide uncertainty quantification for relevant parameters to predict the potential effects of different regulatory ozone standards. We also describe some variations in the results under different modeling assumptions. We analyze databases from the National Morbidity, Mortality, and Air Pollution Study [10] that is designed to study the public health effects of air pollutants. The Health Effects Institutes began this study with researchers from Johns Hopkins and Harvard University in 1996. The NMMAPS database contains 108 U.S. urban areas from 1987 to 2000 where one can build a multisite time-series model of ozone and mortality simultaneously with meteorology information such as temperature or dew point and air pollution ( O 3 , PM 10 , SO 2 , NO 2 , and CO ) (see also [11]).
The focus of this paper is on statistical approaches to describe how the new ozone regulations affect mortality. We shall describe our results that mortality decreases as limits of acceptable ozone level get lower through the statistical rollback approach.

2. Statistical Rollback Approach

The rollback transformation [9] is one of the methods to adjust current ozone processes to follow new ozone regulatory standards. There are two main conditions to consider in air quality data. The first one is baseline conditions characterized by unadjusted air quality data monitored at fixed locations during recent years. The other one is attainment conditions generated by fitting an air quality adjustment procedure (AQAP) to the baseline. QQ-plots of high ozone levels in 1992 under baseline conditions against low ozone levels in 1999 under attainment conditions are shown in Figure 1. The figure exhibits different behaviors at high quantiles.

2.1. Quantile Matching Approach

Suppose that current ozone level, denoted by x, is a random variable having a distribution function F θ . We assume that adjusted ozone level corresponding to x under a new scenario follows the same distribution with different parameters, F θ * , and there exists a mapping m such that a random variable m ( x ) has the distribution function F θ * . Thus, the adjusted ozone level can be obtained using m ( x ) . That is, θ * is the adjusted model parameter based on the relationship between current levels and adjusted levels under the new scenario. Moreover, let F ( x ; θ ) and G ( z ; γ ) be the ozone distribution functions based on the baseline condition and attainment condition, respectively, (e.g., G F ). That is, m ( x ) = F θ * 1 ( F θ ( x ) ) or m ( x ) = G 1 ( F ( x ) ) . Note that F and G can be estimated by baseline and attainment years’ ozone data, respectively. Attainment year based on its AQI values or the condition G F is required. G may require adjustment based on ρ ( s ) = x ( b ) x ( s ) , where x ( b ) is current design values and x ( s ) is design values under scenario s.
Let f θ be a density function of current ozone levels x. Then adjusted rollback values z corresponding to x are obtained as follows:
  • Based on current ozone levels x = { x i } i = 1 n , estimate parameters θ (i.e., θ ^ ).
  • Compute q i x , the quantile of x i such that
    q i x = x i f θ ^ ( z ) d z ,
    for i = 1 , , n .
  • Estimate parameters θ * under new scenario (i.e., θ ^ * ).
  • Determine the corresponding z i satisfying
    q i x = z i f θ ^ * ( z ) d z ,
    for i = 1 , , n .
  • { z i } i = 1 n are adjusted (rollback) values of x .
We extend this approach through the Bayesian framework by putting prior distributions of parameters, π ( θ ) . By generating samples from posterior distribution, π ( θ | x ) f θ π ( θ ) , we can generate rollback ensembles in the same fashion.
In general, the distribution of x can be expressed as a mixture of the form
f ( x ) = ( 1 w ) h ( x | θ 1 ) x < u w g ( x | θ 2 ) x u ,
where u is a threshold, h can be any parametric distribution such as Weibull, truncated Gamma or Gaussian distribution, and g is a Generalized Pareto distribution (GPD) with parameters θ 2 = ( σ , ξ ) .
The quantile matching approach is simple and straightforward but causes various problems in adjusting multiple ozone processes in large dimensions. Now we consider more feasible approaches in the next section.

2.2. Weibull Approach in Rollback

Let y t ( c ) denote hourly (or daily average, daily maximum, daily 8-h maximum) ozone level for each time t and city c and y t * ( c ) denote rollbacked process corresponding to baseline process y t ( c ) . That is,
y t * = γ ( y t ) η ,
where we omit city index for simplicity. We fit y t to a parametric distribution such as Weibull distribution.
We obtain following information for each year i and each city c:
( A Q I i ( c ) , κ i ( c ) , δ i ( c ) ) ,
where κ and δ are weibull parameters. Note that γ and η are functions of κ and δ .
Now we consider a model of the form
E κ i ( c ) δ i ( c ) = g 1 ( A Q I i ( c ) ; θ 1 ) g 2 ( A Q I i ( c ) ; θ 2 ) and C o v κ i ( c ) δ i ( c ) = Σ ( c ) .
For example, we consider a simple linear model:
E κ i ( c ) δ i ( c ) = 1 A Q I i ( c ) 0 0 0 0 1 A Q I i ( c ) α 1 ( c ) β 1 ( c ) α 2 ( c ) β 2 ( c )
and
C o v κ i ( c ) δ i ( c ) = Σ ( c ) = σ 11 ( c ) σ 12 ( c ) σ 21 ( c ) σ 22 ( c ) .
For new standard regulation A Q I * , we can obtain κ * ( c ) and δ * ( c ) for each city c using
κ * ( c ) = α 1 ( c ) + β 1 ( c ) A Q I * and δ * ( c ) = α 2 ( c ) + β 2 ( c ) A Q I * .
It is observed that Σ ( c ) introduces some uncertainty in κ * ( c ) and δ * ( c ) . In general, hierarchical structures for ( α 1 ( c ) , β 1 ( c ) , α 2 ( c ) , β 2 ( c ) ) can be considered. Usually, δ provides a better fit based on A C L V 1 . Thus, we may need to adopt the relationship between κ and δ :
δ = A C L V 1 log ( n ) 1 / κ ,
where
A C L V 1 = T A Q I × C L V 1 A Q I .
Otherwise, we may select an average value or baseline value of κ or δ .
Let P A Q I s ( c ) be the 4th largest 8-h daily maximum concentration of z t , which is
δ * x t δ κ / κ * .
Then,
y t * = T A Q I a ( c ) P A Q I s ( c ) z t .

2.3. Log-Normal Approach in Rollback

Let y t be a baseline process and y t * be rollbacked process. Assume that
y t * = γ ( y t ) η ,
where y t l o g N o r m a l ( μ , σ 2 ) . That is,
log y t N ( μ , σ 2 ) and log y t * N ( λ + η μ , η 2 σ 2 ) ,
where λ = log γ .
Similarily, we fit y t to a log-normal distribution for each time t and city c. It is observed that y t can be any covariate in the risk model (e.g., daily average, daily maximum, and daily 8-h maximum). Then, can obtain following information for each year i and each city c:
( A Q I i ( c ) , μ i ( c ) , σ i ( c ) ) .
Now we consider a simple linear model:
E μ i ( c ) σ i ( c ) = 1 A Q I i ( c ) 0 0 0 0 1 A Q I i ( c ) α 3 ( c ) β 3 ( c ) α 4 ( c ) β 4 ( c )
and
C o v μ i ( c ) σ i ( c ) = Σ ( c ) = σ 33 ( c ) σ 34 ( c ) σ 43 ( c ) σ 44 ( c ) .
For new standard regulation A Q I * , we can determine μ * ( c ) and σ * ( c ) for each city c using
λ + η μ = α 3 ( c ) + β 3 ( c ) A Q I * and η σ = α 4 ( c ) + β 4 ( c ) A Q I * ,
where μ and σ are estimated based on baseline year. In general, the following relationship is considered:
log T A Q I A Q I + μ + σ Φ 1 1 n 1 = λ + η μ + η σ Φ 1 1 n 1 .
It is observed that η = κ / κ * and λ = log γ = log ( δ * ) κ κ * log ( δ ) = log ( δ * ) η log ( δ ) . Furthermore, well-fitted models of μ ( c ) and δ ( c ) can provide a good estimation based on following relationships
λ + η μ = α 3 ( c ) + β 3 ( c ) A Q I * and λ = log ( δ * ) η log ( δ ) .
That is,
η ^ = α 3 ( c ) + β 3 ( c ) A Q I * log ( δ * ) μ log ( δ )
and
λ ^ = log ( δ * ) η ^ log ( δ ) .
We can estimate one parameter (e.g., δ * or μ * ) first and then tune other parameter (e.g., κ * or σ * ) to make rollbacked AQI set to be new standard regulation A Q I * based on the range of the parameter. In the 8-h daily maximum rollback, we can directly solve the following equation by plugging in estimated parameters of δ * or μ * :
A Q I * = γ ( A Q I ) η ,
where η = κ / κ * and γ = δ * / δ κ κ * . That is,
κ ^ * = κ log A Q I δ log A Q I * δ * .
Furthermore, the following empirical relationships can be useful.
daily 1 h maximum = 2.5 × daily average
and
daily 8 h maximum = 1.33 × daily average .
As previous approaches (except for the 8-h maximum covariate) do not guarantee rollbacked values satisfying exact new standard regulation A Q I * , hourly scale rollback will be more appropriate for adjustment after rollback. Notice that EPA standard regulation is based on the average of three consecutive years’ AQI values.

3. Application to NMMAPS Data

3.1. Statistical Modeling

We apply an overdispersed Poisson model in generalized linear models (see [8,12]) to the NMMAPS data. Denote Y t c by the number of daily non-accidental deaths in community c on day t. The Poisson process with intensity function μ t c can be expressed as
Y t c P o i s s o n ( μ t c ) with V a r ( Y t c ) = ϕ c μ t c ,
where the parameter ϕ c describes overdispersion for community c. Note that all overdispersed Poisson models are assumed to be mutually independent over time. We can also model the intensity function μ t c with some essential covariates such as ozone levels at different lags, seasonality, long-term trends, weather, and co-pollutants for three age groups (<65, 65–74, and ≥75 year). Natural cubic splines are useful tools for getting smoothing functions of time to account for seasonality and long-term trends in which influenza epidemics, for example, can affect mortality. The interaction term between smoothing functions of time and age-specific indicators (<65, 65–74, and ≥75 years) is considered as it can adjust the possible seasonal mortality patterns by age group. We also control for some potential covariates related to weather by smoothing functions of dew point, average dew points of the previous three days, temperature, the average temperature of the previous three days as follows:
log μ t c = β c x t c + α c D O W t + γ 1 c n s ( t i m e , 7 / y e a r ) + γ 2 c n s ( T t c , 6 )
+ γ 3 c n s ( T t 1 , t 3 c , 6 ) + γ 4 c n s ( D t c , 3 ) + γ 5 c n s ( D t 1 , t 3 c , 3 )
+ interaction terms for age groups and time ,
where μ t c is the expected number of deaths, x t c is the average daily O 3 concentrations of the current and the previous day in community c on the day t, D O W t is the days of the week (categorical) on day t. We define n s ( t i m e , 7 / y e a r ) , n s ( T t c , 6 ) , n s ( T t 1 , t 3 c , 6 ) , n s ( D t c , 3 ) and n s ( D t 1 , t 3 c , 3 ) as the natural cubic spline function of time with 7 degrees of freedom per year, temperature with 6 degrees of freedom, the average temperature of the previous three days with 6 degrees of freedom, dew point with 3 degrees of freedom and the average dew points of the previous three days with 3 degree of freedom, respectively. The last term in (4) indicates the interaction terms of age-specific indicators and natural cubic spline functions of time.
Relative mortality rates associated with exposure to ozone over the past few days can be estimated in a specific community by constrained or unconstrained distributed-lag models as daily ozone levels are readily available. Those models can be more flexible in the sense that they can be more suitable for exploring the time lag between exposures to ozone and deaths than single-lag models (see [7]).
For example, the constrained distributed-lag models (CDL) and unconstrained distributed-lag models (UDL) can be expressed as
β c x t c = β 0 c x t + β 1 c x ¯ t : t 3 c + β 2 c x ¯ t : t 6 c or β c x t c = j = 0 6 β j c x t j c ,
respectively.
The national average relative mortality rate caused by ozone effects can be estimated through a Bayesian hierarchical model in which both variability within-community and across-community can be accounted. The national average relative rate will be estimated by integrating estimates for the relative rate from the distributed-lag models for each specific community. Through the two-stage model, first, variation across communities is considered over the short-term ozone effects. The national average relative rate is then estimated. The ozone effects O 3 on mortality is modeled for each community c as follows:
β ^ c | β c , Σ ^ c M V N ( β c , Σ ^ c ) ,
where β c is the true relative rate in a specific community, β ^ c is its estimate, and Σ ^ c is the estimated covariance matrix corresponding to β ^ c . We put the multivariate normal distribution on β c ,
β c | μ , Ω M V N ( μ , Ω ) ,
where μ is the true national average relative rate and Ω is the true covariance matrix of relative rate in a specific community, β c . Note that there are still several sensitivity issues related to the modeling: (1) co-pollutant such as PM 10 can be included as a potential confounder; (2) days with high temperatures can be excluded to control for the effects of heat waves; (3) the degrees of freedom (df) in the smooth functions of time needs to be specified to control for long-term trends and seasonality; and (4) various ozone exposure metrics such as daily average, 1-h maximum, and 8-h maximum can be considered (see [7]).

3.2. Inferences

Estimating the impact of new ozone regulatory standards on the total nonaccidental deaths is not achieved through the national average relative rate μ as the rollback adjustment is applied to each community separately. Instead, we directly calculate the expected total nonaccidental deaths of the original and adjusted (rollbacked) ozone process, respectively.
Denote g c by a rollback function for community c, which adjusts observed ozone concentrations to meet new ozone regulation standards. Our focus is then on the difference in expected nonaccidental deaths between before and after rollback transformation. Let E [ log ( μ t c ) ] = x t c β ^ c + M t c θ ^ c , where M t c is the design matrix of covariates except for ozone (i.e., D O W t , n s ( t i m e , 7 / y e a r ) , n s ( T t c , 6 ) , and n s ( D t c , 6 ) , etc.), and x t c is the design matrix of unconstrained (or constrained) distributed-lag ozones for community c (see [7]). If the amount of ozone reduction r is the same for each community c, then the reduction ratio of the expected total death is 1 exp r β ^ c for each community c at time t. Thus, the reduction rate of the national expected total death is 1 exp r μ ^ at time t. As the μ is obscure, the reduction rate of the national expected total death can be better estimated by weighting mean effects across the 98 cities, where weights are proportional to populations.
For simplicity, the unconstrained distributed-lag model is only considered here although the constrained distributed-lag models can be also readily applied via reparameterization. Let z t be rollbacked ozone concentrations of x t for community c (i.e., z t c = g c ( x t c ) ). After the rollback transformation, the expected total deaths for community c during the ozone season is as follows
t T O 3 exp j = 0 6 β j g c ( x t j c ) + M t c θ ^ c ,
where T O 3 is the set of time indices for the ozone season. The log ratio of national expected death before and after the rollback transformation during ozone season can be expressed as
log c C t T O 3 w c , t exp j = 0 6 β j ( z t j c x t j c ) ,
where C is the set of indices of communities in NMMAPS data and w c , t is the weight of community c at time t,
w c , t = exp j = 0 6 β j x t j c + M t c θ ^ c c C t T O 3 exp j = 0 6 β j x t j c + M t c θ ^ c .
We can estimate the posterior distribution of the log ratio in (9) through the MCMC approach. As π ( β | Y ) can be approximated by π ( β | β ^ ) and a prior distribution of ( μ , Ω ) well, we can generate β from the posterior distribution π ( β | Y ) . One of the efficient algorithms in application is TLNISE (two level Normal independent sampling estimation) algorithm (see [13]). Without a Bayesian approach, restricted maximum likelihood (REML) method can provide similar estimates.
Finally, with samples from π ( β | Y ) the posterior distribution of the log ratio of national expected death in (9) can be estimated as follows
log c C t T O 3 exp j = 0 6 β j ( i ) z t j c + M t c θ ^ c c C t T O 3 exp j = 0 6 β j ( i ) x t j c + M t c θ ^ c ,
where β ( i ) indicates each sample from π ( β | Y ) or π ( β | β ^ ) for i = 1 , , N .

4. Results

We investigate how the rollback transformation predicts that ozone series change under new ozone regulation standards possibly affecting public health. The transformation is based on the current AQI and new regulation standards in the attainment years. Figure 2 shows Q-Q plots of ozone levles under new regulation standards against the current AQI using various rallback functions and statistical rollbacks. A common AQI is applied to all cities (common rollback) or a different AQI is applied to each city (city-specific rollback).
The regular meteorological model in the NMMAPS is based on nonlinear functions mainly consisting of temperature and dewpoint at lag 0 and the average of them at lags 13. However, meteorological confounding in a distributed-lag model for ozone may exist at lags 46 so that we account for temperature and dewpoint at lag 0 through 6 in the “distributed-lag” meteorological model through nonlinear splines with 4 and 3 df, respectively.
We use the reduction rate of the expected total death to assess the effect of the proposed new ozone regulatory. Posterior means and 95% credible intervals for total mortality per 1000 deaths are shown in Table 1 and Table 2. They are based on unconstrained distributed-lag models with a common rollback and city-specific rollback approach, respectively. Most credible intervals have positive lower limits in the common rollback approach except for Weibull rollback at level 75 regularization, providing good evidence that there are reduced rates. As expected, the mortality rate decreases as lower regulation increases. Similar results are observed in CDL. The ranges of the estimated reduction rates are approximately 0.94–1.9 for 75 ppb, 1.47–2.9 for 70 ppb and 4.9–5.0 for 60 ppb for all models and rollback functions. Statistical rollbacks tend to show similar results as quadratic rollbacks for overall regularization levels. The ranges of the estimated reduction rates are roughly 1.0–2.0 for 75 ppb, 1.4–2.8 for 70 ppb and 2.4–4.8 for 60 ppb for all rollback functions in the city-specific rollback approach. Posterior variances of the relative risks for each specific city vary with rollback functions. Rather than we include daily average ozone concentration as a covariate, we include daily maximum ozone concentration and daily 8-h maximum ozone concentration so that they reduce the mortality rate more in Table 3. We also conduct statistical inference with simple MLE and pooled MLE and then compare them with Bayesian inference. Both MLE and pooled MLE give a slightly higher reduction rate than one based on the Bayesian inference in Table 4. Several issues associated with high temperature on the mortality reduction. These are investigated by presenting mortality reduction with and without high temperaturein Table 5.
Several issues associated with high temperature on the mortality reduction. These are investigated by presenting mortality reduction with and without high temperature in Table 5.

5. Concluding Remarks

As the rollback approach provides an adjusted ozone process in a deterministic manner, we introduce a statistical rollback approach based on the quantile matching method to allow for uncertainty and thus provide flexible rollback sets. The proposed methods are applied to epidemiologic data (NMMAPs) and we assess the impact of new ozone regulation standards on public health under various settings. Another possible model is to consider different time lags in the Poisson process model with constrained or unconstrained distributed lags together. Using respiratory deaths rather than nonaccidental deaths may also provide useful results. One of the issues that need to be addressed is collinearity between ozone and temperature, for example. Other studies revealed strong effects in the Northeast and Industrial Midwest, less strong but still significant effects in the Southeast and possibly Southern California. Thus, we can also consider a regional or spatial structure in the method. Finally, the parametric rollback approach using quantile matching may allow uncertainty in ozone process itself, which can be extended to a fully Bayesian framework.

Author Contributions

Y.K. concieved the idea and developed the method presented herein and wrote the paper. J.L. performed data analysis. Both authors have read and approved the final manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (No. 2018R1D1A1B07043352).

Acknowledgments

We are grateful to the Editor-in-chief, Associate Editor and the anonymous referees for their helpful comments.

Conflicts of Interest

There are no potential conflict of interest.

References

  1. US Environmental Protection Agency. National Ambient Air Quality Standards for ozone, final rule. Fed. Regist. 1997, 62, 38855–38896. [Google Scholar]
  2. Zhang, J.; Wei, Y.; Fang, Z. Ozone pollution: A major health hazard worldwide. Front. Immunol. 2019, 10, 2518. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Samet, J.M. Ozone and respiratory health. The story continues. Am. J. Respir. Crit. Care Med. 2015, 192, 272–273. [Google Scholar] [CrossRef] [PubMed]
  4. Dockery, D.W.; Pope, C.A.; Xu, X.; Spengler, J.D.; Ware, J.H.; Fay, M.E.; Ferris, B.G.; Speizer, F.E. An association between air pollution and mortality in six U.S. cities. N. Engl. J. Med. 1993, 329, 1753–1759. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Ostro, B.; Chestnut, L. Assessing the health benefits of reducing particulate matter air pollution in the United States. Environ. Res. 1998, 76, 94–106. [Google Scholar] [CrossRef] [PubMed]
  6. Levy, J.I.; Carrothers, T.J.; Tuomisto, J.T.; Hammitt, J.K.; Evans, J.S. Assessing the public health benefits of reduced ozone concentrations. Environ. Health Perspect. 2001, 109, 1215–1226. [Google Scholar] [CrossRef] [PubMed]
  7. Bell, M.L.; McDermott, A.; Zeger, S.L.; Samet, J.M.; Dominici, F. Ozone and short-term mortality in 95 US urban communities. J. Am. Med. Assoc. 2004, 292, 2372–2378. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Bell, M.L.; Peng, R.D.; Dominici, F. The Exposure-Response Curve for Ozone and Risk of Mortality and the Adequacy of Current Ozone Regulations. Environ. Health Perspect. 2006, 114, 532–536. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Johnson, T. A Guide to Selected Algorithms, Distributions, and Databases Used in Exposure Models Developed by the Office Air Quality Planning and Standards; US Environmental Protection Agency: Washington, DC, USA, 2002. [Google Scholar]
  10. National Morbidity Mortality and Air Pollution Study (NMMAPS) Database. Available online: http://www.ihapss.jhsph.edu/ (accessed on 28 October 2004).
  11. Samet, J.M.; Zeger, S.L.; Dominici, F.; Curriero, F.; Coursac, I.; Dockery, D.W.; Schwartz, J.; Zanobetti, A. The National Morbidity, Mortality, and Air Pollution Study. Part II: Morbidity and mortality from air pollution in the United States. Res. Rep. Health Eff. Inst. 2000, 94, 5–70. [Google Scholar] [PubMed]
  12. McCullagh, P.; Nelder, J.A. Generalized Linear Models; Chapman & Hall/CRC: London, UK, 1989. [Google Scholar]
  13. Everson, P.J.; Morris, C.N. Inference for multivariate normal hierarchical models. J. R. Stat. Soc. Ser. B 2000, 62, 399–412. [Google Scholar] [CrossRef]
Figure 1. QQ-plots of high ozone levels in 1992 under baseline conditions against low ozone levels in 1999 under attainment conditions at six different locations.
Figure 1. QQ-plots of high ozone levels in 1992 under baseline conditions against low ozone levels in 1999 under attainment conditions at six different locations.
Applsci 11 02388 g001
Figure 2. Q-Q plots based on various rallback functions (left) and statistical rollbacks (right) based on ozone data in New York (top) and Detroit (bottom) during 2003.
Figure 2. Q-Q plots based on various rallback functions (left) and statistical rollbacks (right) based on ozone data in New York (top) and Detroit (bottom) during 2003.
Applsci 11 02388 g002
Table 1. Posterior means and 95% confidence intervals of total mortality reduction per 1000 deaths (1998–2000): common rollback.
Table 1. Posterior means and 95% confidence intervals of total mortality reduction per 1000 deaths (1998–2000): common rollback.
ModelReg.Prop.Prop. w/ BGQuadraticWeibullStat.WStat.LN
level 751.871.090.95–0.620.980.94
(1.07, 2.62)(0.62, 1.53)(0.54, 1.35)(–1.18, –0.11)(0.52, 1.37)(0.49, 1.33)
CDLlevel 702.861.661.460.601.511.47
(1.65, 4.09)(0.95, 2.39)(0.84, 2.08)(0.10, 1.09)(0.88, 2.14)(0.83, 2.10)
level 604.922.862.512.982.612.54
(2.84, 7.02)(1.64, 4.09)(1.43, 3.61)(1.53, 4.30)(1.48, 3.75)(1.46, 3.65)
level 751.911.110.98–0.621.030.98
(1.15, 2.72)(0.66, 1.59)(0.58, 1.39)(–1.17, –0.07)(0.63, 1.44)(0.59, 1.40)
UDLlevel 702.921.701.490.591.581.52
(1.66, 4.14)(0.97, 2.41)(0.84, 2.10)(0.10, 1.09)(0.92, 2.18)(0.85, 2.13)
level 605.012.922.563.052.792.60
(2.93, 7.19)(1.66, 4.23)(1.49, 3.70)(1.73, 4.38)(1.72, 3.93)(1.56, 3.78)
Table 2. Posterior means and 95% confidence intervals of total mortality reduction per 1000 deaths (1998–2000): city-specific rollback.
Table 2. Posterior means and 95% confidence intervals of total mortality reduction per 1000 deaths (1998–2000): city-specific rollback.
ModelReg.Prop.Prop. w/ BGQuadraticWeibullStat.WStat.LN
level 751.961.200.98–0.521.020.99
(1.06, 2.92)(0.61, 1.79)(0.51, 1.46)(–1.33, 0.25)(0.54, 1.51)(0.50, 1.48)
CDLlevel 702.751.661.380.441.421.39
(1.47, 4.02)(0.87, 2.46)(0.74, 2.03)(–0.36, 1.26)(0.78, 2.09)(0.75, 2.05)
level 604.622.762.402.672.562.49
(2.64, 6.69)(1.58, 4.06)(1.36, 3.52)(1.38, 3.99)(1.52, 3.70)(1.45, 3.62)
level 751.961.190.98–0.511.010.99
(0.93, 2.91)(0.52, 1.79)(0.45, 1.46)(–1.36, 0.36)(0.49, 1.49)(0.50, 1.48)
UDLlevel 702.781.691.400.451.441.41
(1.50, 3.96)(0.88, 2.42)(0.74, 1.99)(–0.42, 1.23)(0.79, 2.08)(0.74, 2.03)
level 604.812.882.502.792.622.58
(2.66, 7.05)(1.55, 4.26)(1.36, 3.70)(1.41, 4.31)(1.49, 3.84)(1.44, 3.79)
Table 3. Posterior means and 95% confidence intervals of total mortality reduction per 1000 deaths (1998–2000): city-specific rollback and CDL model.
Table 3. Posterior means and 95% confidence intervals of total mortality reduction per 1000 deaths (1998–2000): city-specific rollback and CDL model.
CovariateReg.Prop.Prop. w/ BGQuadraticWeibullStat.WStat.LN
Dailylevel 751.961.200.98–0.521.020.99
Ave (1.06, 2.92)(0.61, 1.79)(0.51, 1.46)(–1.33, 0.25)(0.54, 1.51)(0.50, 1.48)
level 604.622.762.402.672.562.49
(2.64, 6.69)(1.58, 4.06)(1.36, 3.52)(1.38, 3.99)(1.52, 3.70)(1.45, 3.62)
Dailylevel 752.321.911.662.482.122.09
Max (1.54, 3.14)(1.26, 2.59)(1.09, 2.27)(1.62, 3.40)(1.59, 2.79)(1.55, 2.73)
level 606.195.104.255.134.874.59
(4.11, 8.31)(3.33, 6.87)(2.79, 5.75)(3.41, 6.91)(3.38, 6.35)(3.09, 6.06)
Dailylevel 752.101.611.361.901.561.52
8 h Max (1.26, 2.92)(0.96, 2.23)(0.79, 1.94)(1.13, 2.67)(0.96, 2.15)(0.94, 2.12)
level 605.534.233.434.343.673.55
(3.37, 7.50)(2.52, 5.76)(1.99, 4.69)(2.52, 5.91)(2.26, 4.94)(2.14, 4.81)
Table 4. Posterior means and 95% confidence intervals of total mortality reduction per 1000 deaths (1998–2000): city-specific rollback and CDL model.
Table 4. Posterior means and 95% confidence intervals of total mortality reduction per 1000 deaths (1998–2000): city-specific rollback and CDL model.
CovariateReg.Prop.Prop. w/ BGQuadraticWeibullStat.WStat.LN
Bayesianlevel 751.961.200.98–0.521.020.99
(1.06, 2.92)(0.61, 1.79)(0.51, 1.46)(–1.33, 0.25)(0.54, 1.51)(0.50, 1.48)
level 604.622.762.402.672.562.49
(2.64, 6.69)(1.58, 4.06)(1.36, 3.52)(1.38, 3.99)(1.52, 3.70)(1.45, 3.62)
MLElevel 752.231.311.09–0.301.211.19
(1.14, 3.21)(0.62, 1.91)(0.54, 1.58)(–1.15, 0.48)(0.66, 1.70)(0.64, 1.69)
level 605.072.932.523.032.562.52
(3.00, 7.28)(1.65, 4.29)(1.42, 3.68)(1.69, 4.37)(1.49, 3.73)(1.39, 3.69)
Pooledlevel 751.921.160.95–0.501.151.10
MLE (1.79, 2.05)(1.09, 1.24)(0.89, 1.02)(–0.61, –0.38)(1.08, 1.22)(1.04, 1.18)
level 604.652.782.402.692.492.44
(4.37, 4.92)(2.62, 2.93)(2.27, 2.54)(2.48, 2.88)(2.38, 2.62)(2.32, 2.59)
Table 5. Posterior means and 95% confidence intervals of total mortality reduction per 1000 Deaths (1998–2000): without high temperature (>85).
Table 5. Posterior means and 95% confidence intervals of total mortality reduction per 1000 Deaths (1998–2000): without high temperature (>85).
Prop.Prop. w/ BGQuadraticWeibullStat.WStat.LN
With High Temp1.961.200.98–0.521.020.99
(1.06, 2.92)(0.61, 1.79)(0.51, 1.46)(–1.33, 0.25)(0.54, 1.51)(0.50, 1.48)
Without High Temp1.390.830.67–0.470.750.71
(0.33, 2.41)(0.15, 1.46)(0.15, 1.17)(–1.36, 0.45)(0.22, 1.26)(0.18, 1.22)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kim, Y.; Lee, J. A Risk Assessment for Ozone Regulation Based on Statistical Rollback. Appl. Sci. 2021, 11, 2388. https://0-doi-org.brum.beds.ac.uk/10.3390/app11052388

AMA Style

Kim Y, Lee J. A Risk Assessment for Ozone Regulation Based on Statistical Rollback. Applied Sciences. 2021; 11(5):2388. https://0-doi-org.brum.beds.ac.uk/10.3390/app11052388

Chicago/Turabian Style

Kim, Yongku, and Jeongjin Lee. 2021. "A Risk Assessment for Ozone Regulation Based on Statistical Rollback" Applied Sciences 11, no. 5: 2388. https://0-doi-org.brum.beds.ac.uk/10.3390/app11052388

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop