The Kumaraswamy Weibull distribution with application to failure data

For the first time, we introduce and study some mathematical properties of the Kumaraswamy Weibull distribution that is a quite flexible model in analyzing positive data. It contains as special sub-models the exponentiated Weibull, exponentiated Rayleigh, exponentiated exponential, Weibull and also the new Kumaraswamy exponential distribution. We provide explicit expressions for the moments and moment generating function. We examine the asymptotic distributions of the extreme values. Explicit expressions are derived for the mean deviations, Bonferroni and Lorenz curves, reliability and Renyi entropy. The moments of the order statistics are calculated. We also discuss the estimation of the parameters by maximum likelihood. We obtain the expected information matrix. We provide applications involving two real data sets on failure times. Finally, some multivariate generalizations of the Kumaraswamy Weibull distribution are discussed. (C) 2010 The Franklin Institute. Published by Elsevier Ltd. All rights reserved.

A log-extended Weibull regression model

A bathtub-shaped failure rate function is very useful in survival analysis and reliability studies. The well-known lifetime distributions do not have this property. For the first time, we propose a location-scale regression model based on the logarithm of an extended Weibull distribution which has the ability to deal with bathtub-shaped failure rate functions. We use the method of maximum likelihood to estimate the model parameters and some inferential procedures are presented. We reanalyze a real data set under the new model and the log-modified Weibull regression model. We perform a model check based on martingale-type residuals and generated envelopes and the statistics AIC and BIC to select appropriate models. (C) 2009 Elsevier B.V. All rights reserved.

The generalized log-gamma mixture model with covariates: local influence and residual analysis

In a sample of censored survival times, the presence of an immune proportion of individuals who are not subject to death, failure or relapse, may be indicated by a relatively high number of individuals with large censored survival times. In this paper the generalized log-gamma model is modified for the possibility that long-term survivors may be present in the data. The model attempts to separately estimate the effects of covariates on the surviving fraction, that is, the proportion of the population for which the event never occurs. The logistic function is used for the regression model of the surviving fraction. Inference for the model parameters is considered via maximum likelihood. Some influence methods, such as the local influence and total local influence of an individual are derived, analyzed and discussed. Finally, a data set from the medical area is analyzed under the log-gamma generalized mixture model. A residual analysis is performed in order to select an appropriate model.

A generalized modified Weibull distribution for lifetime modeling

A four parameter generalization of the Weibull distribution capable of modeling a bathtub-shaped hazard rate function is defined and studied. The beauty and importance of this distribution lies in its ability to model monotone as well as non-monotone failure rates, which are quite common in lifetime problems and reliability. The new distribution has a number of well-known lifetime special sub-models, such as the Weibull, extreme value, exponentiated Weibull, generalized Rayleigh and modified Weibull distributions, among others. We derive two infinite sum representations for its moments. The density of the order statistics is obtained. The method of maximum likelihood is used for estimating the model parameters. Also, the observed information matrix is obtained. Two applications are presented to illustrate the proposed distribution. (C) 2008 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.

General results for the beta-modified Weibull distribution

We study in detail the so-called beta-modified Weibull distribution, motivated by the wide use of the Weibull distribution in practice, and also for the fact that the generalization provides a continuous crossover towards cases with different shapes. The new distribution is important since it contains as special sub-models some widely-known distributions, such as the generalized modified Weibull, beta Weibull, exponentiated Weibull, beta exponential, modified Weibull and Weibull distributions, among several others. It also provides more flexibility to analyse complex real data. Various mathematical properties of this distribution are derived, including its moments and moment generating function. We examine the asymptotic distributions of the extreme values. Explicit expressions are also derived for the chf, mean deviations, Bonferroni and Lorenz curves, reliability and entropies. The estimation of parameters is approached by two methods: moments and maximum likelihood. We compare by simulation the performances of the estimates from these methods. We obtain the expected information matrix. Two applications are presented to illustrate the proposed distribution.

The exponentiated generalized gamma distribution with application to lifetime data

A four-parameter extension of the generalized gamma distribution capable of modelling a bathtub-shaped hazard rate function is defined and studied. The beauty and importance of this distribution lies in its ability to model monotone and non-monotone failure rate functions, which are quite common in lifetime data analysis and reliability. The new distribution has a number of well-known lifetime special sub-models, such as the exponentiated Weibull, exponentiated generalized half-normal, exponentiated gamma and generalized Rayleigh, among others. We derive two infinite sum representations for its moments. We calculate the density of the order statistics and two expansions for their moments. The method of maximum likelihood is used for estimating the model parameters and the observed information matrix is obtained. Finally, a real data set from the medical area is analysed.

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.

Performance of a multi-level TDR system exposed to tropical field conditions - A time-space comparison with tensiometry

Time-domain reflectometry (TDR) is an important technique to obtain series of soil water content measurements in the field. Diode-segmented probes represent an improvement in TDR applicability, allowing measurements of the soil water content profile with a single probe. In this paper we explore an extensive soil water content dataset obtained by tensiometry and TDR from internal drainage experiments in two consecutive years in a tropical soil in Brazil. Comparisons between the variation patterns of the water content estimated by both methods exhibited evidences of deterioration of the TDR system during this two year period at field conditions. The results showed consistency in the variation pattern for the tensiometry data, whereas TDR estimates were inconsistent, with sensitivity decreasing over time. This suggests that difficulties may arise for the long-term use of this TDR system under tropical field conditions. (c) 2008 Elsevier B.V. All rights reserved.

An integrated molecular map of yellow passion fruit based on simultaneous maximum-likelihood estimation of linkage and linkage phases

The development of genetic maps for auto-incompatible species, such as the yellow passion fruit (Passiflora edulis Sims f.flavicarpa Deg.) is restricted due to the unfeasibility of obtaining traditional mapping populations based on inbred lines. For this reason, yellow passion fruit linkage maps were generally constructed using a strategy known as two-way pseudo-testeross, based on monoparental dominant markers segregating in a 1:1 fashion. Due to the lack of information from these markers in one of the parents, two individual (parental) maps were obtained. However, integration of these maps is essential, and biparental markers can be used for such an operation. The objective of our study was to construct an integrated molecular map for a full-sib population of yellow passion fruit combining different loci configuration generated from amplified fragment length polymorphisms (AFLPs) and microsatellite markers and using a novel approach based on simultaneous maximum-likelihood estimation of linkage and linkage phases, specially designed for outcrossing species. Of the total number of loci, approximate to 76%, 21%, 0.7%, and 2.3% did segregate in 1:1, 3:1, 1:2:1, and 1:1:1:1 ratios, respectively. Ten linkage groups (LGs) were established with a logarithm of the odds (LOD) score >= 5.0 assuming a recombination fraction : <= 0.35. On average, 24 markers were assigned per LG, representing a total map length of 1687 cM, with a marker density of 6.9 cM. No markers were placed as accessories on the map as was done with previously constructed individual maps.

Estimating Amino Acid Requirement of Brazilian Freshwater Fish from Muscle Amino Acid Profile

Information on nutritional requirement of some Brazilian farmed fish species, especially essential amino acids (EAA) requirements, is scarce. The estimation of amino acids requirements based on amino acid composition of fish is a fast and reliable alternative. Matrinxa, Brycon amazonicus, and curimbata, Prochilodus lineatus, are two important Brazilian fish with potential for aquaculture. The objective of the present study was to estimate amino acid requirements of these species and analyze similarities among amino acid composition of different fish species by cluster analysis. To estimate amino acid requirement, the following formula was used: amino acid requirement = [(amount of an individual amino acid in fish muscle tissue) x (average totalEAA requirement among channel catfish, Ictalurus punctatus, Nile tilapia, Oreochromis niloticus, and common carp, Cyprinus carpio)]/(average fish muscle totalEAA). Most values found lie within the range of requirements determined for other omnivorous fish species, in exception of leucine requirement estimated for both species, and arginine requirement estimated for matrinxa alone. Rather than writing off the need for regular dose-response assays under the ideal protein concept to determine EAA requirements of curimbata and matrinxa, results set solid base for the study of tropical species dietary amino acids requirements.

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.

Ovarian function in cows submitted to acute stress during proestrus

Studies in which ACTH was administrated in heifers after the occurrence of luteolysis showed a delay in the onset of estrus and the estrus duration was shortened. This study evaluated the effect of acute stress by road transportation on estrous behavior and ovulation, monitored by serum progesterone and cortisol concentrations in cows at the periovulatory period, using a crossover design. Eleven crossbred cows, divided into Control and Transport groups had their estrus cycle synchronized with GnRH, an intravaginal progesterone device, and cloprostenol. Thirty hours after withdrawal of the device, the animals of the Transport group were transported for 60 min by truck and those from the Control group remained at pasture. Ovarian ultrasound examination was performed every 12 h from device withdrawal until ovulation in every cow. From the day after removal of the device until ovulation estrous behavior was monitored 24 h a day. Blood samples for serum cortisol and progesterone concentrations were taken at -90, -60, 0, 30, 60 and 180 min in relation to the end of transportation. Transportation during the estrous period induced stress in cows as reflected by changes in serum concentrations of progesterone and cortisol. However, we did not detect impairment in estrus expression, estrus duration or ovulation (P>0.05). (C) 2010 Elsevier B.V. All rights reserved.

mRNA profile of Nellore calves after primary infection with Haemonchus placei

Haemonchus parasites are responsible for many losses in animal production. However, few studies are available, especially of zebu cattle. In this study, we investigated mRNA differences of immune response genes in naive Nellore calves infected with Haemonchus placei, relating these differences to patterns of cellular infiltrate. Calves were infected with 15,000 H. placei 13 larvae and after 7 days lymph node and abomasum tissues were collected. IL-2, IL-4, IL-8, IL-12, IL-13, IFN-gamma, MCP-1, lysozyme, pepsinogen and TNF-alpha genes were evaluated by qPCR. Mast cells, eosinophils and globular leukocytes were counted by abomasum histology. In the infected group, IL-4, IL-13 and TNF-alpha were up-regulated in the abomasal lymph node. In the abomasum, IL-13 increased and TNF-alpha was down-regulated (p < 0.05). No differences were detected for mast cells and eosinophil counts in abomasal tissue (p > 0.05). We conclude that for this infection time, there was Th2 polarization but that cellular infiltrate in abomasal tissue takes longer to develop. (C) 2010 Elsevier B.V. All rights reserved.

Nutrition in bean plants and anthracnose intensity in function of silicon and copper application.

Nutrition in bean plants and anthracnose intensity in function of silicon and copper application. The objective of this work was to evaluate the effect of calcium silicate and copper sulfate on anthracnose intensity and nutrition of bean plants. The experiment was conducted using an experimental design in randomized blocks following a 4 x 4 factorial arrangement , (four levels of calcium silicate and four levels of copper sulfate) and two additional treatments (plants without inoculation and plants sprinkled with Benomyl). Four evaluations of the incidence and severity of anthracnose were done, in addition to measuring, total leaf area. At the end of the evaluations, incidence: and data were integrated over time, obtaining the area under disease progress curve (AUDPC). Contents of N, P, K, Ca, Mg, B, Cu, Fe, Mn, Zn, Si and lignin were determined in the aerial Part. A linear decrease of the intensity AUDPC was observed with the increase of the doses of calcium silicate. The severity AUDPC was influenced by the doses of copper, obtaining a reduction of 35% on the higher dosage. The supply of silicon and copper altered the content of the K, mg, S, Zn, Ca and Si in the aerial part of the bean plants.