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.

Different Approaches for Modeling Grouped Survival Data: A Mango Tree Study

Interval-censored survival data, in which the event of interest is not observed exactly but is only known to occur within some time interval, occur very frequently. In some situations, event times might be censored into different, possibly overlapping intervals of variable widths; however, in other situations, information is available for all units at the same observed visit time. In the latter cases, interval-censored data are termed grouped survival data. Here we present alternative approaches for analyzing interval-censored data. We illustrate these techniques using a survival data set involving mango tree lifetimes. This study is an example of grouped survival data.

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.

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.

Evaluation of algorithms used to order markers on genetic maps

When building genetic maps, it is necessary to choose from several marker ordering algorithms and criteria, and the choice is not always simple. In this study, we evaluate the efficiency of algorithms try (TRY), seriation (SER), rapid chain delineation (RCD), recombination counting and ordering (RECORD) and unidirectional growth (UG), as well as the criteria PARF (product of adjacent recombination fractions), SARF (sum of adjacent recombination fractions), SALOD (sum of adjacent LOD scores) and LHMC (likelihood through hidden Markov chains), used with the RIPPLE algorithm for error verification, in the construction of genetic linkage maps. A linkage map of a hypothetical diploid and monoecious plant species was simulated containing one linkage group and 21 markers with fixed distance of 3 cM between them. In all, 700 F(2) populations were randomly simulated with and 400 individuals with different combinations of dominant and co-dominant markers, as well as 10 and 20% of missing data. The simulations showed that, in the presence of co-dominant markers only, any combination of algorithm and criteria may be used, even for a reduced population size. In the case of a smaller proportion of dominant markers, any of the algorithms and criteria (except SALOD) investigated may be used. In the presence of high proportions of dominant markers and smaller samples (around 100), the probability of repulsion linkage increases between them and, in this case, use of the algorithms TRY and SER associated to RIPPLE with criterion LHMC would provide better results. Heredity (2009) 103, 494-502; doi:10.1038/hdy.2009.96; published online 29 July 2009

Central and peripheral mediation of human force sensation following eccentric or concentric contractions

Fatigue was induced in the triceps brachii of the experimental arm by a regimen of either eccentric or concentric muscle actions. Estimates of force were assessed using a contralateral limb-matching procedure, in which target force levels (25 %, 50 % or 75 % of maximum) were defined by the unfatigued control arm. Maximum isometric force-generating capacity was reduced by 31 % immediately following eccentric contractions, and remained depressed at 24 (25 %) and 48 h (13 %) post-exercise. A less marked reduction (8.3 %) was observed immediately following concentric contractions. Those participants who performed prior eccentric contractions, consistently (at all force levels), and persistently (throughout the recovery period), overestimated the level of force applied by the experimental arm. In other words, they believed that they were generating more force than they actually achieved. When the forces applied by the experimental and the control arm, were each expressed as a proportion of the maximum force that could be attained at that time, the estimates matched extremely closely. This outcome is that which would be expected if the estimates of force were based on a sense of effort. Following eccentric exercise, the amplitude of the EMG activity recorded from the experimental arm was substantially greater than that recorded from the control arm. Cortically evoked potentials recorded from the triceps brachii (and extensor carpi radialis) of the experimental arm were also substantially larger than those elicited prior to exercise. The sense of effort was evidently not based upon a corollary of the central motor command. Rather, the relationship between the sense of effort and the motor command appears to have been altered as a result of the fatiguing eccentric contractions. It is proposed that the sense of effort is associated with activity in neural centres upstream of the motor cortex.

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.

Remotely sensed estimates of vegetation structural characteristics in restored wetlands, Southern California

Traditional field sampling approaches for ecological studies of restored habitat can only cover small areas in detail, con be time consuming, and are often invasive and destructive. Spatially extensive and non-invasive remotely sensed data can make field sampling more focused and efficient. The objective of this work was to investigate the feasibility and accuracy of hand-held and airborne remotely sensed data to estimate vegetation structural parameters for an indicator plant species in a restored wetland. High spatial resolution, digital, multispectral camera images were captured from an aircraft over Sweetwater Marsh (San Diego County, California) during each growing season between 1992-1996. Field data were collected concurrently, which included plant heights, proportional ground cover and canopy architecture type, and spectral radiometer measurements. Spartina foliosa (Pacific cordgrass) is the indicator species for the restoration monitoring. A conceptual model summarizing the controls on the spectral reflectance properties of Pacific cordgrass was established. Empirical models were developed relating the stem length, density, and canopy architecture of cordgrass to normalized-difference-vegetation-index values. The most promising results were obtained from empirical estimates of total ground cover using image data that had been stratified into high, middle, and low marsh zones. As part of on-going restoration monitoring activities, this model is being used to provide maps of estimated vegetation cover.

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.