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.

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.

Isolation and characterization of actinobacteria ectosymbionts from Acromyrmex subterraneus brunneus (Hymenoptera, Formicidae)

The ectosymbiont actinobacterium Pseudonocardia was isolated from the integument of Acromyrmex leaf-cutter ants and seems to play a crucial role in maintaining asepsis of the nest. Currently, there has been an intensive search for Pseudonocardia associated with several attine species, but few studies have indicated that other actinobacteria may be associated with these ants as well. We therefore characterized the culturable actinobacteria community associated with the integument of the fungus-growing ant Acromyrmex subterraneus brunneus Forel, 1893 (Hymenoptera: Formicidae). Ectosymbionts were isolated using four different media and characterized by morphological and molecular (16S rDNA) methods. A total of 20 strains were isolated, of which 17 were characterized as Streptomyces spp., and one isolate each as Pseudonocardia, Kitassatospora and Propionicimonas. Unlike other Acromyrmex species, A. subterraneus brunneus is associated with a diversity of actinobacteria. Even though Pseudonocardia is present on this leaf-cutting ant`s integument, the number and diversity of Streptomyces spp. found differs from those of previous studies with other attine ants and suggest that different culturing approaches are needed to characterize the true diversity of microbes colonizing the integument of attine ants. Moreover, understanding the diversity of the culturable actinobacteria associated with A. subterraneus brunneus should increase our knowledge of the evolutionary relationship of this intricate symbiotic association. (C) 2010 Elsevier GmbH. All rights reserved.

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.

Effects of light intensity and dilution rate on the semicontinuous cultivation of Arthrospira (Spirulina) platensis. A kinetic Monod-type approach

Semicontinuous cultures were carried out at different dilution rates (D) and light intensities (I) to determine the maximum productivity of Arthrospira platensis cultivated in helicoidal photobioreactor up to the achievement of pseudo-steady-state conditions. At I = 108 mu mol photons m(-2) s(-1), the semicontinuous regime ensured the highest values of maximum cell concentration (X(m) = 5772 +/- 113 mg L(-1)) and productivity (P(XS) = 1319 +/- 25 mg L(-1) d(-1)) at the lowest (D = 0.1 day(-1)) and the highest (D = 0.3 day(-1)) dilution rates, respectively. A kinetic model derived from that of Monod was proposed to determine the relationship between the product of light intensity to dilution rate (ID) and the cell productivity, which were shown to exert a combined influence on this parameter. This result put into evidence that pseudosteady-state conditions could be modified according to circumstances, conveniently varying one or other of the two independent variables. (C) 2010 Elsevier Ltd. All rights reserved.

Adsorption of Ni(2+), Zn(2+) and Pb(2+) onto dry biomass of Arthrospira (Spirulina) platensis and Chlorella vulgaris. I. Single metal systems

Adsorption of Ni(2+), Zn(2+) or Pb(2+) by dry biomass of Arthrospira (Spirulina) platensis and Chlorella vulgaris was studied as a function of contact time and initial metal concentration. The zero point of charge calculated for these biosorbents (pH(zpc) 4.0 and 3.4, respectively) and additional pH tests suggested the use of pH in the range 5.0-5.5 for the experiments. The equilibrium isotherms were evaluated in terms of maximum sorption capacity and sorption affinity. The pseudo first and second order kinetic models were considered to interpret the experimental data, and the latter best described the adsorption system. Both the Freundlich and Langmuir models were shown to well describe the sorption isotherms, thus suggesting an intermediate mono/multilayer sorption mechanism. Compared to A. platensis (q(e) = 0.354, 0.495 and 0.508 mmol g(-1) for Ni(2+), Pb(2)+ and Zn(2+), respectively), C. vulgaris behaved as a better biosorbent because of higher equilibrium sorption capacity (q(e) = 0.499, 0.634 and 0.664 mmol g(-1), respectively). The removal efficiency decreased with increasing metal concentration, pointing out a passive adsorption process involving the active sites on the surface of the biomasses. The FT-IR spectroscopy evidenced that ions removal occurred mainly by interaction between metal and carboxylate groups present on both the cell walls. (C) 2011 Elsevier B.V. All rights reserved.

Positive selection on an acrosomal sperm protein, M7 lysin, in three species of the mussel genus Mytilus

Marine invertebrate sperm proteins are particularly interesting because they are characterized by positive selection and are likely to be involved in prezyogotic isolation and, thus, speciation. Here, we present the first survey of inter and intraspecific variation of a bivalve sperm protein among a group of species that regularly hybridize in nature. M7 lysin is found in sperm acrosomes of mussels and dissolves the egg vitelline coat, permitting fertilization. We sequenced multiple alleles of the mature protein-coding region of M7 lysin from allopatric populations of mussels in the Mytilus edulis species group (M. edulis, M. galloprovincialis, and M. trossulus). A significant McDonald-Kreitman test showed an excess of fixed amino acid replacing substitutions between species, consistent with positive selection. In addition, Kolmogorov-Smirnov tests showed significant heterogeneity in polymorphism to divergence ratios for both synonymous variation and combined synonymous and non-synonymous variation within M. galloprovincialis. These results indicate that there has been adaptive evolution at M7 lysin and, furthermore, shows that positive selection on sperm proteins can occur even when post-zygotic reproductive isolation is incomplete.

Sampling phylogenetic tree space with the generalized Gibbs sampler

The generalized Gibbs sampler (GGS) is a recently developed Markov chain Monte Carlo (MCMC) technique that enables Gibbs-like sampling of state spaces that lack a convenient representation in terms of a fixed coordinate system. This paper describes a new sampler, called the tree sampler, which uses the GGS to sample from a state space consisting of phylogenetic trees. The tree sampler is useful for a wide range of phylogenetic applications, including Bayesian, maximum likelihood, and maximum parsimony methods. A fast new algorithm to search for a maximum parsimony phylogeny is presented, using the tree sampler in the context of simulated annealing. The mathematics underlying the algorithm is explained and its time complexity is analyzed. The method is tested on two large data sets consisting of 123 sequences and 500 sequences, respectively. The new algorithm is shown to compare very favorably in terms of speed and accuracy to the program DNAPARS from the PHYLIP package.

On the mechanism of microbial cell disruption in high-pressure homogenisation

The tensions produced in the wall of a rigid, thin-walled, liquid-filled sphere as it moves with an axisymmetric straining flow are examined. This problem has not been previously addressed. A generalised correlation for the maximum wall tension, expressed in dimensionless form as a Weber number (We), is developed in terms of the acceleration number (Ac) and Reynolds number (Re) of the straining flow. At low Reynolds number We is dominated by viscous forces, while inertial forces due to internal pressure gradients caused by sphere acceleration dominate at higher Re. The generalised correlation has been used to examine the case of a typical yeast cell (a thin-walled, liquid-filled sphere) passing through a typical high-pressure homogeniser (a straining-flow device). At 56 MPa homogenising pressure, a 6 mu m yeast cell experiences tensions in the inertially dominated regime (Re = 100). The correlation gives We = 0.206, corresponding to a maximum wall tension of 8 Nm(-1). This is equivalent to an applied compressive force of 150 mu N and compares favourably with the force required to break yeast cells under compressive micromanipulation (40-90 mu N). Inertial forces may therefore be an important and previously unrecognised. mechanism of microbial cell disruption during high-pressure homogenisation. Further work is required to examine the likelihood of cell deformation in the high-strain-rate short-residence-time environment of the homogeniser, and the effect that such deformation may have on the contribution of inertial forces to disruption. (C) 1998 Published by Elsevier Science Ltd. All rights reserved.

Distinguishing species and populations of Rhipicephaline ticks with its 2 ribosomal RNA

Most populations and some species of ticks of the genera Boophilus (5 spp.) and Rhipicephalus (ca. 75 spp.) cannot be distinguished phenotypically. Moreover, there is doubt about the validity of species in these genera. I studied the entire second internal transcribed spacer (ITS 2) rRNA of 16 populations of rhipicephaline ticks to address these problems: Boophilus,microplus from Australia, Kenya, South Africa and Brazil (4 populations); Boophilus decoloratus from Kenya; Rhipicephalus appendiculatus from Kenya, Zimbabwe and Zambia (7 populations); Rhipicephalus zambesiensis from Zimbabwe (3 populations); and Rhipicephalus evertsi from Kenya. Each of the 16 populations had a unique ITS 2, but most of the nucleotide variation occurred among species and genera. ITS 2 rRNA can be used to distinguish the populations and species of Boophilus and Rhipicephalus studied here. Little support was found for the hypothesis that B. microplus from Australia and South Africa are different species. ITS 2 appears useful for phylogenetic inference in the Rhipicephalinae because in genetic distance, maximum likelihood, and maximum parsimony analyses, most branches leading to species had >95% bootstrap support. Rhipicephalus appendiculatus and R, zambeziensis are closely related, yet their ITS 2 sequences could be distinguished unambiguously. This lends weight to a previous proposal that Rhipicephalus sanguineus and Rhipicephalus turanicus, and Rhipicephalus pumlilio and Rhipicephalus camicasi, respectively, are conspecific, because each of these pairs of species had identical sequences for ca. 250 bp of ITS 2 rRNA.

On modifications to the long-term survival mixture model in the presence of competing risks

A mixture model for long-term survivors has been adopted in various fields such as biostatistics and criminology where some individuals may never experience the type of failure under study. It is directly applicable in situations where the only information available from follow-up on individuals who will never experience this type of failure is in the form of censored observations. In this paper, we consider a modification to the model so that it still applies in the case where during the follow-up period it becomes known that an individual will never experience failure from the cause of interest. Unless a model allows for this additional information, a consistent survival analysis will not be obtained. A partial maximum likelihood (ML) approach is proposed that preserves the simplicity of the long-term survival mixture model and provides consistent estimators of the quantities of interest. Some simulation experiments are performed to assess the efficiency of the partial ML approach relative to the full ML approach for survival in the presence of competing risks.

Hemichordates and deuterostome evolution: robust molecular phylogenetic support for a hemichordate plus echinoderm clade

Hemichordates were traditionally allied to the chordates, but recent molecular analyses have suggested that hemichordates are a sister group to the echinoderms, a relationship that has important consequences for the interpretation of the evolution of deuterostome body plans. However, the molecular phylogenetic analyses to date have not provided robust support for the hemichordate + echinoderm clade. We use a maximum likelihood framework, including the parametric bootstrap, to reanalyze DNA data from complete mitochondrial genomes and nuclear 18S rRNA. This approach provides the first statistically significant support for the hemichordate + echinoderm clade from molecular data. This grouping implies that the ancestral deuterostome had features that included an adult with a pharynx and a dorsal nerve cord and an indirectly developing dipleurula-like larva.