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.

Characterisation of Culex quinquefasciatus (Diptera: Culicidae) larval habitats at ground level and temporal fluctuations of larval abundance in Córdoba, Argentina

The aims of this study were to characterise the ground-level larval habitats of the mosquito Culex quinquefasciatus, to determine the relationships between habitat characteristics and larval abundance and to examine seasonal larval-stage variations in Córdoba city. Every two weeks for two years, 15 larval habitats (natural and artificial water bodies, including shallow wells, drains, retention ponds, canals and ditches) were visited and sampled for larval mosquitoes. Data regarding the water depth, temperature and pH, permanence, the presence of aquatic vegetation and the density of collected mosquito larvae were recorded. Data on the average air temperatures and accumulated precipitation during the 15 days prior to each sampling date were also obtained. Cx. quinquefasciatus larvae were collected throughout the study period and were generally most abundant in the summer season. Generalised linear mixed models indicated the average air temperature and presence of dicotyledonous aquatic vegetation as variables that served as important predictors of larval densities. Additionally, permanent breeding sites supported high larval densities. In Córdoba city and possibly in other highly populated cities at the same latitude with the same environmental conditions, control programs should focus on permanent larval habitats with aquatic vegetation during the early spring, when the Cx. quinquefasciatus population begins to increase.

She’s a femme fatale: low-density larval development produces good disease vectors

Two hypotheses for how conditions for larval mosquitoes affect vectorial capacity make opposite predictions about the relationship of adult size and frequency of infection with vector-borne pathogens. Competition among larvae produces small adult females. The competition-susceptibility hypothesis postulates that small females are more susceptible to infection and predicts frequency of infection should decrease with size. The competition-longevity hypothesis postulates that small females have lower longevity and lower probability of becoming competent to transmit the pathogen and thus predicts frequency of infection should increase with size. We tested these hypotheses for Aedes aegypti in Rio de Janeiro, Brazil, during a dengue outbreak. In the laboratory, longevity increases with size, then decreases at the largest sizes. For field-collected females, generalised linear mixed model comparisons showed that a model with a linear increase of frequency of dengue with size produced the best Akaike’s information criterion with a correction for small sample sizes (AICc). Consensus prediction of three competing models indicated that frequency of infection increases monotonically with female size, consistent with the competition-longevity hypothesis. Site frequency of infection was not significantly related to site mean size of females. Thus, our data indicate that uncrowded, low competition conditions for larvae produce the females that are most likely to be important vectors of dengue. More generally, ecological conditions, particularly crowding and intraspecific competition among larvae, are likely to affect vector-borne pathogen transmission in nature, in this case via effects on longevity of resulting adults. Heterogeneity among individual vectors in likelihood of infection is a generally important outcome of ecological conditions impacting vectors as larvae.

Vigor and field emergence in "high oleic" sunflower seeds, in southeast of Buenos Aires province

Sunflower crop was based, as yet, on high linoleic cultivars, but in the last years request for oil with higher content of oleic acid has increased, due to their dietary characteristics. At the beginning, high oleic cultivars were used to be sown in warm regions, but then the concern about growing it in temperate areas, as the south-east of Buenos Aires Province, was posed. In this region, early sowings are recommended, so that grain filling matches with a period of appropriate hydric and light conditions, as to result in greater yields. However, early sowings are limited by low soil temperature, that delays seedling emergence, resulting in heterogeneous stand establishment. The aim of this work was to evaluate seed performance of four high oleic cultivars in the southern area of Buenos Aires Province, by means of vigor tests and field trials. Germination, cold, tetrazolium viability, tetrazolium viability with cold, accelerated ageing tests and three field sowings at different soil temperatures were performed. Data were analyzed by Anova using generalised linear models, and tests and cultivars were contrasted among themselves. Similar seedling emergence under optimal and suboptimal temperatures for high oleic and high linoleic cultivars was recorded. The success of seedling establishment does not appear to be related to the acidic composition of seeds.

Computational modeling for irrigated agriculture planning. Part I: general description and linear programming

Linear programming models are effective tools to support initial or periodic planning of agricultural enterprises, requiring, however, technical coefficients that can be determined using computer simulation models. This paper, presented in two parts, deals with the development, application and tests of a methodology and of a computational modeling tool to support planning of irrigated agriculture activities. Part I aimed at the development and application, including sensitivity analysis, of a multiyear linear programming model to optimize the financial return and water use, at farm level for Jaíba irrigation scheme, Minas Gerais State, Brazil, using data on crop irrigation requirement and yield, obtained from previous simulation with MCID model. The linear programming model outputted a crop pattern to which a maximum total net present value of R$ 372,723.00 for the four years period, was obtained. Constraints on monthly water availability, labor, land and production were critical in the optimal solution. In relation to the water use optimization, it was verified that an expressive reductions on the irrigation requirements may be achieved by small reductions on the maximum total net present value.

Simulating maize phenology as a function of air temperature with a linear and a nonlinear model

The objective of this study was to adapt a nonlinear model (Wang and Engel - WE) for simulating the phenology of maize (Zea mays L.), and to evaluate this model and a linear one (thermal time), in order to predict developmental stages of a field-grown maize variety. A field experiment, during 2005/2006 and 2006/2007 was conducted in Santa Maria, RS, Brazil, in two growing seasons, with seven sowing dates each. Dates of emergence, silking, and physiological maturity of the maize variety BRS Missões were recorded in six replications in each sowing date. Data collected in 2005/2006 growing season were used to estimate the coefficients of the two models, and data collected in the 2006/2007 growing season were used as independent data set for model evaluations. The nonlinear WE model accurately predicted the date of silking and physiological maturity, and had a lower root mean square error (RMSE) than the linear (thermal time) model. The overall RMSE for silking and physiological maturity was 2.7 and 4.8 days with WE model, and 5.6 and 8.3 days with thermal time model, respectively.

Computational modeling for irrigated agriculture planning. Part II: risk analysis

Techniques of evaluation of risks coming from inherent uncertainties to the agricultural activity should accompany planning studies. The risk analysis should be carried out by risk simulation using techniques as the Monte Carlo method. This study was carried out to develop a computer program so-called P-RISCO for the application of risky simulations on linear programming models, to apply to a case study, as well to test the results comparatively to the @RISK program. In the risk analysis it was observed that the average of the output variable total net present value, U, was considerably lower than the maximum U value obtained from the linear programming model. It was also verified that the enterprise will be front to expressive risk of shortage of water in the month of April, what doesn't happen for the cropping pattern obtained by the minimization of the irrigation requirement in the months of April in the four years. The scenario analysis indicated that the sale price of the passion fruit crop exercises expressive influence on the financial performance of the enterprise. In the comparative analysis it was verified the equivalence of P-RISCO and @RISK programs in the execution of the risk simulation for the considered scenario.

Modeling and simulation of forced-air cooling of strawberries using variable convective coefficient

In the forced-air cooling process of fruits occurs, besides the convective heat transfer, the mass transfer by evaporation. The energy need in the evaporation is taken from fruit that has its temperature lowered. In this study it has been proposed the use of empirical correlations for calculating the convective heat transfer coefficient as a function of surface temperature of the strawberry during the cooling process. The aim of this variation of the convective coefficient is to compensate the effect of evaporation in the heat transfer process. Linear and exponential correlations are tested, both with two adjustable parameters. The simulations are performed using experimental conditions reported in the literature for the cooling of strawberries. The results confirm the suitability of the proposed methodology.

Substitution-Newton-Raphson method applied to the modeling of a vapour compression refrigeration system using different representations of the thermodynamic properties of R-134a

This paper gives a detailed presentation of the Substitution-Newton-Raphson method, suitable for large sparse non-linear systems. It combines the Successive Substitution method and the Newton-Raphson method in such way as to take the best advantages of both, keeping the convergence features of the Newton-Raphson with the low requirements of memory and time of the Successive Substitution schemes. The large system is solved employing few effective variables, using the greatest possible part of the model equations in substitution fashion to fix the remaining variables, but maintaining the convergence characteristics of the Newton-Raphson. The methodology is exemplified through a simple algebraic system, and applied to a simple thermodynamic, mechanical and heat transfer modeling of a single-stage vapor compression refrigeration system. Three distinct approaches for reproducing the thermodynamic properties of the refrigerant R-134a are compared: the linear interpolation from tabulated data, the use of polynomial fitted curves and the use of functions derived from the Helmholtz free energy.

Mathematical modeling of dehydration of 'Fuji' and 'Gala' apples slices using infrared

The objective of this work was to determine and model the infrared dehydration curves of apple slices - Fuji and Gala varieties. The slices were dehydrated until constant mass, in a prototype dryer with infrared heating source. The applied temperatures ranged from 50 to 100 °C. Due to the physical characteristics of the product, the dehydration curve was divided in two periods, constant and falling, separated by the critical moisture content. A linear model was used to describe the constant dehydration period. Empirical models traditionally used to model the drying behavior of agricultural products were fitted to the experimental data of the falling dehydration period. Critical moisture contents of 2.811 and 3.103 kgw kgs-1 were observed for the Fuji and Gala varieties, respectively. Based on the results, it was concluded that the constant dehydration rates presented a direct relationship with the temperature; thus, it was possible to fit a model that describes the moisture content variation in function of time and temperature. Among the tested models, which describe the falling dehydration period, the model proposed by Midilli presented the best fit for all studied conditions.

Mathematical modeling of microwave dried celery leaves and determination of the effective moisture diffusivities and activation energy

Celery (Apium graveolens L. var. secalinum Alef) leaves with 50±0.07 g weight and 91.75±0.15% humidity (~11.21 db) were dried using 8 different microwave power densities ranging between 1.8-20 W g-1, until the humidity fell down to 8.95±0.23% (~0.1 db). Microwave drying processes were completed between 5.5 and 77 min depending on the microwave power densities. In this study, measured values were compared with predicted values obtained from twenty thin layer drying theoretical, semi-empirical and empirical equations with a new thin layer drying equation. Within applied microwave power density; models whose coefficient and correlation (R²) values are highest were chosen as the best models. Weibull distribution model gave the most suitable predictions at all power density. At increasing microwave power densities, the effective moisture diffusivity values ranged from 1.595 10-10 to 6.377 10-12 m2 s-1. The activation energy was calculated using an exponential expression based on Arrhenius equation. The linear relationship between the drying rate constant and effective moisture diffusivity gave the best fit.

Testes para verificar a igualdade de parâmetros e a identidade de modelos de regressão não-linear em dados de experimento com delineamento em blocos casualizados

Considerou-se o ajustamento de equações de regressão não-linear e o teste da razão de verossimilhança, com aproximações pelas estatísticas qui-quadrado e F, para testar as hipóteses de igualdade de qualquer subconjunto de parâmetros e de identidade dos modelos para dados com repetições provenientes de experimento com delineamento em blocos completos casualizados. Concluiu-se que as duas aproximações podem ser utilizadas, mas a aproximação pela estatística F deve ser preferida, principalmente para pequenas amostras.

Atributos físicos de um Latossolo Vermelho correlacionados linear e espacialmente com a consorciação de guandu com milheto

O sistema cultivo mínimo, por possibilitar pouca movimentação de solo, menor número de operações agrícolas sem incorporação dos resíduos vegetais, apresenta vantagens em razão do menor custo de preparo e da redução das perdas de solo e água. No ano agrícola de 2006/2007, na Fazenda de Ensino e Pesquisa da Faculdade de Engenharia de Ilha Solteira, SP, Brasil - FEIS/UNESP, situada nas condições do Cerrado Brasileiro, objetivou-se analisar a produtividade de massa de matéria seca da consorciação de forragem (guandu+milheto) (MSF), em função de atributos físicos do solo, tais como resistência à penetração (RP), umidade gravimétrica (UG), umidade volumétrica (UV) e densidade do solo (DS) nas profundidades de 0,0-0,10 m; 0,10-0,20 m e 0,20-0,30 m. Para tanto, foi instalado um ensaio, contendo 117 pontos amostrais, em um Latossolo Vermelho distroférrico, sob pivô central, numa área experimental de 1600 m² sob cultivo mínimo. A análise estatística constou de análise descritiva inicial dos atributos e análise das correlações lineares simples entre eles, e, finalmente, de análise geoestatística. Do ponto de vista da correlação espacial, o atributo que mais bem explica a produtividade de massa de matéria seca da consorciação é a densidade do solo na camada de 0,20-0,30 m, com uma correlação inversa, indicando que as espécies se desenvolvem bem em solos adensados.

Regressão linear com duplo truncamento na distribuição da variável dependente

Do problema do ajuste de uma regressão linear, quando a distribuição da variável dependente tem duplo truncamento, utilizando a função de máxima verossimilhança e um processo iterativo.