 Published:
We’re sorry, something doesn't seem to be working properly.
Please try refreshing the page. If that doesn't work, please contact support so we can address the problem.
A tensor product quasiPoisson model for estimating health effects of multiple ambient pollutants on mortality
Environmental Healthvolume 18, Article number: 38 (2019)
Abstract
Background
People are exposed to mixtures of highly correlated gaseous, liquid and solid pollutants. However, in previous studies, the assessment of air pollution effects was mainly based on singlepollutant models or was simultaneously included as multiple pollutants in a model. It is essential to develop appropriate methods to accurately estimate the health effects of multiple pollutants in the presence of a high correlation between pollutants.
Methods
The flexible tensor product smooths of multiple pollutants was applied for the first time in a quasiPoisson model to estimate the health effects of SO_{2}, NO_{2} and PM_{10} on daily allcause deaths during 2005–2012 in Guangzhou, China. The results were compared with those from three other conventional models, including the singlepollutant model and the threepollutant model with and without firstorder interactions.
Results
The tensor product model revealed a complex interaction among three pollutants and significant combined effects of PM_{10}, NO_{2} and SO_{2}, which revealed a 2.53% (95%CI: 1.03–4.01%) increase in mortality associated with an interquartilerange (IQR) increase in the concentrations of all three pollutants. The combined effect estimated by the singlepollutant model was 5.63% (95% CI: 3.96–7.34%). Although the conventional threepollutant models produced combined effect estimates (2.20, 95%CI, 1.18–3.23%; 2.78, 95%CI: 1.35–4.23%) similar to those of the tensor product model, they distorted the estimates and inflated the variances of the estimates when attributing the combined health effects to individual pollutants.
Conclusions
The singlepollutant model or conventional multipollutant model may yield misleading results in the presence of collinearity. The tensor product quasiPoisson regression provides a novel approach to the assessment of the health impacts of multiple pollutants by flexibly fitting the interaction effects and avoiding the collinearity problem.
Background
An increasing number of studies have demonstrated the association between daily ambient air pollution levels and adverse health outcomes worldwide [1,2,3,4]. Some studies generally estimated the health effect of each pollutant separately via a singlepollutant model controlling for time trends and meteorological factors [5,6,7]. However, pollution occurs as a mixture of highly correlated gaseous, liquid and solid pollutants. Some studies observed that the effect using the singlepollutant model was weakened after other pollutants were simultaneously included [8]. Therefore, the results of the singlepollutant model tends to overestimate the health effects of pollutants. Moreover, many studies observed synergistic effects of combinations of pollutants on health outcomes [9, 10]. To advance our understanding of the biological mechanisms of pollutant toxicity, it is essential to move from the singlepollutant model to the multipollutant model [11]. This would provide important information for guiding regulatory policies for public health.
Many researchers have recognized the limitations of the singlepollutant model and have taken some efforts to overcome such limitations. One of commonly used methods is to simultaneously include two or three pollutant variables in a generalized additive model (GAM); then, the effects of one pollutant can be estimated after adjusting for other pollutants in the model [12,13,14]. This simple resolution ignores the high correlation of the occurrence of pollutants, which likely leads to a collinearity problem, such as variance inflation or bias even reverse estimation of parameters. For example, Chen et al. [15] and Kan et al. [16] showed that PM_{10} was negatively associated with daily respiratory mortality when the NO_{2} or SO_{2} was considered simultaneously, suggesting that the PM_{10} had a protective effect. To alleviate this problem, the Health Effects Institute supports the development of innovative statistical methods for studying the health effects of multiple pollutants in the presence of collinearity.
Several different approaches have been attempted to analyze the effects of multiple pollutants. One approach is the principal component analysis (PCA) that can avoid unstable parameter estimates in the condition of collinearity [17, 18]. However, PCA only maximizes the variance interpreted by the linear combination of independent variables without directly considering the relationship between the dependent variable and independent variables. Partial leastsquares regression (PLSR) combining PCA with canonical correlation analysis captures not only the information of independent variables but also the relationship with the dependent variable, yet this approach produces linear combinations of the original independent variables that are difficult to interpret in practice. Using Monte Carlo simulations, Månsson and Shukur [19] showed that the estimators produced by Poisson ridge regression were more accurate (i.e. had less standard errors) than those produced by the Poisson loglinear model. However, the increase in accuracy was at the expense of increasing bias in the condition of strong collinearity. Bayesian model averaging (BMA) averages the effects of models for all possible combinations of independent variables; these models are weighted by the posterior probability of the model [20]. However, there are different views on the interpretation of effect estimates from BMA in the case of extensive collinearity. Bobb et al. [21] proposed Bayesian kernel machine regression (BKMR) that regressed the heath outcome on a flexible kernel function of multiple pollutants. BKMR can be highly sensitive to the specification of the prior distribution and can solve the problem of collinearity using hierarchical variable selection.
To date, there is no standard model for evaluating health effects of multiple pollutants. Sun et al. [22] concluded that no method was unanimous across all simulation scenarios that differed in sample size, the number of pollutants and the strength of exposureresponse association. The choice of the appropriate method should depend on the goal of the study. Wood proposed a general method for constructing lowrank tensor product smooths of multiple variables that could quantitatively estimate the nonlinear effects of multiple independent variables on a dependent variable [23, 24]. The tensor product smooths has been widely applied in geometry, physics, mechanics and quantum theory [25]. To the best of our knowledge, the tensor product smooths for GAM was never employed to evaluate the health impacts of multiple pollutants. This study employed firstly this approach to assess nonlinear effects and potential interaction effects of SO_{2}, NO_{2} and PM_{10} on mortality in Guangzhou, China, during 2005–2012.
Methods
Data
Daily meteorological data for mean temperature, relative humidity and atmospheric pressure were collected at the Guangzhou Weather Station and were downloaded from the China Meteorological Data Sharing Service System. Air pollution data, including particulate matter < 10 μm in aerodynamic diameter (PM_{10}), nitrogen dioxide (NO_{2}) and sulfur dioxide (SO_{2}), were provided by the Guangzhou Bureau of Environmental Protection. The daily concentrations of PM_{10}, NO_{2} and SO_{2} were the average levels collected from seven fixedsite air quality monitoring stations that were all located in central urban areas and had complete data during the study period of 2005–2012 [26]. There were 40 missing observations for PM_{10}, NO_{2}, and SO_{2}, accounting for 1.5% of the study days. We did not replace the missing values. Because the levels of air pollution and health effects of pollution on mortality may be different in suburban areas, to maintain consistency of the study population, the mortality data were also limited to death registration records of the whole population at six central urban districts.
Statistical analysis
First, we constructed a basic quasiPoisson regression model without air pollutant variables. Natural cubic splines (ns) of daily mean relative humidity, atmospheric pressure and temperature were used to control the nonlinear and nonmonotonic confounding effects of meteorological factors on mortality. Time can be considered as a proxy of some other unobserved timevarying confounders. A natural cubic spline function of time was used to control for the longterm trend and seasonal variations of daily mortality and potential confounding effects of unobserved timevarying covariates. The partial autocorrelation function (PACF) of residuals and the Akaike’s Information Criterion (AIC) were used to guide the selection of the number degrees of freedom (df) for ns. When a df of 7 was used for time variables, the PACF was free of patterns and was not autocorrelated [16, 27]. The value of df equal to 3 was used for meteorological variables because this value has been shown to control well their impacts on daily mortality changes [28,29,30]. We used a 14day moving average of mean temperature to control the lag effect because previous studies in Guangzhou observed that the effects of mean temperature on mortality persisted for approximately 14 days [31, 32]. Day of the week (DOW) and a holiday (HOD) indicator were included in the model as indicator variables.
After the basic model was constructed, we added air pollution variable(s) to the model. The reason was that many studies, including our data, showed that the adverse impacts of each specific pollutant were statistically significant at lag 0 and lag 1, and these impacts almost disappeared at lag 3 and during longer lags [33, 34]. To adequately capture the lag effects, a 3day moving average of air pollution concentrations (i.e., lag 0–2) was used, and the concentrations were introduced into the model in the following four ways:
Model I (singlepollutant model) incorporated each of pollutant variables into the basic model, separately:
where β(PM_{10t}/NO_{2t}/SO_{2t}) signifies that respective singlepollutant models for PM_{10}, NO_{2} and SO_{2} are established; β is the coefficient of a pollutant variable.
Model II (threepollutant noninteraction model), simultaneously included three pollutant variables without considering firstorder interactions and is specified as follows:
Model III (threepollutant interaction model) adds linear terms of the firstorder interaction between pollutants to model II:
where β_{1}, β_{2} and β_{3} are the regression coefficients of pollutants, and β_{12}, β_{13} and β_{23} are the regression coefficients of firstorder interaction terms.
Model IV (tensor product model) constructs the tensor product smooths (te) of three pollutants (PM_{10}, NO_{2} and SO_{2})and is specified as follows:
where te(PM_{10}, NO_{2}, SO_{2}) is the tensor product smooth function of three pollutants. For models IIV, μ_{t} is the expected number of deaths on day t, α is the intercept, ns(.) denotes the unknown smooth functions modeled by natural cubic splines, and θ_{1} θ_{4} represent the coefficients vector of smooth functions.
Construction of the tensor product smooths has been described previously in detail [23]. Mathematically, supposing that there are two vector spaces V and W, where the V has a basis e_{i} (i = 1 … m), and the W has a basis f_{j} (j = 1 … n), the tensor product of V ⊗ W is generated by \( \sum \limits_{i=1}^m\sum \limits_{j=1}^n{e}_i\otimes {f}_j \), where the product operation ⊗ is Kronecker product. In this paper, we briefly introduce the construction process of tensor product smooths of three covariates (i.e., PM_{10}, NO_{2}, SO_{2}). It is easy to generalize the tensor product smooths to more variables in the same way. Basis functions for representing the smooth function of each pollutant are given by \( f\left({PM}_{10}\right)=\sum \limits_{i=1}^I{\alpha}_i{\phi}_i\left({\mathrm{PM}}_{10}\right) \), \( f\left({\mathrm{NO}}_2\right)=\sum \limits_{k=1}^K{\omega}_k{\eta}_k\left({\mathrm{NO}}_2\right) \) and \( f\left({\mathrm{SO}}_2\right)=\sum \limits_{l=1}^L{b}_l{\varphi}_l\left({\mathrm{SO}}_2\right) \),
where the ϕ_{i}(PM_{10}), φ_{l}(SO_{2}) and η_{k}(NO_{2}) are the basis functions and α_{i}, b_{l} and ω_{k} are parameters. Next, the smooth of PM_{10}, denoted by f(PM_{10}), is converted into a smooth function of PM_{10} and SO_{2}, denoted by f(PM_{10}, SO_{2}) (to make f(PM_{10}) change smoothly with SO_{2}), which can be achieved by allowing α_{i} to vary smoothly with SO_{2}. Thus, we obtain
which results in
Similarly, we can now create a smooth function of PM_{10}, SO_{2} and NO_{2} by allowing f(PM_{10}, SO_{2}) to vary smoothly with NO_{2}. Therefore, with following the same reasoning as before, we obtain
where v_{ilk} are the coefficients of basis functions of tensor product smooths.
According to the above description, it is clear that the complexity of the tensor product smooth function depends on the type of basis function and its dimensionality; therefore, it is essential to choose an appropriate basis function describing the relationship between an adverse outcome and air pollutants. In this study, we used the cubic spline basis function that is widely used to flexibly fit a potential nonlinear relationship.
To avoid the problem of collinearity and controlling the tradeoff between goodness of fit and model smoothness, we define the following penalized likelihood function to estimate parameters:
(and)
where l (v) is the likelihood function, and J(f(PM_{10}, NO_{2}, SO_{2})) represents tensor product penalties. \( {\lambda}_{{\mathrm{PM}}_{10}} \), \( {\lambda}_{{\mathrm{NO}}_2} \) and \( {\lambda}_{{\mathrm{SO}}_2} \) are smoothing parameters. Such parameters can be selected by Generalized Cross Validation (GCV). The estimates, denoted by \( {\widehat{v}}_{ilk} \), are obtained by penalized iteratively reweighted least squares which maximizes the penalized likelihood function. Furthermore, stratified analyses by educational attainment, sex and age were conducted respectively.
Estimation of the combined effects
The combined effects of three air pollutants were measured by the combined rate ratio (RR). For the singlepollutant model and the threepollutant noninteraction model, RR was calculated as the exponentiated sum of the products of each pollutant’s regression coefficient and the pollutant’s increment. For the threepollutant interaction model, RR was calculated as the exponentiated sum of the following two parts: 1) the product of each pollutant’s regression coefficient and the pollutant’s increment and 2) the product of each interaction term coefficient and the difference value of the interaction term. The standard error of RR was estimated based on the variancecovariance matrix of pollutants’ coefficients. The detailed derivation process is described in the supplemental materials and uses the method proposed by Winquistet al.[13]. For the tensor production model, the RR is calculated as follows:
where
Therefore, we can get RR = exp (te^{T}ν).
The values \( \boldsymbol{te}\left({{\mathrm{PM}}_{10}}^{{\mathrm{j}}^{\mathrm{th}}},{{\mathrm{NO}}_2}^{{\mathrm{j}}^{\mathrm{th}}},{{\mathrm{SO}}_2}^{{\mathrm{j}}^{\mathrm{th}}}\right) \) and \( \boldsymbol{te}\left({{\mathrm{PM}}_{10}}^{i^{\mathrm{th}}},{{\mathrm{NO}}_2}^{i^{\mathrm{th}}},{{\mathrm{SO}}_2}^{i^{\mathrm{th}}}\right) \) are the values of the tensor product smooths at the j^{th} and i^{th} percentiles of pollutants’ concentrations, respectively. The vector ν is the coefficient vector of the basis function of tensor product smooths. The vector te^{T} is the transposed difference vector of basis functions of tensor product smooths at the j^{th} percentiles and the i^{th} percentiles of pollutants’ concentrations. The 95% confidence interval (CI) of RR was estimated by the following formula:
where \( \widehat{\Sigma} \) is the variancecovariance of the coefficients of the basis function and α =0.05.
To improve the readability of the results and better understand the excess burden due to air pollution, we calculated the excess risk (ER) of mortality by (RR1)*100%, i.e., the percentage increase in mortality associated with air pollution.
To clarify the robustness of the lag selection for air pollutants, we conducted a sensitivity analysis to determine the effects of air pollution at lag 0–1 by using a 2day moving average in models IIV.
All analyses were performed using R software version 3.3.1, by using the “mgcv” used to fit the GAM.
Results
The total number of allcause deaths was 193,715 during 2005–2012, with 66.34 cases per day. The daily mean temperature had an average of 22.5 ^{°}C and ranged from 5.1 ^{°}C to 34.2 ^{°}C. The mean daily average concentrations of PM_{10,} SO_{2} and NO_{2} were 74.9 μg/m^{3}, 41.0 μg/m^{3} and 62.4 μg/m^{3}, respectively. The interquartile ranges (IQR) of three pollutants were 48.2 μg/m^{3}, 33.7 μg/m^{3} and 33.3 μg/m^{3}, respectively (Table 1). The time series plot of air ambient pollutants and allcause deaths revealed a consistent seasonality; there were generally more deaths and higher concentrations of pollutants in the cool season (from May to October) than in the warm season (from November to April) (Fig. 1). The Spearman correlations of PM_{10} and NO_{2}, PM_{10} and SO_{2}, and NO_{2} and SO_{2} were 0.83, 0.60 and 0.65, respectively (P < 0.01).
The tensor product model can produce a fourdimensional plot for the RR of mortality associated with the levels of three pollutants. To intuitively show the combined effects of air pollutants, the concentration of one of pollutants was fixed at the reference value (the 25th percentile was used in this study); then, a threedimensional graph of RR of the other two pollutants can be drawn (Fig. 2). Generally, as the concentration of pollutants increased, a stronger combined effect of pollutants on mortality was observed, but the relationship was not completely linear or monotonic. The effect of one pollutant varied with the concentration of another pollutant. For example, the effect of PM_{10} increased significantly with the concentration of SO_{2} at low concentrations, but the curve of PM_{10} effect was relatively flat for high concentrations of SO_{2}, which indicated that there was an interaction between pollutants.
Table 2 and Fig. 3 show the effects of air pollution estimated by four different models. The singlepollutant model (Model I) revealed that an IQR increase in the concentration of PM_{10}, NO_{2} and SO_{2} was associated with increases in mortality of 2.03% (95% CI: 1.08–2.99%), 1.63% (95% CI: 0.69–2.59%) and 1.86% (95% CI: 0.94–2.78%), respectively. The level of NO_{2} was negatively associated with the risk of mortality if PM_{10} and SO_{2} were included simultaneously in the noninteraction model (models IIIII). The combined effect of three pollutants estimated by the singlepollutant model (ER = 5.63, 95%CI: 3.96–7.34%) was much greater than those obtained from other three multipollutant models. The combined effect estimated by the threepollutant noninteraction model (ER = 2.20, 95%CI: 1.18–3.23%) was slightly smaller than those from the interaction model (ER = 2.78, 95%CI: 1.35–4.23%) and the tensor product model (ER = 2.53, 95%CI: 1.03–4.01%). Fig. 3(ac) shows that the tensor product model, compared to the threepollutant noninteraction model and interaction model, produced more precise effect estimates (i.e., narrower 95%CI), especially for the main effects of single pollutants.
Table 3 shows the results of stratified analyses by the tensor product model. There were a total number of 193,715 deaths during 2005–2012 in Guangzhou; 24.4, 56.4 and 92.6% of deaths were of those aged under 65, males and individuals with low educational attainment (illiterate or only primary education), respectively. Females, elderly people (65 years of age or above) and individuals with low educational attainment were more vulnerable to ambient air pollution exposure than were males, younger individuals and those with high educational attainment.
The sensitivity analysis for all models revealed that the effects at lag 0–1 were slightly larger than those at lag 0–2 (Additional file 1: Figure S1).
Discussion
Most of previous studies of the health effects of air pollution are based on a singlepollutant model or include multiple pollutants in a model at the same time, thus ignoring the possibility of high correlation among pollutants that leads to the inflation of variance and the instability of model parameter estimation. The commonly used solution is stepwise regression and removing one of the highly correlated variables; however, this method may cause the loss of some important information. In this study, we proposed a quasiPoisson regression with tensor product smooth to estimate the health effects of multiple pollutants. Overall, we observed that the combined effect of three pollutants (PM_{10}, SO_{2} and NO_{2}) gradually increased along with concentration of pollutants. Moreover, there was a complex interaction among the pollutants. The results of the stratified analysis showed that elderly individuals, females and people with low educational attainment were more vulnerable to air pollution.
The singlepollutant model, a common method of assessing health effects of air pollution, revealed that the levels of each single pollutant (PM_{10}, NO_{2} and SO_{2}) had a significantly positive relationship with the risk of allcause deaths, with ER values associated with IQR increase in the concentration of lag 0–2 of 2.03, 1.63 and 1.86%, respectively. We observed that the effect estimates of PM_{10} and NO_{2}were greater at lag 0–1, with ER values per IQR increase of 2.13, 1.73 and 1.79% (i.e., 0.44, 0.52 and 0.53% per 10ug/m^{3} increase). Consistently, Wu et al. reported that ER of mortality associated with IQR increase in PM_{10} in Guangzhou was 1.77% (95%CI: 0.37–3.18%) in 2006–2009 and 2.03% (95%CI: 0.79–3.30%) in 2010–2013, and the effect estimates were very similar after adjusting for NO_{2} or SO_{2} [35]. The estimate of the PM_{10} effect was also very similar to that reported in 90 American cities and in 38 Chinese cities (i.e., 0.5 and 0.44% per 10ug/m^{3} increase), while the effects of SO_{2} and NO_{2} we observed were smaller than the estimates in four Asian cities (i.e., 1.00 and 1.23% per 10μg/m^{3} increase) [34, 36].
The threepollutant models produced abnormal health impact estimates of single pollutants (e.g., a negative effect of NO_{2}) and w19pt?>The above phenomenon was due to strong collinearity among pollutants with Spearman correlation coefficients of 0.60–0.83, which caused variance inflation of the regression parameters and bias in the statistical inference [11]. As to the combined effects of three pollutants, it was clear that the additive effect of a single pollutant based on the singlepollutant model was much higher than the combined effect estimated by other three multipollutant models. The reasons are that the single pollutant model does not take into account the high correlation between pollutants, and one air pollutant may be a substitute for another air pollutant or the mixture of air pollutants [37].
The most innovative point of this study is that the tensor product quasiPoisson regression model was, for the first time, applied to assess the combined effect of several pollutants on adverse health outcome in the presence of multicollinearity and to explore the complex interaction of several pollutants. Several methods that were proposed in the past made additional assumptions of the existence of a linear or specific functional relationship between atmospheric pollutants and the response variable or only explored the first order interaction; in contrast, the GAM with tensor product smooth not only relaxes those assumptions but also flexibly explores the nonlinear effect and the complex interaction of multiple pollutants. Specially, this method gives a highly visual depiction of the relationship between allcause deaths and air pollutants, which may allow analysts to discover more undetectable characteristics of data than would be possible with the use of a simple model, and provide a new approach to comprehensive evaluation of health effects of multiple pollutants.
In this study, the tensor product model was also used to identify subpopulations that were more susceptible to exposure to a mixture of multiple pollutants. Some studies that assessed the health effect of single pollutants revealed that females were more susceptible than males to adverse effects of ambient air pollution [38]. Kan et al. observed that the effects of PM_{10} on mortality were almost twice as strong in females [39]. The reasons for sexual differences were unclear and deserved to be further investigated. We observed greater combined effects of multiple pollutants in females than in males. Residents with low educational attainment were more sensitive to exposure to ambient air pollution than those with high educational attainment, which was consistent with previous studies [39, 40]. Most individuals with low educational attainment experience higher financial risk and, therefore, potentially also experience inferior living conditions and inadequate healthcare. Additionally, individuals with low educational attainment were more likely to be occupationally exposed to air pollution [41].
Our study has some limitations. Although the GAM with tensor product smooth is flexible enough to explore the relationship between adverse health outcome and air pollutants, the “curse of dimensionality” occurs when the number of variables in the tensor product smooth increases, which is inherent in many flexible models [23]. The total dimensionality is the product of dimensionalities of each pollutant; thus, a large sample size is needed to accurately estimate parameters, and the computation is timeconsuming. Moreover, it is difficult to interpret results of the model in highdimensional spaces. We assessed health effects of three main pollutants to demonstrate the methods; however, other pollutants, such as ozone and PM_{2.5} that might be related to mortality were not considered in this study because of data unavailability during the study period.
Conclusions
Ambient air pollution has a significant impact on mortality. Our findings indicate that there exists a high correlation and a complex interaction among air pollutants. The singlepollutant model is unable to accurately estimate the combined effects of multiple pollutants. The simple multipollutant model with or without the firstorder interaction term may yield misleading results when attributing health effects to individual pollutants. This study provides a novel method, called a tensor product quasiPoisson model, for assessing the health effects of multiple pollutants, which can be applied in various epidemiologic studies to determine the combined effects of environmental, social and individual risk factors, especially if such factors are highly correlated and have nonlinear health effects.
Abbreviations
 BKMR:

Bayesian kernel machine regression
 BMA:

Bayesian model averaging
 CI:

Confidence interval
 df:

Degree of freedom
 ER:

Excess risk
 GAM:

Generalized additive model
 HOD:

Holiday
 IQR:

Interquartilerange
 NO_{2} :

Nitrogen dioxide
 ns:

Nature cubic spline
 PACF:

Partial autocorrelation function
 PCA:

Principal component analysis
 PLSR:

Partial leastsquare regression
 PM_{10} :

Particulate matter less than 10 μm in aerodynamic diameter
 RR:

Rate ratio
 SO_{2} :

Sulfur dioxide
 te:

Tensor product smooths
References
 1.
Pope CA, Burnett RT, Krewski D, Jerrett M, Shi Y, Calle EE, et al. Cardiovascular mortality and exposure to airborne fine particulate matter and cigarette smoke shape of the exposureresponse relationship. Circulation. 2009;120:941–8.
 2.
Zanobetti A, Gold DR, Stone PH, Suh HH, Schwartz J, Coull BA, et al. Reduction in heart rate variability with traffic and air pollution in patients with coronary artery disease. Environ Health Perspect. 2010;118:324–30.
 3.
Chen R, Kan H, Chen B, Huang W, Bai Z, Song G, et al. Association of particulate air pollution with daily mortality: the China air pollution and health effects study. Am J Epidemiol. 2012;175:1173–81.
 4.
Hart JE, Bertrand KA, DuPre N, James P, Vieira VM, VoPham T, et al. Exposure to hazardous air pollutants and risk of incident breast cancer in the nurses' health study II. Environ Health. 2018;17:28.
 5.
Dominici F, Burnett RT. Risk models for particulate air pollution. J Toxicol Environ Health A. 2003;66:1883–90.
 6.
Liu Y, Xie S, Yu Q, Huo X, Ming X, Wang J, et al. Shortterm effects of ambient air pollution on pediatric outpatient visits for respiratory diseases in Yichang city, China. Environ Pollut. 2017;227:116–24.
 7.
Dominici F, Mcdermott A, Zeger SL, Samet JM. On the use of generalized additive models in timeseries studies of air pollution and health. Am J Epidemiol. 2002;156:193–203.
 8.
Vedal S, Brauer M, White R, Petkau J. Air pollution and daily mortality in a city with low levels of pollution. Environ Health Perspect. 2003;111:45–52.
 9.
Jacoby DB. Synergism between rhinovirus infection and oxidant pollutant exposure enhances airway epithelial cell cytokine production. Environ Health Perspect. 2002;110:665–70.
 10.
Vempilly J, Abejie B, Diep V, Gushiken M, Rawat M, Tyner TR. The synergetic effect of ambient PM_{2.5} exposure and rhinovirus infection in airway dysfunction in asthma: a pilot observational study from the Central Valley of California. Exp Lung Res. 2013;39:434–40.
 11.
Francesca D, Peng RD, Barr CD, Bell ML. Protecting human health from air pollution: shifting from a singlepollutant to a multipollutant approach. Epidemiology. 2010;21:187–94.
 12.
Boezen HM, Sc VDZ, Postma DS, Vonk JM, Gerritsen J, Hoek G, et al. Effects of ambient air pollution on upper and lower respiratory symptoms and peak expiratory flow in children. Lancet. 1999;353:874–8.
 13.
Winquist A, Kirrane E, Klein M, et al. Joint effects of ambient air pollutants on pediatric asthma emergency department visits in Atlanta, 19982004. Epidemiology. 2014;25:666–73.
 14.
Tango T. Effect of air pollution on lung cancer: a Poisson regression model based on vital statistics. Environ Health Perspect. 1994;102:41–5.
 15.
Chen G, Song G, Jiang L, Zhang Y, Zhao N, Chen B, et al. Shortterm effects of ambient gaseous pollutants and particulate matter on daily mortality in Shanghai, China. J Occup Health. 2008;50:41–7.
 16.
Kan H, Chen B, Zhao N, London SJ, Song G, Chen G, et al. Part 1. A timeseries study of ambient air pollution and daily mortality in Shanghai, China. Res Rep Health Eff Inst. 2010;154:17–78.
 17.
Arif AA, Shah SM. Association between personal exposure to volatile organic compounds and asthma among US adult population. Int Arch Occup Environ Health. 2007;80:711–9.
 18.
Holberg CJ, O'Rourke MK, Lebowitz MD. Multivariate analysis of ambient environmental factors and respiratory effects. Int J Epidemiol. 1987;16:399–410.
 19.
Månsson K, Shukur G. A Poisson ridge regression estimator. Econ Model. 2011;28:1475–81.
 20.
Koop G, Tole L. Measuring the health effects of air pollution: to what extent can we really say that people are dying from bad air? J Environ Econ Manag. 2004;47:30–54.
 21.
Bobb JF, Valeri L, Claus HB, Christiani DC, Wright RO, Mazumdar M, et al. Bayesian kernel machine regression for estimating the health effects of multipollutant mixtures. Biostatistics. 2015;16:493–508.
 22.
Sun Z, Tao Y, Shi L, Ferguson KK, Meeker JD, Kyun PS, et al. Statistical strategies for constructing health risk models with multiple pollutants and their interactions: possible choices and comparisons. Environ Health. 2013;12:85.
 23.
Wood SN. Lowrank scaleinvariant tensor product smooths for generalized additive mixed models. Biometrics. 2006;62:1025–36.
 24.
Zhao Y, Zhang Q, Levesley J. Multilevel sparse grids collocation for linear partial differential equations, with tensor product smooth basis functions. Comput Math Appl. 2018;75:883–99.
 25.
Pulmannová S. Tensor products of quantum structures and their applications in quantum measurements. Int J Theor Phys. 2003;42:907–19.
 26.
Li L, Yang J, Guo C, Chen PY, Ou CQ, Guo Y. Particulate matter modifies the magnitude and time course of the nonlinear temperaturemortality association. Environ Pollut. 2015;196:423–30.
 27.
Yang J, Ou CQ, Song YF, Li L, Chen PY, Liu QY. Estimating years of life lost from cardiovascular mortality related to air pollution in Guangzhou, China. Sci Total Environ. 2016;573:1566–72.
 28.
Muggeo VM, Hajat S. Modelling the nonlinear multiplelag effects of ambient temperature on mortality in Santiago and Palermo: a constrained segmented distributed lag approach. Occup Environ Med. 2009;66:584–91.
 29.
Braga A, Zanobetti A, Schwartz J. The time course of weatherrelated deaths. Epidemiology. 2001;12:662–7.
 30.
Dominici F, Peng RD, Bell ML, Pham L, McDermott A, Zeger SL, et al. Fine particulate air pollution and hospital admission for cardiovascular and respiratory diseases. JAMA. 2006;295:1127–34.
 31.
Guo Y, Li S, Zhang Y, Armstrong B, Jaakkola JJ, Tong S, et al. Extremely cold and hot temperatures increase the risk of ischaemic heart disease mortality: epidemiological evidence from China. Heart. 2013;99:195–203.
 32.
Yang J, Liu HZ, Ou CQ, Lin GZ, Zhou Q, Shen GC, et al. Global climate change: impact of diurnal temperature range on mortality in Guangzhou, China. Environ Pollut. 2013;175:131–6.
 33.
Guo Y, Tong S, Li S, Barnett AG, Yu W, Zhang Y, et al. Gaseous air pollution and emergency hospital visits for hypertension in Beijing, China: a timestratified casecrossover study. Environ Health. 2010;9:57.
 34.
Samet JM, Zeger SL, Dominici F, Curriero F, Coursac I, Dockery DW, et al. 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 71–79.
 35.
Rongshan W, Liuju Z, Xiaoliang H, Hongbing X, Shuo L, Baihuan F, et al. Temporal variations in ambient particulate matter reduction associated shortterm mortality risks in Guangzhou, China: a timeseries analysis (2006–2016). Sci Total Environ. 2018;645:491–8.
 36.
Yin P, He G, Fan M, Chiu KY, Fan M, Liu C, et al. Particulate air pollution and mortality in 38 of China’s largest cities: time series analysis. BMJ. 2017;356:j667.
 37.
Roberts S. A new model for investigating the mortality effects of multiple air pollutants in air pollution mortality timeseries studies. J Toxicol Environ Health A. 2006;69:417–35.
 38.
Clougherty JE. A growing role for gender analysis in air pollution epidemiology. Cien Saude Colet. 2011;16:2221–38.
 39.
Kan H, London SJ, Chen G, Zhang Y, Song G, Zhao N, et al. Season, sex, age, and education as modifiers of the effects of outdoor air pollution on daily mortality in Shanghai, China: the public health and air pollution in Asia (PAPA) study. Environ Health Perspect. 2008;116:1183–8.
 40.
Li L, Yang J, Song YF, Chen PY, Ou CQ. The burden of COPD mortality due to ambient air pollution in Guangzhou, China. Sci Rep. 2016;6:25900.
 41.
Finkelstein MM, Jerrett M, DeLuca P, Finkelstein N, Verma DK, Chapman K, et al. Relation between income, air pollution and mortality: a cohort study. CMAJ. 2003;169:397–402.
Acknowledgements
We appreciate Guangzhou Centre for Disease Control and Prevention and Guangzhou Bureau of Environmental Protection for providing allcause death data and air pollution data, respectively.
Funding
This work was supported by National Nature Science Foundation of China [81573249] and Nature Science Foundation of Guangdong Province [2016A030313530]. The funders had no role in this study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
Raw data of allcause mortality data, daily air pollution concentrations and daily metrological data were requested and obtained from Guangzhou Centre for Disease Control, Prevention and Guangzhou Bureau of Environmental Protection and China Meteorological Data Sharing Service System, respectively. Raw data will not be shared because the authors are not authorized for distribution of data.
Author information
Author notes
Affiliations
Contributions
CQ initiated the study. LL and SQ collected and cleaned data. LJ and SQ performed the statistical analysis. LJ and CQ drafted the manuscript. All authors made contributions to interpretation of data and revising the manuscript. All authors read and approved the final manuscript.
Corresponding author
Correspondence to ChunQuan Ou.
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file
Additional file 1:
FigureS1. Excess risk of mortality and 95% confidence intervals (%) associated with an IQR increment in moving average lag of 0–1 and 0–2 of air pollutants. Lag01and lag0–2 refer to moving average of the current day and preceding 1 day and 2 days. IIV denote the models IIV. (TIFF 551 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Tensor product
 QuasiPoisson regression
 Collinearity
 Health effects
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Please note that comments may be removed without notice if they are flagged by another user or do not comply with our community guidelines.