Principal components and generalized linear modeling in the correlation between hospital admissions and air pollution

OBJECTIVE To analyze the association between concentrations of air pollutants and admissions for respiratory causes in children. METHODS Ecological time series study. Daily figures for hospital admissions of children aged < 6, and daily concentrations of air pollutants (PM10, SO2, NO2, O3 and CO) were analyzed in the Região da Grande Vitória, ES, Southeastern Brazil, from January 2005 to December 2010. For statistical analysis, two techniques were combined: Poisson regression with generalized additive models and principal model component analysis. Those analysis techniques complemented each other and provided more significant estimates in the estimation of relative risk. The models were adjusted for temporal trend, seasonality, day of the week, meteorological factors and autocorrelation. In the final adjustment of the model, it was necessary to include models of the Autoregressive Moving Average Models (p, q) type in the residuals in order to eliminate the autocorrelation structures present in the components. RESULTS For every 10:49 μg/m3 increase (interquartile range) in levels of the pollutant PM10 there was a 3.0% increase in the relative risk estimated using the generalized additive model analysis of main components-seasonal autoregressive – while in the usual generalized additive model, the estimate was 2.0%. CONCLUSIONS Compared to the usual generalized additive model, in general, the proposed aspect of generalized additive model − principal component analysis, showed better results in estimating relative risk and quality of fit.

The effects of air pollution on cardiovascular diseases: lag structures

OBJECTIVE: To assess the lag structure between air pollution exposure and elderly cardiovascular diseases hospital admissions, by gender. METHODS: Health data of people aged 64 years or older was stratified by gender in São Paulo city, Southeastern Brazil, from 1996 to 2001. Daily levels of air pollutants (CO, PM10, O3, NO2, and SO2) , minimum temperature, and relative humidity were also analyzed. It were fitted generalized additive Poisson regressions and used constrained distributed lag models adjusted for long time trend, weekdays, weather and holidays to assess the lagged effects of air pollutants on hospital admissions up to 20 days after exposure. RESULTS: Interquartile range increases in PM10 (26.21 mug/m³) and SO2 (10.73 mug/m³) were associated with 3.17% (95% CI: 2.09-4.25) increase in congestive heart failure and 0.89% (95% CI: 0.18-1.61) increase in total cardiovascular diseases at lag 0, respectively. Effects were higher among female group for most of the analyzed outcomes. Effects of air pollutants for different outcomes and gender groups were predominately acute and some "harvesting" were found. CONLUSIONS: The results show that cardiovascular diseases in São Paulo are strongly affected by air pollution.

Air pollution and neonatal deaths in São Paulo, Brazil

Air pollution has been associated with health effects on different age groups. The present study was designed to assess the impact of daily changes in air pollutants (NO2, SO2, CO, O3, and particle matter (PM10)) on total number of daily neonatal deaths (those that occur between the first and the 28th days of life) in São Paulo, from January 1998 to December 2000, since adverse outcomes such as neonatal deaths associated with air pollution in Brazil have not been evaluated before. Generalized additive Poisson regression models were used and nonparametric smooth functions (loess) were adopted to control long-term trend, temperature, humidity, and short-term trends. A linear term was used for holidays. The association between air pollutants and neonatal deaths showed a short time lag. Interquartile range increases in PM10 (23.3 µg/m³) and SO2 (9.2 µg/m³) were associated with increases of 4% (95% CI, 2-6) and 6% (95% CI, 4-8), respectively. Instead of adopting a two-pollutant model we created an index to represent PM10 and SO2 effects. For an interquartile range increase in the index an increase of 6.3% (95% CI, 6.1-6.5) in neonatal deaths was observed. These results agree with previous studies performed by our group showing the deleterious effects of air pollutants during the perinatal period. The method reported here represents an alternative approach to analyze the relationship between highly correlated pollutants and public health problems, reinforcing the idea of the synergic effects of air pollutants in public health.

Effect of air pollution on pediatric respiratory emergency room visits and hospital admissions

In order to assess the effect of air pollution on pediatric respiratory morbidity, we carried out a time series study using daily levels of PM10, SO2, NO2, ozone, and CO and daily numbers of pediatric respiratory emergency room visits and hospital admissions at the Children's Institute of the University of São Paulo Medical School, from August 1996 to August 1997. In this period there were 43,635 hospital emergency room visits, 4534 of which were due to lower respiratory tract disease. The total number of hospital admissions was 6785, 1021 of which were due to lower respiratory tract infectious and/or obstructive diseases. The three health end-points under investigation were the daily number of emergency room visits due to lower respiratory tract diseases, hospital admissions due to pneumonia, and hospital admissions due to asthma or bronchiolitis. Generalized additive Poisson regression models were fitted, controlling for smooth functions of time, temperature and humidity, and an indicator of weekdays. NO2 was positively associated with all outcomes. Interquartile range increases (65.04 µg/m³) in NO2 moving averages were associated with an 18.4% increase (95% confidence interval, 95% CI = 12.5-24.3) in emergency room visits due to lower respiratory tract diseases (4-day moving average), a 17.6% increase (95% CI = 3.3-32.7) in hospital admissions due to pneumonia or bronchopneumonia (3-day moving average), and a 31.4% increase (95% CI = 7.2-55.7) in hospital admissions due to asthma or bronchiolitis (2-day moving average). The study showed that air pollution considerably affects children's respiratory morbidity, deserving attention from the health authorities.

Association between NOx exposure and deaths caused by respiratory diseases in a medium-sized Brazilian city

Exposure to nitrogen oxides (NOx) emitted by burning fossil fuels has been associated with respiratory diseases. We aimed to estimate the effects of NOx exposure on mortality owing to respiratory diseases in residents of Taubaté, São Paulo, Brazil, of all ages and both sexes. This time-series ecological study from August 1, 2011 to July 31, 2012 used information on deaths caused by respiratory diseases obtained from the Health Department of Taubaté. Estimated daily levels of pollutants (NOx, particulate matter, ozone, carbon monoxide) were obtained from the Centro de Previsão de Tempo e Estudos Climáticos Coupled Aerosol and Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modeling System. These environmental variables were used to adjust the multipollutant model for apparent temperature. To estimate association between hospitalizations owing to asthma and air pollutants, generalized additive Poisson regression models were developed, with lags as much as 5 days. There were 385 deaths with a daily mean (±SD) of 1.05±1.03 (range: 0-5). Exposure to NOx was significantly associated with mortality owing to respiratory diseases: relative risk (RR)=1.035 (95% confidence interval [CI]: 1.008-1.063) for lag 2, RR=1.064 (95%CI: 1.017-1.112) lag 3, RR=1.055 (95%CI: 1.025-1.085) lag 4, and RR=1.042 (95%CI: 1.010-1.076) lag 5. A 3 µg/m3 reduction in NOx concentration resulted in a decrease of 10-18 percentage points in risk of death caused by respiratory diseases. Even at NOx concentrations below the acceptable standard, there is association with deaths caused by respiratory diseases.

Air pollution effects on myocardial infarction

OBJECTIVE: Myocardial infarction is an acute and severe cardiovascular disease that generally leads to patient admissions to intensive care units and few cases are initially admitted to infirmaries. The objective of the study was to assess whether estimates of air pollution effects on myocardial infarction morbidity are modified by the source of health information. METHODS: The study was carried out in hospitals of the Brazilian Health System in the city of São Paulo, Southern Brazil. A time series study (1998-1999) was performed using two outcomes: infarction admissions to infirmaries and to intensive care units, both for people older than 64 years of age. Generalized linear models controlling for seasonality (long and short-term trends) and weather were used. The eight-day cumulative effects of air pollutants were assessed using third degree polynomial distributed lag models. RESULTS: Almost 70% of daily hospital admissions due to myocardial infarction were to infirmaries. Despite that, the effects of air pollutants on infarction were higher for intensive care units admissions. All pollutants were positively associated with the study outcomes but SO2 presented the strongest statistically significant association. An interquartile range increase on SO2 concentration was associated with increases of 13% (95% CI: 6-19) and 8% (95% CI: 2-13) of intensive care units and infirmary infarction admissions, respectively. CONCLUSIONS: It may be assumed there is a misclassification of myocardial infarction admissions to infirmaries leading to overestimation. Also, despite the absolute number of events, admissions to intensive care units data provides a more adequate estimate of the magnitude of air pollution effects on infarction admissions.

The AIDS epidemic in the Amazon region: a spatial case-control study in Rondonia, Brazil

OBJECTIVE To analyze spatial changes in the risk of AIDS and the relationship between AIDS incidence and socioeconomic variables in the state of Rondonia, Amazon region. METHODS A spatial, population case-control study in Rondonia, Brazil, based on 1,780 cases reported to the Epidemiological Surveillance System and controls based on demographic data from 1987 to 2006. The cases were grouped into five consecutive four-year periods. A generalized additive model was adjusted to the data; the dependent variable was the status of the individuals (case or control), and the independent variables were a bi-dimensional spline of the geographic coordinates and some municipality-level socioeconomic variables. The observed values of the Moran’s I test were compared to a reference distribution of values generated under conditions of spatial randomness. RESULTS AIDS risk shows a marked spatial and temporal pattern. The disease incidence is related to socioeconomic variables at the municipal level in Rondônia, such as urbanization and human capital. The highest incidence rates of AIDS are in municipalities along the BR-364 highway and calculations of the Moran’s I test show positive spatial correlation associated with proximity of the municipality to the highway in the third and fourth periods (p = 0.05). CONCLUSIONS Incidence of the disease is higher in municipalities of greater economic wealth and urbanization, and in those municipalities bisected by Rondônia’s main roads. The rapid development associated with the opening up of once remote regions may be accompanied by an increase in these risks to health.

Factors associated with vaccination coverage in children < 5 years in Angola

OBJECTIVE To analyze vaccination coverage and factors associated with a complete immunization scheme in children < 5 years old. METHODS This cross-sectional household census survey evaluated 1,209 children < 5 years old living in Bom Jesus, Angola, in 2010. Data were obtained from interviews, questionnaires, child immunization histories, and maternal health histories. The statistical analysis used generalized linear models, in which the dependent variable followed a binary distribution (vaccinated, unvaccinated) and the association function was logarithmic and had the children’s individual, familial, and socioeconomic factors as independent variables. RESULTS Vaccination coverage was 37.0%, higher in children < 1 year (55.0%) and heterogeneous across neighborhoods; 52.0% of children of both sexes had no immunization records. The prevalence rate of vaccination significantly varied according to child age, mother’s level of education, family size, ownership of household appliances, and destination of domestic waste. CONCLUSIONS Vulnerable groups with vaccination coverage below recommended levels continue to be present. Some factors indicate inequalities that represent barriers to full immunization, indicating the need to implement more equitable policies. The knowledge of these factors contributes to planning immunization promotion measures that focus on the most vulnerable groups.

Occurrence of severe dengue in Rio de Janeiro: an ecological study

Introduction This study aimed to analyze the relationship between the incidence of severe dengue during the 2008 epidemic in Rio de Janeiro, Brazil, and socioeconomic indicators, as well as indicators of health service availability and previous circulation of the dengue virus serotype-3 (DENV-3). Methods In this ecological study, the units of analysis were the districts of Rio de Janeiro. The data were incorporated into generalized linear models, and the incidence of severe dengue in each district was the outcome variable. Results The districts with more cases of dengue fever in the 2001 epidemic and a higher percentage of residents who declared their skin color or race as black had higher incidence rates of severe dengue in the 2008 epidemic [incidence rate ratio (IRR)= 1.21; 95% confidence interval (95%CI)= 1.05-1.40 and IRR= 1.34; 95%CI= 1.16-1.54, respectively]. In contrast, the districts with Family Health Strategy (FHS) clinics were more likely to have lower incidence rates of severe dengue in the 2008 epidemic (IRR= 0.81; 95%CI= 0.70-0.93). Conclusions At the ecological level, our findings suggest the persistence of health inequalities in this region of Brazil that are possibly due to greater social vulnerability among the self-declared black population. Additionally, the protective effect of FHS clinics may be due to the ease of access to other levels of care in the health system or to a reduced vulnerability to dengue transmission that is afforded by local practices to promote health.

Factors associated with involuntary hospital admissions in technology-dependent children

Abstract OBJECTIVE To identify the factors associated with involuntary hospital admissions of technology-dependent children, in the municipality of Ribeirão Preto, São Paulo State, Brazil. METHOD A cross-sectional study, with a quantitative approach. After an active search, 124 children who qualified under the inclusion criteria, that is to say, children from birth to age 12, were identified. Data was collected in home visits to mothers or the people responsible for the children, through the application of a questionnaire. Analysis of the data followed the assumptions of the Generalized Linear Models technique. RESULTS 102 technology-dependent children aged between 6 months and 12 years participated in the study, of whom 57% were male. The average number of involuntary hospital admissions in the previous year among the children studied was 0.71 (±1.29). In the final model the following variables were significantly associated with the outcome: age (OR=0.991; CI95%=0.985-0.997), and the number of devices (OR=0.387; CI95%=0.219-0.684), which were characterized as factors of protection and quantity of medications (OR=1.532; CI95%=1.297-1.810), representing a risk factor for involuntary hospital admissions in technology-dependent children. CONCLUSION The results constitute input data for consideration of the process of care for technology-dependent children by supplying an explanatory model for involuntary hospital admissions for this client group.

The effect of environment on development and survival of pupae of the necrophagous fly Ophyra albuquerquei Lopes (Diptera, Muscidae)

The effect of environment on development and survival of pupae of the necrophagous fly Ophyra albuquerquei Lopes (Diptera, Muscidae). Species of Ophyra Robineau-Desvoidy, 1830 are found in decomposing bodies, usually in fresh, bloated and decay stages. Ophyra albuquerquei Lopes, for example, can be found in animal carcasses. The influence of environmental factors has not been evaluated in puparia of O. albuquerquei. Thus, the focus of this work was motivated by the need for models to predict the development of a necrophagous insect as a function of abiotic factors. Colonies of O. albuquerquei were maintained in the laboratory to obtain pupae. On the tenth day of each month 200 pupae, divided equally into 10 glass jars, were exposed to the environment and checked daily for adult emergence of each sample. We concluded that the high survival rate observed suggested that the diets used for rearing the larvae and maintaining the adults were appropriate. Also, the data adjusted to robust generalized linear models and there were no interruptions of O. albuquerquei pupae development within the limits of temperatures studied in southern Rio Grande do Sul, given the high survival presented.

A time series analysis of meteorological factors and hospital outpatient admissions for cardiovascular disease in the Northern district of Guizhou Province, China

Findings on the effects of weather on health, especially the effects of ambient temperature on overall morbidity, remain inconsistent. We conducted a time series study to examine the acute effects of meteorological factors (mainly air temperature) on daily hospital outpatient admissions for cardiovascular disease (CVD) in Zunyi City, China, from January 1, 2007 to November 30, 2009. We used the generalized additive model with penalized splines to analyze hospital outpatient admissions, climatic parameters, and covariate data. Results show that, in Zunyi, air temperature was associated with hospital outpatient admission for CVD. When air temperature was less than 10°C, hospital outpatient admissions for CVD increased 1.07-fold with each increase of 1°C, and when air temperature was more than 10°C, an increase in air temperature by 1°C was associated with a 0.99-fold decrease in hospital outpatient admissions for CVD over the previous year. Our analyses provided statistically significant evidence that in China meteorological factors have adverse effects on the health of the general population. Further research with consistent methodology is needed to clarify the magnitude of these effects and to show which populations and individuals are vulnerable.

Degree of multicollinearity and variables involved in linear dependence in additive-dominant models

The objective of this work was to assess the degree of multicollinearity and to identify the variables involved in linear dependence relations in additive-dominant models. Data of birth weight (n=141,567), yearling weight (n=58,124), and scrotal circumference (n=20,371) of Montana Tropical composite cattle were used. Diagnosis of multicollinearity was based on the variance inflation factor (VIF) and on the evaluation of the condition indexes and eigenvalues from the correlation matrix among explanatory variables. The first model studied (RM) included the fixed effect of dam age class at calving and the covariates associated to the direct and maternal additive and non-additive effects. The second model (R) included all the effects of the RM model except the maternal additive effects. Multicollinearity was detected in both models for all traits considered, with VIF values of 1.03 - 70.20 for RM and 1.03 - 60.70 for R. Collinearity increased with the increase of variables in the model and the decrease in the number of observations, and it was classified as weak, with condition index values between 10.00 and 26.77. In general, the variables associated with additive and non-additive effects were involved in multicollinearity, partially due to the natural connection between these covariables as fractions of the biological types in breed composition.

Building predictive models of soil particle-size distribution

Is it possible to build predictive models (PMs) of soil particle-size distribution (psd) in a region with complex geology and a young and unstable land-surface? The main objective of this study was to answer this question. A set of 339 soil samples from a small slope catchment in Southern Brazil was used to build PMs of psd in the surface soil layer. Multiple linear regression models were constructed using terrain attributes (elevation, slope, catchment area, convergence index, and topographic wetness index). The PMs explained more than half of the data variance. This performance is similar to (or even better than) that of the conventional soil mapping approach. For some size fractions, the PM performance can reach 70 %. Largest uncertainties were observed in geologically more complex areas. Therefore, significant improvements in the predictions can only be achieved if accurate geological data is made available. Meanwhile, PMs built on terrain attributes are efficient in predicting the particle-size distribution (psd) of soils in regions of complex geology.

Random regression models to estimate genetic parameters for milk production of Guzerat cows using orthogonal Legendre polynomials

The objective of this work was to compare random regression models for the estimation of genetic parameters for Guzerat milk production, using orthogonal Legendre polynomials. Records (20,524) of test-day milk yield (TDMY) from 2,816 first-lactation Guzerat cows were used. TDMY grouped into 10-monthly classes were analyzed for additive genetic effect and for environmental and residual permanent effects (random effects), whereas the contemporary group, calving age (linear and quadratic effects) and mean lactation curve were analized as fixed effects. Trajectories for the additive genetic and permanent environmental effects were modeled by means of a covariance function employing orthogonal Legendre polynomials ranging from the second to the fifth order. Residual variances were considered in one, four, six, or ten variance classes. The best model had six residual variance classes. The heritability estimates for the TDMY records varied from 0.19 to 0.32. The random regression model that used a second-order Legendre polynomial for the additive genetic effect, and a fifth-order polynomial for the permanent environmental effect is adequate for comparison by the main employed criteria. The model with a second-order Legendre polynomial for the additive genetic effect, and that with a fourth-order for the permanent environmental effect could also be employed in these analyses.