Generalized additive models with principal component analysis: an application to time series of respiratory disease and air pollution data

作     者：de Souza, Juliana B. Reisen, Valderio A. Franco, Glaura C. Ispany, Marton Bondon, Pascal Santos, Jane Meri 

作者机构：Univ Fed Espirito Santo Vitoria Spain Cent Supelec Gif Sur Yvette France Univ Fed Minas Gerais Belo Horizonte MG Brazil Univ Debrecen Debrecen Hungary Ctr Natl Rech Sci Gif Sur Yvette France Univ Paris Saclay Paris France 

出 版 物：《JOURNAL OF THE ROYAL STATISTICAL SOCIETY SERIES C-APPLIED STATISTICS》 (皇家统计学会志，C辑：应用统计学)

年 卷 期：2018年第67卷第2期

页      面：453-480页

核心收录：

学科分类：0202[经济学-应用经济学] 02[经济学] 020208[经济学-统计学] 07[理学] 0714[理学-统计学（可授理学、经济学学位）] 

基　　金：National Council for Scientific and Technological Development (the Conselho Nacional de Desenvolvimento Cientifico e Tecnologico) Brazilian Federal Agency for the Support and Evaluation of Graduate Education (the Coordenacao de Aperfeicoamento de Pessoal de Nivel Superior) Espirito Santo State Research Foundation (Fundacao de Amparo a Pesquisa do Espirito Santo) Minas Gerais State Research Foundation (the Fundacao de Amparo a Pesquisa do Estado de Minas Gerais) CentraleSupelec 

主　　题：Generalized additive model Multicollinearity Principal component analysis Relative risk Serial correlation Vector auto-regressive model 

摘      要：Environmental epidemiological studies of the health effects of air pollution frequently utilize the generalized additive model (GAM) as the standard statistical methodology, considering the ambient air pollutants as explanatory covariates. Although exposure to air pollutants is multi-dimensional, the majority of these studies consider only a single pollutant as a covariate in the GAM model. This model restriction may be because the pollutant variables do not only have serial dependence but also interdependence between themselves. In an attempt to convey a more realistic model, we propose here the hybrid generalized additive model-principal component analysis-vector auto-regressive (GAM-PCA-VAR) model, which is a combination of PCA and GAMs along with a VAR process. The PCA is used to eliminate the multicollinearity between the pollutants whereas the VAR model is used to handle the serial correlation of the data to produce white noise processes as covariates in the GAM. Some theoretical and simulation results of the methodology proposed are discussed, with special attention to the effect of time correlation of the covariates on the PCA and, consequently, on the estimates of the parameters in the GAM and on the relative risk, which is a commonly used statistical quantity to measure the effect of the covariates, especially the pollutants, on population health. As a main motivation to the methodology, a real data set is analysed with the aim of quantifying the association between respiratory disease and air pollution concentrations, especially particulate matter PM10, sulphur dioxide, nitrogen dioxide, carbon monoxide and ozone. The empirical results show that the GAM-PCA-VAR model can remove the auto-correlations from the principal components. In addition, this method produces estimates of the relative risk, for each pollutant, which are not affected by the serial correlation in the data. This, in general, leads to more pronounced values of the estimated risk compared