Cladistic analysis of continuous modularized traits provides phylogenetic signals in Homo evolution

Evolutionary novelties in the skeleton are usually expressed as changes in the timing of growth of features intrinsically integrated at different hierarchical levels of development(1). As a consequence, most of the shape- traits observed across species do vary quantitatively rather than qualitatively(2), in a multivariate space(3) and in a modularized way(4,5). Because most phylogenetic analyses normally use discrete, hypothetically independent characters(6), previous attempts have disregarded the phylogenetic signals potentially enclosed in the shape of morphological structures. When analysing low taxonomic levels, where most variation is quantitative in nature, solving basic requirements like the choice of characters and the capacity of using continuous, integrated traits is of crucial importance in recovering wider phylogenetic information. This is particularly relevant when analysing extinct lineages, where available data are limited to fossilized structures. Here we show that when continuous, multivariant and modularized characters are treated as such, cladistic analysis successfully solves relationships among main Homo taxa. Our attempt is based on a combination of cladistics, evolutionary- development- derived selection of characters, and geometric morphometrics methods. In contrast with previous cladistic analyses of hominid phylogeny, our method accounts for the quantitative nature of the traits, and respects their morphological integration patterns. Because complex phenotypes are observable across different taxonomic groups and are potentially informative about phylogenetic relationships, future analyses should point strongly to the incorporation of these types of trait.

On the utility of female genitalia in the systematics of Curculionidae: examples from Entimini and Hylobiina (Coleoptera : Curculionidae)

The utility of the female genitalia and associated sclerites (tergite and sternite VIII) in the systematics of the Curculionidae is discussed. Examples from the tribe Entimini (Entiminae) and subtribe Hylobiina (Molytinae: Hylobiini) are given. The female characters prove to be informative for establishing the phylogenetic relationships among genera of Entimini. They are essential in determining species groups within the genus Amiticus Pascoe, Hylobiina.

A novel locus for split-hand/foot malformation associated with tibial hemimelia (SHFLD syndrome) maps to chromosome region 17p13.1-17p13.3

Split-hand/foot malformation (SHFM) associated with aplasia of long bones, SHFLD syndrome or Tibial hemimelia-ectrodactyly syndrome is a rare condition with autosomal dominant inheritance, reduced penetrance and an incidence estimated to be about 1 in 1,000,000 liveborns. To date, three chromosomal regions have been reported as strong candidates for harboring SHFLD syndrome genes: 1q42.2-q43, 6q14.1 and 2q14.2. We characterized the phenotype of nine affected individuals from a large family with the aim of mapping the causative gene. Among the nine affected patients, four had only SHFM of the hands and no tibial defects, three had both defects and two had only unilateral tibial hemimelia. In keeping with previous publications of this and other families, there was clear evidence of both variable expression and incomplete penetrance, the latter bearing hallmarks of anticipation. Segregation analysis and multipoint Lod scores calculations (maximum Lod score of 5.03 using the LINKMAP software) using all potentially informative family members, both affected and unaffected, identified the chromosomal region 17p13.1-17p13.3 as the best and only candidate for harboring a novel mutated gene responsible for the syndrome in this family. The candidate gene CRK located within this region was sequenced but no pathogenic mutation was detected.

Animal models for genetic neuromuscular diseases

The neuromuscular disorders are a heterogeneous group of genetic diseases, caused by mutations in genes coding sarcolemmal, sarcomeric, and citosolic muscle proteins. Deficiencies or loss of function of these proteins leads to variable degree of progressive loss of motor ability. Several animal models, manifesting phenotypes observed in neuromuscular diseases, have been identified in nature or generated in laboratory. These models generally present physiological alterations observed in human patients and can be used as important tools for genetic, clinic, and histopathological studies. The mdx mouse is the most widely used animal model for Duchenne muscular dystrophy (DMD). Although it is a good genetic and biochemical model, presenting total deficiency of the protein dystrophin in the muscle, this mouse is not useful for clinical trials because of its very mild phenotype. The canine golden retriever MD model represents a more clinically similar model of DMD due to its larger size and significant muscle weakness. Autosomal recessive limb-girdle MD forms models include the SJL/J mice, which develop a spontaneous myopathy resulting from a mutation in the Dysferlin gene, being a model for LGMD2B. For the human sarcoglycanopahties (SG), the BIO14.6 hamster is the spontaneous animal model for delta-SG deficiency, whereas some canine models with deficiency of SG proteins have also been identified. More recently, using the homologous recombination technique in embryonic stem cell, several mouse models have been developed with null mutations in each one of the four SG genes. All sarcoglycan-null animals display a progressive muscular dystrophy of variable severity and share the property of a significant secondary reduction in the expression of the other members of the sarcoglycan subcomplex and other components of the Dystrophin-glycoprotein complex. Mouse models for congenital MD include the dy/dy (dystrophia-muscularis) mouse and the allelic mutant dy(2J)/dy(2J) mouse, both presenting significant reduction of alpha 2-laminin in the muscle and a severe phenotype. The myodystrophy mouse (Large(myd)) harbors a mutation in the glycosyltransferase Large, which leads to altered glycosylation of alpha-DG, and also a severe phenotype. Other informative models for muscle proteins include the knockout mouse for myostatin, which demonstrated that this protein is a negative regulator of muscle growth. Additionally, the stress syndrome in pigs, caused by mutations in the porcine RYR1 gene, helped to localize the gene causing malignant hypertermia and Central Core myopathy in humans. The study of animal models for genetic diseases, in spite of the existence of differences in some phenotypes, can provide important clues to the understanding of the pathogenesis of these disorders and are also very valuable for testing strategies for therapeutic approaches.

Estimating the number of ozone peaks in Mexico City using a non-homogeneous Poisson model

In this paper, we consider the problem of estimating the number of times an air quality standard is exceeded in a given period of time. A non-homogeneous Poisson model is proposed to analyse this issue. The rate at which the Poisson events occur is given by a rate function lambda(t), t >= 0. This rate function also depends on some parameters that need to be estimated. Two forms of lambda(t), t >= 0 are considered. One of them is of the Weibull form and the other is of the exponentiated-Weibull form. The parameters estimation is made using a Bayesian formulation based on the Gibbs sampling algorithm. The assignation of the prior distributions for the parameters is made in two stages. In the first stage, non-informative prior distributions are considered. Using the information provided by the first stage, more informative prior distributions are used in the second one. The theoretical development is applied to data provided by the monitoring network of Mexico City. The rate function that best fit the data varies according to the region of the city and/or threshold that is considered. In some cases the best fit is the Weibull form and in other cases the best option is the exponentiated-Weibull. Copyright (C) 2007 John Wiley & Sons, Ltd.

A bivariate regression model for matched paired survival data: local influence and residual analysis

The use of bivariate distributions plays a fundamental role in survival and reliability studies. In this paper, we consider a location scale model for bivariate survival times based on the proposal of a copula to model the dependence of bivariate survival data. For the proposed model, we consider inferential procedures based on maximum likelihood. Gains in efficiency from bivariate models are also examined in the censored data setting. For different parameter settings, sample sizes and censoring percentages, various simulation studies are performed and compared to the performance of the bivariate regression model for matched paired survival data. Sensitivity analysis methods such as local and total influence are presented and derived under three perturbation schemes. The martingale marginal and the deviance marginal residual measures are used to check the adequacy of the model. Furthermore, we propose a new measure which we call modified deviance component residual. The methodology in the paper is illustrated on a lifetime data set for kidney patients.

The log-exponentiated Weibull regression model for interval-censored data

In interval-censored survival data, the event of interest is not observed exactly but is only known to occur within some time interval. Such data appear very frequently. In this paper, we are concerned only with parametric forms, and so a location-scale regression model based on the exponentiated Weibull distribution is proposed for modeling interval-censored data. We show that the proposed log-exponentiated Weibull regression model for interval-censored data represents a parametric family of models that include other regression models that are broadly used in lifetime data analysis. Assuming the use of interval-censored data, we employ a frequentist analysis, a jackknife estimator, a parametric bootstrap and a Bayesian analysis 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. Furthermore, for different parameter settings, sample sizes and censoring percentages, various 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 straightforwardly extended to a modified deviance residual in log-exponentiated Weibull regression models for interval-censored data. (C) 2009 Elsevier B.V. All rights reserved.

The complementary exponential geometric distribution: Model, properties, and a comparison with its counterpart

In this paper, we proposed a new two-parameter lifetime distribution with increasing failure rate, the complementary exponential geometric distribution, which is complementary to the exponential geometric model proposed by Adamidis and Loukas (1998). The new distribution arises on a latent complementary risks scenario, in which the lifetime associated with a particular risk is not observable; rather, we observe only the maximum lifetime value among all risks. The properties of the proposed distribution are discussed, including a formal proof of its probability density function and explicit algebraic formulas for its reliability and failure rate functions, moments, including the mean and variance, variation coefficient, and modal value. The parameter estimation is based on the usual maximum likelihood approach. We report the results of a misspecification simulation study performed in order to assess the extent of misspecification errors when testing the exponential geometric distribution against our complementary one in the presence of different sample size and censoring percentage. The methodology is illustrated on four real datasets; we also make a comparison between both modeling approaches. (C) 2011 Elsevier B.V. All rights reserved.

Log-Burr XII regression models with censored data

In survival analysis applications, the failure rate function may frequently present a unimodal shape. In such case, the log-normal or log-logistic distributions are used. In this paper, we shall be concerned only with parametric forms, so a location-scale regression model based on the Burr XII distribution is proposed for modeling data with a unimodal failure rate function as an alternative to the log-logistic regression model. Assuming censored data, we consider a classic analysis, a Bayesian analysis and a jackknife estimator for the parameters of the proposed model. For different parameter settings, sample sizes and censoring percentages, various simulation studies are performed and compared to the performance of the log-logistic and log-Burr XII regression models. Besides, we use sensitivity analysis to detect influential or outlying observations, and residual analysis is used to check the assumptions in the model. Finally, we analyze a real data set under log-Buff XII regression models. (C) 2008 Published by Elsevier B.V.

Deviance residuals in generalised log-gamma regression models with censored observations

In this article, we compare three residuals based on the deviance component in generalised log-gamma regression models with censored observations. For different parameter settings, sample sizes and censoring percentages, various simulation studies are performed and the empirical distribution of each residual is displayed and compared with the standard normal distribution. For all cases studied, the empirical distributions of the proposed residuals are in general symmetric around zero, but only a martingale-type residual presented negligible kurtosis for the majority of the cases studied. These studies suggest that the residual analysis usually performed in normal linear regression models can be straightforwardly extended for the martingale-type residual in generalised log-gamma regression models with censored data. A lifetime data set is analysed under log-gamma regression models and a model checking based on the martingale-type residual is performed.

Improved likelihood inference for the shape parameter in Weibull regression

We obtain adjustments to the profile likelihood function in Weibull regression models with and without censoring. Specifically, we consider two different modified profile likelihoods: (i) the one proposed by Cox and Reid [Cox, D.R. and Reid, N., 1987, Parameter orthogonality and approximate conditional inference. Journal of the Royal Statistical Society B, 49, 1-39.], and (ii) an approximation to the one proposed by Barndorff-Nielsen [Barndorff-Nielsen, O.E., 1983, On a formula for the distribution of the maximum likelihood estimator. Biometrika, 70, 343-365.], the approximation having been obtained using the results by Fraser and Reid [Fraser, D.A.S. and Reid, N., 1995, Ancillaries and third-order significance. Utilitas Mathematica, 47, 33-53.] and by Fraser et al. [Fraser, D.A.S., Reid, N. and Wu, J., 1999, A simple formula for tail probabilities for frequentist and Bayesian inference. Biometrika, 86, 655-661.]. We focus on point estimation and likelihood ratio tests on the shape parameter in the class of Weibull regression models. We derive some distributional properties of the different maximum likelihood estimators and likelihood ratio tests. The numerical evidence presented in the paper favors the approximation to Barndorff-Nielsen`s adjustment.

Inferential Implications of Over-Parametrization: A Case Study in Incomplete Categorical Data

P>In the context of either Bayesian or classical sensitivity analyses of over-parametrized models for incomplete categorical data, it is well known that prior-dependence on posterior inferences of nonidentifiable parameters or that too parsimonious over-parametrized models may lead to erroneous conclusions. Nevertheless, some authors either pay no attention to which parameters are nonidentifiable or do not appropriately account for possible prior-dependence. We review the literature on this topic and consider simple examples to emphasize that in both inferential frameworks, the subjective components can influence results in nontrivial ways, irrespectively of the sample size. Specifically, we show that prior distributions commonly regarded as slightly informative or noninformative may actually be too informative for nonidentifiable parameters, and that the choice of over-parametrized models may drastically impact the results, suggesting that a careful examination of their effects should be considered before drawing conclusions.Resume Que ce soit dans un cadre Bayesien ou classique, il est bien connu que la surparametrisation, dans les modeles pour donnees categorielles incompletes, peut conduire a des conclusions erronees. Cependant, certains auteurs persistent a negliger les problemes lies a la presence de parametres non identifies. Nous passons en revue la litterature dans ce domaine, et considerons quelques exemples surparametres simples dans lesquels les elements subjectifs influencent de facon non negligeable les resultats, independamment de la taille des echantillons. Plus precisement, nous montrons comment des a priori consideres comme peu ou non-informatifs peuvent se reveler extremement informatifs en ce qui concerne les parametres non identifies, et que le recours a des modeles surparametres peut avoir sur les conclusions finales un impact considerable. Ceci suggere un examen tres attentif de l`impact potentiel des a priori.

Bayesian estimation of regression parameters in elliptical measurement error models

The main object of this paper is to discuss the Bayes estimation of the regression coefficients in the elliptically distributed simple regression model with measurement errors. The posterior distribution for the line parameters is obtained in a closed form, considering the following: the ratio of the error variances is known, informative prior distribution for the error variance, and non-informative prior distributions for the regression coefficients and for the incidental parameters. We proved that the posterior distribution of the regression coefficients has at most two real modes. Situations with a single mode are more likely than those with two modes, especially in large samples. The precision of the modal estimators is studied by deriving the Hessian matrix, which although complicated can be computed numerically. The posterior mean is estimated by using the Gibbs sampling algorithm and approximations by normal distributions. The results are applied to a real data set and connections with results in the literature are reported. (C) 2011 Elsevier B.V. All rights reserved.

Testing hypotheses in the Birnbaum-Saunders distribution under type-II censored samples

The two-parameter Birnbaum-Saunders distribution has been used successfully to model fatigue failure times. Although censoring is typical in reliability and survival studies, little work has been published on the analysis of censored data for this distribution. In this paper, we address the issue of performing testing inference on the two parameters of the Birnbaum-Saunders distribution under type-II right censored samples. The likelihood ratio statistic and a recently proposed statistic, the gradient statistic, provide a convenient framework for statistical inference in such a case, since they do not require to obtain, estimate or invert an information matrix, which is an advantage in problems involving censored data. An extensive Monte Carlo simulation study is carried out in order to investigate and compare the finite sample performance of the likelihood ratio and the gradient tests. Our numerical results show evidence that the gradient test should be preferred. Further, we also consider the generalized Birnbaum-Saunders distribution under type-II right censored samples and present some Monte Carlo simulations for testing the parameters in this class of models using the likelihood ratio and gradient tests. Three empirical applications are presented. (C) 2011 Elsevier B.V. All rights reserved.

Distribution of ABCB1 polymorphisms among Brazilians: impact of population admixture

Introduction: Interethnic admixture is a source of cryptic population structure that may lead to spurious genotype-phenotype associations in pharmacogenomic studies. We studied the impact of population stratification on the distribution of ABCB1 polymorphisms (1236C > T, 2677G > T/A and 3435C > T) among Brazilians, a highly admixed population with Amerindian, European and African ancestral roots. Methods: Individual DNA from 320 healthy adults was genotyped with a panel of ancestry informative markers, and the proportions of African component of ancestry (ACA) were estimated. ABCB1 genotypes were determined by the single base extension/termination method. We describe the association between ABCB1 polymorphisms and ACA by fitting a linear proportional odds logistic regression model to the data. Results: The distribution of the ABCB1 2677G > T/A and 3435C > T, but not the 1236C > T, SNPs displayed a significant trend for decreasing frequency of the T alleles and TT genotypes from White to Intermediate to Black individuals. The same trend was observed in the frequency of the T/nonG/T haplotype at the 1236, 2677 and 3435 loci. When the population sample was proportioned in quartiles, according to the individual ACA estimates, the frequency of the T allele and TT genotype at each locus declined progressively from the lowest (< 0.25 ACA) to the highest (> 0.75 ACA) quartile. Linear proportional odds logistic regression analysis confirmed that the odds of having the T allele at each locus decreases in a continuous manner with the increase of the ACA, throughout the ACA range (0.13-0.94) observed in the overall population sample. A significant association was also detected between the individual ACA estimates and the presence of the T/nonG/T haplotype in the overall population. Conclusion: Self-identification according to the racial/color categories proposed by the Brazilian Census is insufficient to properly control for population stratification in pharmacogenomic studies of ABCB1.