Distributed lag interaction models with two pollutants

Distributed lag models (DLMs) have been widely used in environmental epidemiology to quantify the lagged effects of air pollution on a health outcome of interest such as mortality and morbidity. Most previous DLM approaches only consider one pollutant at a time. In this article, we propose distributed lag interaction model (DLIM) to characterize the joint lagged effect of two pollutants. One natural way to model the interaction surface is by assuming that the underlying basis functions are tensor products of the basis functions that generate the main-effect distributed lag functions. We extend Tukey's one-degree-of-freedom interaction structure to the two-dimensional DLM context. We also consider shrinkage versions of the two to allow departure from the specified Tukey's interaction structure and achieve bias-variance tradeoff. We derive the marginal lag effects of one pollutant when the other pollutant is fixed at certain quantiles. In a simulation study, we show that the shrinkage methods have better average performance in terms of mean squared error (MSE) across different scenarios. We illustrate the proposed methods by using the National Morbidity, Mortality, and Air Pollution Study (NMMAPS) data to model the joint effects of PM10 and O3 on mortality count in Chicago, Illinois, from 1987 to 2000.

[1]  F. Dominici,et al.  Does the Effect of PM10 on Mortality Depend on PM Nickel and Vanadium Content? A Reanalysis of the NMMAPS Data , 2007, Environmental health perspectives.

[2]  S L Zeger,et al.  Bayesian Distributed Lag Models: Estimating Effects of Particulate Matter Air Pollution on Daily Mortality , 2009, Biometrics.

[3]  Arnab Maity,et al.  Testing in semiparametric models with interaction, with applications to gene-environment interactions. , 2009, Journal of the Royal Statistical Society. Series B, Statistical methodology.

[4]  Andrew O. Finley,et al.  Improving the performance of predictive process modeling for large datasets , 2009, Comput. Stat. Data Anal..

[5]  V. Muggeo.  Bivariate distributed lag models for the analysis of temperature‐by‐pollutant interaction effect on mortality , 2007 .

[6]  A Gasparrini,et al.  Distributed lag non-linear models , 2010, Statistics in medicine.

[7]  Andrew Thomas,et al.  The BUGS project: Evolution, critique and future directions , 2009, Statistics in medicine.

[8]  Joe L. Mauderly,et al.  Is There Evidence for Synergy Among Air Pollutants in Causing Health Effects? , 2008, Environmental health perspectives.

[9]  Paul H. C. Eilers,et al.  Direct generalized additive modeling with penalized likelihood , 1998 .

[10]  Shirley Almon.  The Distributed Lag Between Capital Appropriations and Expenditures , 1965 .

[11]  Isabella Annesi-Maesano,et al.  Estimating the health effects of exposure to multi-pollutant mixture. , 2012, Annals of epidemiology.

[12]  S. Zeger,et al.  Are the acute effects of particulate matter on mortality in the National Morbidity, Mortality, and Air Pollution Study the result of inadequate control for weather and season? A sensitivity analysis using flexible distributed lag models. , 2005, American journal of epidemiology.

[13]  Bhramar Mukherjee,et al.  Set‐based tests for genetic association in longitudinal studies , 2015, Biometrics.

[14]  S. Roberts.  An Investigation of Distributed Lag Models in the Context of Air Pollution and Mortality Time Series Analysis , 2005, Journal of the Air & Waste Management Association.

[15]  Antonella Zanobetti,et al.  The Temporal Pattern of Mortality Responses to Air Pollution: A Multicity Assessment of Mortality Displacement , 2002, Epidemiology.

[16]  A. Gelfand,et al.  Gaussian predictive process models for large spatial data sets , 2008, Journal of the Royal Statistical Society. Series B, Statistical methodology.

[17]  M. Plummer.  Penalized loss functions for Bayesian model comparison. , 2008, Biostatistics.

[18]  Francesca Dominici,et al.  Revised Analyses of the National Morbidity, Mortality, and Air Pollution Study: Mortality Among Residents Of 90 Cities , 2005, Journal of toxicology and environmental health. Part A.

[19]  Bhramar Mukherjee,et al.  Statistical strategies for constructing health risk models with multiple pollutants and their interactions: possible choices and comparisons , 2013, Environmental Health.

[20]  Joel Schwartz,et al.  REVIEW OF EPIDEMIOLOGICAL EVIDENCE OF HEALTH EFFECTS OF PARTICULATE AIR POLLUTION , 1995 .

[21]  J. Tukey.  One Degree of Freedom for Non-Additivity , 1949 .

[22]  David C Christiani,et al.  Bayesian kernel machine regression for estimating the health effects of multi-pollutant mixtures. , 2015, Biostatistics.

[23]  Bradley P. Carlin,et al.  Bayesian measures of model complexity and fit , 2002 .

[24]  A. Gelman.  Prior distributions for variance parameters in hierarchical models (comment on article by Browne and Draper) , 2004 .

[25]  M P Wand,et al.  Generalized additive distributed lag models: quantifying mortality displacement. , 2000, Biostatistics.

[26]  F. Dominici,et al.  Fine particulate air pollution and hospital admission for cardiovascular and respiratory diseases. , 2006, JAMA.

[27]  D. Dockery,et al.  Health Effects of Fine Particulate Air Pollution: Lines that Connect , 2006, Journal of the Air & Waste Management Association.

[28]  Matthew J. Heaton,et al.  Extending distributed lag models to higher degrees. , 2014, Biostatistics.

[29]  R. Tibshirani.  Regression Shrinkage and Selection via the Lasso , 1996 .

[30]  N. Chatterjee,et al.  Powerful multilocus tests of genetic association in the presence of gene-gene and gene-environment interactions. , 2006, American journal of human genetics.

[31]  Kerrie Mengersen,et al.  Temperature, air pollution and total mortality during summers in Sydney, 1994–2004 , 2008, International journal of biometeorology.

[32]  Christopher D. Barr,et al.  Protecting Human Health From Air Pollution: Shifting From a Single-pollutant to a Multipollutant Approach , 2010, Epidemiology.