Log-modified Weibull regression models with censored data: Sensitivity and residual analysis

This paper proposes a regression model considering the modified Weibull distribution. This distribution can be used to model bathtub-shaped failure rate functions. Assuming censored data, we consider maximum likelihood and Jackknife estimators for the parameters of the model. We derive the appropriate matrices for assessing local influence on the parameter estimates under different perturbation schemes and we also present some ways to perform global influence. Besides, for different parameter settings, sample sizes and censoring percentages, various simulations are performed and the empirical distribution of the modified deviance residual is displayed and compared with the standard normal distribution. These studies suggest that the residual analysis usually performed in normal linear regression models can be straightforwardly extended for a martingale-type residual in log-modified Weibull regression models with censored data. Finally, we analyze a real data set under log-modified Weibull regression models. A diagnostic analysis and a model checking based on the modified deviance residual are performed to select appropriate models. (c) 2008 Elsevier B.V. All rights reserved.

On estimation and influence diagnostics for zero-inflated negative binomial regression models

The zero-inflated negative binomial model is used to account for overdispersion detected in data that are initially analyzed under the zero-Inflated Poisson model A frequentist analysis a jackknife estimator and a non-parametric bootstrap for parameter estimation of zero-inflated negative binomial regression models are considered In addition an EM-type algorithm is developed for performing maximum likelihood estimation Then the appropriate matrices for assessing local influence on the parameter estimates under different perturbation schemes and some ways to perform global influence analysis are derived In order to study departures from the error assumption as well as the presence of outliers residual analysis based on the standardized Pearson residuals is discussed The relevance of the approach is illustrated with a real data set where It is shown that zero-inflated negative binomial regression models seems to fit the data better than the Poisson counterpart (C) 2010 Elsevier B V All rights reserved

Regression models for grouped survival data: Estimation and sensitivity analysis

In this study, regression models are evaluated for grouped survival data when the effect of censoring time is considered in the model and the regression structure is modeled through four link functions. The methodology for grouped survival data is based on life tables, and the times are grouped in k intervals so that ties are eliminated. Thus, the data modeling is performed by considering the discrete models of lifetime regression. The model parameters are estimated by using the maximum likelihood and jackknife methods. To detect influential observations in the proposed models, diagnostic measures based on case deletion, which are denominated global influence, and influence measures based on small perturbations in the data or in the model, referred to as local influence, are used. In addition to those measures, the local influence and the total influential estimate are also employed. Various simulation studies are performed and compared to the performance of the four link functions of the regression models for grouped survival data for different parameter settings, sample sizes and numbers of intervals. Finally, a data set is analyzed by using the proposed regression models. (C) 2010 Elsevier B.V. All rights reserved.

A positive growth response to NaCl applications in Eucalyptus plantations established on K-deficient soils

Contrasting responses of Eucalyptus trees to K fertilizer applications have been reported on soils with low K contents. A complete randomized block experiment was set up in Brazil to test the hypothesis that large atmospheric deposits of NaCl in coastal regions might lead to a partial substitution of K by Na in Eucalyptus physiology and enhance tree growth. Treatments with application of 1.5, 3.0, 4.5 kmol K ha(-1) (K(1.5), K(3.0), 1(4.5, respectively) as KCl, 3.0 kmol K ha(-1) applied as K(2)SO(4), 3.0 kmol Na ha(-1) (Na(3.0)) as NaCl commercialized for cattle feeding, and a mixture of 1.5 kmol K + 1.5 kmol Na ha(-1) (K(1.5) + Na(1.5)) were compared to a control treatment (C) with no K and Na applications. All the plots were fertilized with large amounts of the other nutrients. A positive effect of NaCl applications on the growth of E. grandis trees was observed. NaCl and KCl additions in treatments Na(3.0) and K(3.0) increased above-ground biomass by 56% and 130% three years after planting, respectively, in comparison with the C treatment. By contrast, accumulated litterfall up to age 3 years was not significantly modified. NaCl applications in the Na(3.0) treatment significantly increased Na accumulation in above-ground tree components but did not modify K accumulation, whatever the sampling age. A partial substitution of K by Na in tree physiology, as observed for various agricultural crops, might explain this behaviour. Our results suggest the possibility of applying inexpensive K fertilizers, which are less purified in Na, and explain why high yields are achieved without K fertilizer applications in areas with large dry depositions of marine aerosols. Further investigations are necessary to identify the processes involving Na in Eucalyptus tree physiology. (C) 2009 Elsevier B.V. All rights reserved.

COMPETITION AMONG EUCALYPTUS TREES DEPENDS ON GENETIC VARIATION AND RESOURCE SUPPLY

Genetic variation and environmental heterogeneity fundamentally shape the interactions between plants of the same species. According to the resource partitioning hypothesis, competition between neighbors intensifies as their similarity increases. Such competition may change in response to increasing supplies of limiting resources. We tested the resource partitioning hypothesis in stands of genetically identical (clone-origin) and genetically diverse (seed-origin) Eucalyptus trees with different water and nutrient supplies, using individual-based tree growth models. We found that genetic variation greatly reduced competitive interactions between neighboring trees, supporting the resource partitioning hypothesis. The importance of genetic variation for Eucalyptus growth patterns depended strongly on local stand structure and focal tree size. This suggests that spatial and temporal variation in the strength of species interactions leads to reversals in the growth rank of seed-origin and clone-origin trees. This study is one of the first to experimentally test the resource partitioning hypothesis for intergenotypic vs. intragenotypic interactions in trees. We provide evidence that variation at the level of genes, and not just species, is functionally important for driving individual and community-level processes in forested ecosystems.

Spatial portability of numerical models of leaf wetness duration based on empirical approaches

Leaf wetness duration (LWD) models based on empirical approaches offer practical advantages over physically based models in agricultural applications, but their spatial portability is questionable because they may be biased to the climatic conditions under which they were developed. In our study, spatial portability of three LWD models with empirical characteristics - a RH threshold model, a decision tree model with wind speed correction, and a fuzzy logic model - was evaluated using weather data collected in Brazil, Canada, Costa Rica, Italy and the USA. The fuzzy logic model was more accurate than the other models in estimating LWD measured by painted leaf wetness sensors. The fraction of correct estimates for the fuzzy logic model was greater (0.87) than for the other models (0.85-0.86) across 28 sites where painted sensors were installed, and the degree of agreement k statistic between the model and painted sensors was greater for the fuzzy logic model (0.71) than that for the other models (0.64-0.66). Values of the k statistic for the fuzzy logic model were also less variable across sites than those of the other models. When model estimates were compared with measurements from unpainted leaf wetness sensors, the fuzzy logic model had less mean absolute error (2.5 h day(-1)) than other models (2.6-2.7 h day(-1)) after the model was calibrated for the unpainted sensors. The results suggest that the fuzzy logic model has greater spatial portability than the other models evaluated and merits further validation in comparison with physical models under a wider range of climate conditions. (C) 2010 Elsevier B.V. All rights reserved.

Root Water Extraction under Combined Water and Osmotic Stress

Using a numerical implicit model for root water extraction by a single root in a symmetric radial flow problem, based on the Richards equation and the combined convection-dispersion equation, we investigated some aspects of the response of root water uptake to combined water and osmotic stress. The model implicitly incorporates the effect of simultaneous pressure head and osmotic head on root water uptake, and does not require additional assumptions (additive or multiplicative) to derive the combined effect of water and salt stress. Simulation results showed that relative transpiration equals relative matric flux potential, which is defined as the matric flux potential calculated with an osmotic pressure head-dependent lower bound of integration, divided by the matric flux potential at the onset of limiting hydraulic conditions. In the falling rate phase, the osmotic head near the root surface was shown to increase in time due to decreasing root water extraction rates, causing a more gradual decline of relative transpiration than with water stress alone. Results furthermore show that osmotic stress effects on uptake depend on pressure head or water content, allowing a refinement of the approach in which fixed reduction factors based on the electrical conductivity of the saturated soil solution extract are used. One of the consequences is that osmotic stress is predicted to occur in situations not predicted by the saturation extract analysis approach. It is also shown that this way of combining salinity and water as stressors yields results that are different from a purely multiplicative approach. An analytical steady state solution is presented to calculate the solute content at the root surface, and compared with the outputs of the numerical model. Using the analytical solution, a method has been developed to estimate relative transpiration as a function of system parameters, which are often already used in vadose zone models: potential transpiration rate, root length density, minimum root surface pressure head, and soil theta-h and K-h functions.

Alternative Analytical Expressions for the General van Genuchten-Mualem and van Genuchten-Burdine Hydraulic Conductivity Models

The van Genuchten expressions for the unsaturated soil hydraulic properties, first published in 1980, are used frequently in various vadose zone flow and transport applications assuming a specific relationship between the m and n soil hydraulic parameters. By comparison, probably because of the complexity of the hydraulic conductivity equations, the more general solutions with independent m and n values are rarely used. We expressed the general van Genuchten-Mualem and van Genuchten-Burdine hydraulic conductivity equations in terms of hypergeometric functions, which can be approximated by infinite series that converge rapidly for relatively large values of the van Genuchten-Mualem parameter n but only very slowly when n is close to one. Alternative equations were derived that provide very close approximations of the analytical results. The newly proposed equations allow the use of independent values of the parameters m and n in the soil water retention model of van Genuchten for subsequent prediction of the van Genuchten-Mualem and van Genuchten-Burdine hydraulic conductivity models, thus providing more flexibility in fitting experimental pressure-head-dependent water content, theta(h), and hydraulic conductivity, K(h), or K(theta) data.

Forecasting fuel ethanol consumption in Brazil by time series models: 2006-2012

This article analysed scenarios for Brazilian consumption of ethanol for the period 2006 to 2012. The results show that if the country`s GDP sustains a 4.6% a year growth, domestic consumption of fuel ethanol could increase to 25.16 billion liters in this period, which is a volume relatively close to the forecasted gasoline consumption of 31 billion liters. At a lower GDP growth of 1.22% a year, gasoline consumption would be reduced and domestic ethanol consumption in Brazil would be no higher than 18.32 billion liters. Contrary to the current situation, forecasts indicated that hydrated ethanol consumption could become much higher than anhydrous consumption in Brazil. The former is being consumed in cars moved exclusively by ethanol and flex-fuel cars, successfully introduced in the country at 2003. Flex cars allow Brazilian consumers to choose between gasoline and hydrated ethanol and immediately switch to whichever fuel presents the most favourable relative price.

Cloning, expression, molecular modelling and docking analysis of glutathione transferase from Saccharum officinarum

Sugarcane yield and quality are affected by a number of biotic and abiotic stresses. In response to such stresses, plants may increase the activities of some enzymes such as glutathione transferase (GST), which are involved in the detoxification of xenobiotics. Thus, a sugarcane GST was modelled and molecular docked using the program LIGIN to investigate the contributions of the active site residues towards the binding of reduced glutathione (GSH) and 1-chloro-2,4-dinitrobenzene (CDNB). As a result, W13 and I119 were identified as key residues for the specificity of sugarcane GSTF1 (SoGSTF1) towards CDNB. To obtain a better understanding of the catalytic specificity of sugarcane GST (SoGSTF1), two mutants were designed, W13L and I119F. Tertiary structure models and the same docking procedure were performed to explain the interactions between sugarcane GSTs with GSH and CDNB. An electron-sharing network for GSH interaction was also proposed. The SoGSTF1 and the mutated gene constructions were cloned and expressed in Escherichia coli and the expressed protein purified. Kinetic analyses revealed different Km values not only for CDNB, but also for GSH. The Km values were 0.2, 1.3 and 0.3 mM for GSH, and 0.9, 1.2 and 0.5 mM for CDNB, for the wild type, W13L mutant and I119F mutant, respectively. The V(max) values were 297.6, 224.5 and 171.8 mu mol min(-1) mg(-1) protein for GSH, and 372.3, 170.6 and 160.4 mu mol min(-1) mg(-1) protein for CDNB.

Zn uptake, physiological response and stress attenuation in mycorrhizal jack bean growing in soil with increasing Zn concentrations

The influence of arbuscular mycorrhizal fungi (AMF) inoculation on Canavalia ensiformis growth. nutrient and Zn uptake, and on some physiological parameters in response to increasing soil Zn concentrations was studied. Treatments were applied in seven replicates in a 2 x 4 factorial design, consisting of the inoculation or not with the AMF Glomus etunicatum, and the addition of Zn to soil at the concentrations of 0, 100, 300 and 900 mg kg(-1). AMF inoculation enhanced the accumulation of Zn in tissues and promoted biomass yields and root nodulation. Mycorrhizal plants exhibited relative tolerance to Zn up to 300 mg kg(-1) without exhibiting visual symptoms of toxicity, in contrast to non-mycorrhizal plants which exhibited a significant growth reduction at the same soil Zn concentration. The highest concentration of Zn added to soil was highly toxic to the plants. Leaves of plants grown in high Zn concentration exhibited a Zn-induced proline accumulation and also an increase in soluble amino acid contents; however proline contents were lower in mycorrhizal jack beans. Plants in association or not with the AMF exhibited marked differences in the foliar soluble amino acid profile and composition in response to Zn addition to soil. In general, Zn induced oxidative stress which could be verified by increased lipid peroxidation rates and changes in catalase, ascorbate peroxidase, glutathione reductase and superoxide dismutase activities. In summary, G. etunicatum was able to maintain an efficient symbiosis with jack bean plants in moderately contaminated Zn-soils, improving plant performance under those conditions, which is likely to be due to a combination of physiological and nutritional changes caused by the intimate relation between fungus and plant. The enhanced Zn uptake by AMF inoculated jack bean plants might be of interest for phytoremediation purposes. (C) 2009 Elsevier Ltd. All rights reserved.

Antioxidant response of Nicotiana tabacum cv. Bright Yellow 2 cells to cadmium and nickel stress

Plant cell cultures are a suitable model system for investigation of the physiological mechanisms of tolerance to environmental stress. We have determined the effects of Cd (0.1 and 0.2 mM CdCl(2)) and Ni (0.075 and 0.75 mM NiCl(2)) on Nicotiana tabacum L. cv. Bright Yellow (TBY-2) cell suspension cultures over a 72-h period. Inhibition of growth, loss of cell viability and lipid peroxidation occurred, in general, only when the TBY-2 cells were grown at 0.2 mM CdCl(2) and at 0.75 mM NiCl(2). At 0.1 mM CdCl(2), a significant increase in growth was determined at the end of the experiment. Increases in the activities of all of the four enzymatic antioxidant defence systems tested, were induced by the two concentrations of Cd and Ni, but at different times during the period of metal exposure. Overall, the cellular antioxidant responses to Cd and Ni were similar and were apparently sufficient to avoid oxidative stress at the lower concentrations of Cd and Ni. The activities of glutathione reductase and glutathione S-transferase increased early but transiently, whereas the activities of catalase and guaiacol peroxidase increased in the latter half of the experimental period. Therefore it is likely that the metabolism of reduced glutathione was enhanced during the initial onset of the stress, while catalase and guaiacol-type peroxidase appeared to play a more important role in the antioxidant response once the stress became severe.

Does using stepwise variable selection to build sequential path analysis models make sense?

Causal inference methods - mainly path analysis and structural equation modeling - offer plant physiologists information about cause-and-effect relationships among plant traits. Recently, an unusual approach to causal inference through stepwise variable selection has been proposed and used in various works on plant physiology. The approach should not be considered correct from a biological point of view. Here, it is explained why stepwise variable selection should not be used for causal inference, and shown what strange conclusions can be drawn based upon the former analysis when one aims to interpret cause-and-effect relationships among plant traits.

Landscape unit discrimination for pedological surveys by orbital spectral response

Landscape unit discrimination for pedological surveys by orbital spectral response. The objective of tins study was compare two soil survey methods. The first was performed by methods traditionally used to distinguish landscape units and soil class discrimination. The second was based on soil class distinction through orbital spectral response. In order to establish soil characteristics and their classification, soil samples were collected at two depths in a grid system, with a distance of 500 meters between points. With these samples, physical and chemical analyses were carried out. In the sampling points, the apparent reflectance of the soil, front the orbital image, was determined and, through cluster analysis landscape units were established. In order to evaluate the resemblance reliability between the landscape units established in each method, the Kappa index was used, the value set for the confusion matrix was 0.43, indicating high quality in the comparison, showing that the non-conventional method was as close as the one carried out by photointerpretation.

Growth performance and body composition of pacu Piaractus mesopotamicus (Holmberg 1887) in response to dietary protein and energy levels

Improper dietary protein and energy levels and their ratio will lead to increased fish production cost. This work evaluated effects of dietary protein : energy ratio on growth and body composition of pacu, Piaractus mesopotamicus. Fingerling pacu (15.5 +/- 0.4 g) were fed twice a day for 10 weeks until apparent satiation with diets containing 220, 260, 300, 340 or 380 g kg-1 crude protein (CP) and 10.9, 11.7, 12.6, 13.4 or 14.2 MJ kg-1 digestible energy (DE) in a totally randomized experimental design, 5 x 5 factorial scheme (n = 3). Weight gain, specific growth rate increased and feed conversion ratio (FCR) decreased significantly (P < 0.05) when CP increased from 220 to 271, 268 and 281 g kg-1 respectively. Pacu was able to adjust feed consumption in a wide range of dietary DE concentration. Fish fed 260 CP diets showed best (P < 0.05) protein efficiency ratio and FCR with 11.7-12.6 MJ kg-1; but for the 380 CP-diets group, significant differences were observed only at 14.2 MJ kg-1 dietary energy level, suggesting that pacu favours protein as energy source. DE was the chief influence on whole body chemical composition. Minimum dietary protein requirement of pacu is 270 g kg-1, with an optimum CP : DE of 22.2 g MJ-1.