Influence diagnostics for polyhazard models in the presence of covariates

In this paper, we present various diagnostic methods for polyhazard models. Polyhazard models are a flexible family for fitting lifetime data. Their main advantage over the single hazard models, such as the Weibull and the log-logistic models, is to include a large amount of nonmonotone hazard shapes, as bathtub and multimodal curves. Some influence methods, such as the local influence and total local influence of an individual are derived, analyzed and discussed. A discussion of the computation of the likelihood displacement as well as the normal curvature in the local influence method are presented. Finally, an example with real data is given for illustration.

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.

Technical and economical evaluation of a harvester, working under different spacing and planting arrangement conditions in eucalypts plantations

Plantation spacing selection has the primary objective of assigning each tree enough space for maximum growth and best quality to be attained with a minimum cost. From the harvest standpoint, an increase in stand density directly implies a decrease of individual tree volume, reducing also harvester productive capacity. The objective of this research is to assess the effects of several initial spacings and arrangements in eucalyptus plantations on production capacity, operational capacity and costs of forest harvester. Real operational data were collected from two eucalypt plantations at different initial spacing of 6.0, 7.5, 9.0, 12 and 18 m(2) per tree. Simulation data were obtained from a forest harvester simulator. Using spacing (E), mean tree volume (MV), diameter at breast height (DBH) and height (H) values, a stepwise regression test procedure was run, and correlations computed in order to measure their participation in operational capacity. Operational costs were computed with an accounting method proposed by FAO. Mean tree volume (MV) explained 88% of forest harvester operational capacity. Spacing (E) affected 8.5% of harvester operational capacity; wider spacings were related to higher individual tree volumes. Harvesting operation costs were lower in wider spaced treatments.

Unsaturated hydraulic conductivity of an Oxisol, using a neutron probe.

Unsaturated hydraulic conductivity of an Oxisol, using a neutron probe. The objective of this study was to determine the unsaturated hydraulic conductivity, using a neutron probe, of a clay sandy Oxisol. The Study was carried out in the city of Piracicaba, kite of Sao Paulo, Brazil (22 degrees 42` 43.3 `` S, 47 degrees`37` 10.4 `` W, 546 m). The dimensions of the experimental plot were 45 In x 15 m, in which 40 aluminum tubes were installed in order to access a neutron probe to measure the soil water content at the depths of 0.2, 0.4, 0.6, 0.8 and 1.0 m and, then, calculate the soil water storage of the 0 - 1.0 m soil layer. The distribution of these tubes was made in grids of four columns by ten rows in spacing of 5 x 5 m. The K(theta) functions were determined in the 40 points from regression analyses of theta as function Int and h(z) as a function of Int, being K the hydraulic conductivity, theta the volumetric soil water content, h(z) the soil water storage in the 0 - Z m layer, and t the soil water redistribution time. The neutron probe proved to be an efficient equipment in determining soil water contents, in the instantaneous profile method for determination of the K(theta) function in homogeneous soil.

Definition of priority areas for forest conservation through the ordered weighted averaging method

The general objective of this study was to evaluate the ordered weighted averaging (OWA) method, integrated to a geographic information systems (GIS), in the definition of priority areas for forest conservation in a Brazilian river basin, aiming at to increase the regional biodiversity. We demonstrated how one could obtain a range of alternatives by applying OWA, including the one obtained by the weighted linear combination method and, also the use of the analytic hierarchy process (AHP) to structure the decision problem and to assign the importance to each criterion. The criteria considered important to this study were: proximity to forest patches; proximity among forest patches with larger core area; proximity to surface water; distance from roads: distance from urban areas; and vulnerability to erosion. OWA requires two sets of criteria weights: the weights of relative criterion importance and the order weights. Thus, Participatory Technique was used to define the criteria set and the criterion importance (based in AHP). In order to obtain the second set of weights we considered the influence of each criterion, as well as the importance of each one, on this decision-making process. The sensitivity analysis indicated coherence among the criterion importance weights, the order weights, and the solution. According to this analysis, only the proximity to surface water criterion is not important to identify priority areas for forest conservation. Finally, we can highlight that the OWA method is flexible, easy to be implemented and, mainly, it facilitates a better understanding of the alternative land-use suitability patterns. (C) 2008 Elsevier B.V. All rights reserved.

Genetic Diversity and a PCR-Based Method for Xanthomonas axonopodis Detection in Passion Fruit

Xanthomonas axonopodis pv. passiflorae causes bacterial spot in passion fruit. It attacks the purple and yellow passion fruit as well as the sweet passion fruit. The diversity of 87 isolates of pv. passiflorae collected from across 22 fruit orchards in Brazil was evaluated using molecular profiles and statistical procedures, including an unweighted pair-group method with arithmetical averages-based dendrogram, analysis of molecular variance (AMOVA), and an assigning test that provides information on genetic structure at the population level. Isolates from another eight pathovars were included in the molecular analyses and all were shown to have a distinct repetitive sequence-based polymerase chain reaction profile. Amplified fragment length polymorphism technique revealed considerable diversity among isolates of pv. passiflorae, and AMOVA showed that most of the variance (49.4%) was due to differences between localities. Cluster analysis revealed that most genotypic clusters were homogeneous and that variance was associated primarily with geographic origin. The disease adversely affects fruit production and may kill infected plants. A method for rapid diagnosis of the pathogen, even before the disease symptoms become evident, has value for producers. Here, a set of primers (Xapas) was designed by exploiting a single-nucleotide polymorphism between the sequences of the intergenic 16S-23S rRNA spacer region of the pathovars. Xapas was shown to effectively detect all pv. passiflorae isolates and is recommended for disease diagnosis in passion fruit orchards.

Endophytic Colonization of Potato (Solanum tuberosum L.) by a Novel Competent Bacterial Endophyte, Pseudomonas putida Strain P9, and Its Effect on Associated Bacterial Communities

Pseudomonas putida strain P9 is a novel competent endophyte from potato. P9 causes cultivar-dependent suppression of Phytophthora infestans. Colonization of the rhizoplane and endosphere of potato plants by P9 and its rifampin-resistant derivative P9R was studied. The purposes of this work were to follow the fate of P9 inside growing potato plants and to establish its effect on associated microbial communities. The effects of P9 and P9R inoculation were studied in two separate experiments. The roots of transplants of three different cultivars of potato were dipped in suspensions of P9 or P9R cells, and the plants were planted in soil. The fate of both strains was followed by examining colony growth and by performing PCR-denaturing gradient gel electrophoresis (PCR-DGGE). Colonies of both strains were recovered from rhizoplane and endosphere samples of all three cultivars at two growth stages. A conspicuous band, representing P9 and P9R, was found in all Pseudomonas PCR-DGGE fingerprints for treated plants. The numbers of P9R CFU and the P9R-specific band intensities for the different replicate samples were positively correlated, as determined by linear regression analysis. The effects of plant growth stage, genotype, and the presence of P9R on associated microbial communities were examined by multivariate and unweighted-pair group method with arithmetic mean cluster analyses of PCR-DGGE fingerprints. The presence of strain P9R had an effect on bacterial groups identified as Pseudomonas azotoformans, Pseudomonas veronii, and Pseudomonas syringae. In conclusion, strain P9 is an avid colonizer of potato plants, competing with microbial populations indigenous to the potato phytosphere. Bacterization with a biocontrol agent has an important and previously unexplored effect on plant-associated communities.

Contribution of soil organic carbon to the ion exchange capacity of tropical soils

Highly weathered soils represent about 3 billion ha of the tropical region. Oxisols represent about 60% of the Brazilian territory (more than 5 million km 2), in areas of great agricultural importance. Soil organic carbon (SOC) can be responsible for more than 80% of the cation exchange capacity (CEC) of highly weathered soils, such as Oxisols and Ultisols. The objective of this study was to estimate the contribution of the SOC to the CEC of Brazilian soils from different orders. Surface samples (0.0 to 0.2 m) of 30 uncultivated soils (13 Oxisols, 6 Ultisols, 5 Alfisols, 3 Entisols, I Histosol, 1 Inceptisol. and I Molisol), under native forests and from reforestation sites from Sao Paulo State, Brazil, were collected in order to obtain a large variation of (electro)chemical, physical, and mineralogical soil attributes. Total content of SOC was quantified by titulometric and colorimetric methods. Effective cation exchange capacity (ECEC) was obtained by two methods: the indirect method-summation-estimated the ECECi from the sum of basic cations (Ca+ Mg+ K+ Na) and exchangeable Al; and the direct ECECd obtained by the compulsive exchange method, using unbuffered BaCl2 solution. The contribution of SOC to the soil CEC was estimated by the Bennema statistical method. The amount of SOC var ied from 6.6 g kg(-1) to 213.4 g kg(-1). while clay contents varied from 40 g kg(-1) to 716 g kg(-1). Soil organic carbon contents were strongly associated to the clay contents, suggesting that clay content was the primary variable in controling the variability of SOC contents in the samples. Cation exchange capacity varied from 7.0 mmol(c) kg(-1) to 137.8 mmol(c) kg(-1) and had a positive Correlation with SOC. The mean contribution (per grain) of the SOC (1.64 mmol(c)) for the soil CEC was more than 44 times higher than the contribution of the clay fraction (0.04 mmol(c),). A regression model that considered the SOC content as the only significant variable explained 60% of the variation in the soil total CEC. The importance of SOC was related to soil pedogenetic process, since its contribution to the soil CEC was more evident in Oxisols with predominance of Fe and Al (oxihydr)oxides in the mineral fraction or in Ultisols, that presented illuviated clay. The influence of SOC in the sign and in the magnitude of the net charge of soils reinforce the importance of agricultural management systems that preserve high levels of SOC, in order to improve their sustainability.

Parametric correlation functions to model the structure of permanent environmental (co)variances in milk yield random regression models

The objective of the present study was to estimate milk yield genetic parameters applying random regression models and parametric correlation functions combined with a variance function to model animal permanent environmental effects. A total of 152,145 test-day milk yields from 7,317 first lactations of Holstein cows belonging to herds located in the southeastern region of Brazil were analyzed. Test-day milk yields were divided into 44 weekly classes of days in milk. Contemporary groups were defined by herd-test-day comprising a total of 2,539 classes. The model included direct additive genetic, permanent environmental, and residual random effects. The following fixed effects were considered: contemporary group, age of cow at calving (linear and quadratic regressions), and the population average lactation curve modeled by fourth-order orthogonal Legendre polynomial. Additive genetic effects were modeled by random regression on orthogonal Legendre polynomials of days in milk, whereas permanent environmental effects were estimated using a stationary or nonstationary parametric correlation function combined with a variance function of different orders. The structure of residual variances was modeled using a step function containing 6 variance classes. The genetic parameter estimates obtained with the model using a stationary correlation function associated with a variance function to model permanent environmental effects were similar to those obtained with models employing orthogonal Legendre polynomials for the same effect. A model using a sixth-order polynomial for additive effects and a stationary parametric correlation function associated with a seventh-order variance function to model permanent environmental effects would be sufficient for data fitting.

Random regression models to estimate test-day milk yield genetic parameters Holstein cows in Southeastern Brazil

A total of 152,145 weekly test-day milk yield records from 7317 first lactations of Holstein cows distributed in 93 herds in southeastern Brazil were analyzed. Test-day milk yields were classified into 44 weekly classes of DIM. The contemporary groups were defined as herd-year-week of test-day. The model included direct additive genetic, permanent environmental and residual effects as random and fixed effects of contemporary group and age of cow at calving as covariable, linear and quadratic effects. Mean trends were modeled by a cubic regression on orthogonal polynomials of DIM. Additive genetic and permanent environmental random effects were estimated by random regression on orthogonal Legendre polynomials. Residual variances were modeled using third to seventh-order variance functions or a step function with 1, 6,13,17 and 44 variance classes. Results from Akaike`s and Schwarz`s Bayesian information criterion suggested that a model considering a 7th-order Legendre polynomial for additive effect, a 12th-order polynomial for permanent environment effect and a step function with 6 classes for residual variances, fitted best. However, a parsimonious model, with a 6th-order Legendre polynomial for additive effects and a 7th-order polynomial for permanent environmental effects, yielded very similar genetic parameter estimates. (C) 2008 Elsevier B.V. All rights reserved.

A Log-Linear Regression Model for the Beta-Weibull Distribution

We introduce the log-beta Weibull regression model based on the beta Weibull distribution (Famoye et al., 2005; Lee et al., 2007). We derive expansions for the moment generating function which do not depend on complicated functions. The new regression model represents a parametric family of models that includes as sub-models several widely known regression models that can be applied to censored survival data. We employ a frequentist analysis, a jackknife estimator, and a parametric bootstrap for the parameters of the proposed model. We derive the appropriate matrices for assessing local influences on the parameter estimates under different perturbation schemes and present some ways to assess global influences. Further, for different parameter settings, sample sizes, and censoring percentages, several simulations are performed. In addition, the empirical distribution of some modified residuals are displayed and compared with the standard normal distribution. These studies suggest that the residual analysis usually performed in normal linear regression models can be extended to a modified deviance residual in the proposed regression model applied to censored data. We define martingale and deviance residuals to evaluate the model assumptions. The extended regression model is very useful for the analysis of real data and could give more realistic fits than other special regression models.

Effectiveness of the standard evaluation method for hydraulic nozzles employed in stored grain protection trials

In stored grains, smaller depositions and great variation with respect to theoretical insecticide doses are frequently found. The objective of this work was to study the effectiveness of the standard method (ISO 5682/1-1996) employed to evaluate hydraulic nozzles used in stored corn and wheat grain protection experiments. The transversal volumetric distribution and droplet spectrum of a model TJ-60 8002EVS nozzle were determined in order to calibrate a spraying system for an application rate of 5 L/t and to obtain theoretical concentrations of 10 and 0.5 mg/kg of fenitrothion and esfenvalerate, respectively. After treatment, the corn and wheat grains were processed and deposition was analyzed by gas chromatography. The type of grain did not have any influence on insecticide deposition and was dependent upon insecticide only. The insecticide deposits on the grains only reached 42.1 and 38.2% of the intended theoretical values for fenitrothion and esfenvalerate concentrations, respectively. These results demonstrate the ineffectiveness of the standard evaluation method for hydraulic nozzles employed in stored grain protection experiments.

Biomarkers of metabolic syndrome and its relationship with the zinc nutritional status in obese women

Introduction: Obesity is a chronic disease that induces risk factors for metabolic syndrome and, is associated with disturbances in the metabolism of the zinc. Therefore, the aim of this study was to investigate the existence of relationship between the biomarkers of metabolic syndrome and the zinc nutricional status in obese women. Method: Seventy-three premenopausal women, aged between 20 and 50 years, were divided into two groups: case group, composed of obese (n = 37) and control group, composed of no obese (n = 36). The assessment of the body mass index and waist circumference were carried out using anthropometric measurements. The plasmatic and erythrocytary zinc were analyzed by method atomic absorption spectrophotometry (lambda=213.9 nm). Results: In the study, body mass index and waist circumference were higher in obese women than control group (p < 0.05). The mean plasmatic zinc was 72.2 +/- 9.0 mu g/dl in obese women and 73.4 +/- 8.5 mu g/dl in control group (p > 0.05). The mean erythrocytary zinc was 36.4 +/- 15.0 mu g/gHb and 45.4 +/- 14.3 mu g/gHb in the obese and controls, respectively (p < 0.05). Regression analysis showed that the body mass index (t=-2.85) and waist circumference (t=-2.37) have a negative relationship only with the erythrocytary zinc (R(2)=0.32, p < 0.05). Conclusions: The study shows that there are alterations in biochemical parameters of zinc in obese women, with low zinc concentrations in erythrocytes. Regression analysis demonstrates that the erythrocytary zinc is influenced by biomarkers of the metabolic syndrome, presenting an inverse relationship with the waist circumference and body mass index.