Haplotype distribution of five nuclear genes based on network genealogies and Bayesian inference indicates that Trypanosoma cruzi hybrid strains are polyphyletic

Chagas disease is still a major public health problem in Latin America. Its causative agent, Trypanosoma cruzi, can be typed into three major groups, T. cruzi I, T. cruzi II and hybrids. These groups each have specific genetic characteristics and epidemiological distributions. Several highly virulent strains are found in the hybrid group; their origin is still a matter of debate. The null hypothesis is that the hybrids are of polyphyletic origin, evolving independently from various hybridization events. The alternative hypothesis is that all extant hybrid strains originated from a single hybridization event. We sequenced both alleles of genes encoding EF-1 alpha, actin and SSU rDNA of 26 T. cruzi strains and DHFR-TS and TR of 12 strains. This information was used for network genealogy analysis and Bayesian phylogenies. We found T. cruzi I and T. cruzi II to be monophyletic and that all hybrids had different combinations of T. cruzi I and T. cruzi II haplotypes plus hybrid-specific haplotypes. Bootstrap values (networks) and posterior probabilities (Bayesian phylogenies) of clades supporting the monophyly of hybrids were far below the 95% confidence interval, indicating that the hybrid group is polyphyletic. We hypothesize that T. cruzi I and T. cruzi II are two different species and that the hybrids are extant representatives of independent events of genome hybridization, which sporadically have sufficient fitness to impact on the epidemiology of Chagas disease.

Molecular characterization of the Hepatitis B virus genotypes in Colombia: A Bayesian inference on the genotype F

Hepatitis B is a worldwide health problem affecting about 2 billion people and more than 350 million are chronic carriers of the virus. Nine HBV genotypes (A to I) have been described. The geographical distribution of HBV genotypes is not completely understood due to the limited number of samples from some parts of the world. One such example is Colombia, in which few studies have described the HBV genotypes. In this study, we characterized HBV genotypes in 143 HBsAg-positive volunteer blood donors from Colombia. A fragment of 1306 bp partially comprising HBsAg and the DNA polymerase coding regions (S/POL) was amplified and sequenced. Bayesian phylogenetic analyses were conducted using the Markov Chain Monte Carlo (MCMC) approach to obtain the maximum clade credibility (MCC) tree using BEAST v.1.5.3. Of all samples, 68 were positive and 52 were successfully sequenced. Genotype F was the most prevalent in this population (77%) - subgenotypes F3 (75%) and Fib (2%). Genotype G (7.7%) and subgenotype A2 (15.3%) were also found. Genotype G sequence analysis suggests distinct introductions of this genotype in the country. Furthermore, we estimated the time of the most recent common ancestor (TMRCA) for each HBV/F subgenotype and also for Colombian F3 sequences using two different datasets: (i) 77 sequences comprising 1306 bp of S/POL region and (ii) 283 sequences comprising 681 bp of S/POL region. We also used two other previously estimated evolutionary rates: (i) 2.60 x 10(-4) s/s/y and (ii) 1.5 x 10(-5) s/s/y. Here we report the HBV genotypes circulating in Colombia and estimated the TMRCA for the four different subgenotypes of genotype F. (C) 2010 Elsevier B.V. All rights reserved.

Bayesian inference for a skew-normal IRT model under the centred parameterization

Item response theory (IRT) comprises a set of statistical models which are useful in many fields, especially when there is interest in studying latent variables. These latent variables are directly considered in the Item Response Models (IRM) and they are usually called latent traits. A usual assumption for parameter estimation of the IRM, considering one group of examinees, is to assume that the latent traits are random variables which follow a standard normal distribution. However, many works suggest that this assumption does not apply in many cases. Furthermore, when this assumption does not hold, the parameter estimates tend to be biased and misleading inference can be obtained. Therefore, it is important to model the distribution of the latent traits properly. In this paper we present an alternative latent traits modeling based on the so-called skew-normal distribution; see Genton (2004). We used the centred parameterization, which was proposed by Azzalini (1985). This approach ensures the model identifiability as pointed out by Azevedo et al. (2009b). Also, a Metropolis Hastings within Gibbs sampling (MHWGS) algorithm was built for parameter estimation by using an augmented data approach. A simulation study was performed in order to assess the parameter recovery in the proposed model and the estimation method, and the effect of the asymmetry level of the latent traits distribution on the parameter estimation. Also, a comparison of our approach with other estimation methods (which consider the assumption of symmetric normality for the latent traits distribution) was considered. The results indicated that our proposed algorithm recovers properly all parameters. Specifically, the greater the asymmetry level, the better the performance of our approach compared with other approaches, mainly in the presence of small sample sizes (number of examinees). Furthermore, we analyzed a real data set which presents indication of asymmetry concerning the latent traits distribution. The results obtained by using our approach confirmed the presence of strong negative asymmetry of the latent traits distribution. (C) 2010 Elsevier B.V. All rights reserved.

A synergistic relationship of elevated low-density lipoprotein cholesterol levels and systolic blood pressure with coronary artery calcification

We sought to evaluate this ""response-to-injury"" hypothesis of atherosclerosis by studying the interaction between systolic blood pressure (SBP) and LDL-cholesterol (LDL-C) in predicting the presence of coronary artery calcification (CAC) in asymptomatic men. We Studied 526 men (46 +/- 7 years of age) referred for electron-beam tomography (EBT) exam. The prevalence of CAC was determined across LDL-C tertiles (low: <115 mg/dl; middle: 115-139 mg/dl high: >= 140 mg/dl) within tertiles of SBP (low: <121 mmHg; middle: 121-130 mmHg; high: >= 131 mmHg). CAC was found in 220 (42%) men. There was no linear trend in the presence of CAC across LDL-C tertiles in the low (p = 0.6 for trend) and middle (p = 0.3 for trend) SBP tertile groups, respectively. In contrast, there was a significant trend for increasing CAC with increasing LDL-C (1st: 44%; 2nd: 49%; 3rd: 83%; p < 0.0001 for trend) in the high SBP tertile group. In multivariate logistic analyses (adjusting for age, smoking, triglyceride levels, HDL-cholesterol levels, body mass index, and fasting glucose levels), the odds ratio for any CAC associated with increasing LDL-C was significantly higher in those with highest SBP levels, whereas no such relationship was observed among men with SBP in the lower two tertiles. An interaction term (LDL-C x SBP) incorporated in the multivariate analyses was statistically significant (p = 0.038). The finding of an interaction between SBP and LDL-C relation to CAC in asymptomatic men support the response-to-injury model of atherogenesis. (C) 2007 Elsevier Ireland Ltd. All rights reserved.

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.

Traffic-Related Air Pollution and Perinatal Mortality: A Case-Control Study

BACKGROUND: Ambient levels of air pollution may affect the health of children, as indicated by studies of infant and perinatal mortality. Scientific evidence has also correlated low birth weight and preterm birth, which are important determinants of perinatal death, with air pollution. However, most of these studies used ambient concentrations measured at monitoring sites, which may not consider differential exposure to pollutants found at elevated concentrations near heavy-traffic roadways. OBJECTIVES: Our goal was to examine the association between traffic-related pollution and perinatal mortality. METHODS: We used the information collected for a case-control study conducted in 14 districts in the City of Sao Paulo, Brazil, regarding risk factors for perinatal deaths. We geocoded the residential addresses of cases (fetal and early neonatal deaths) and controls (children who survived the 28th day of life) and calculated a distance-weighted traffic density (DWTD) measure considering all roads contained in a buffer surrounding these homes. RESULTS: Logistic regression revealed a gradient of increasing risk of early neonatal death with higher exposure to traffic-related air pollution. Mothers exposed to the highest quartile of the DWTD compared with those less exposed exhibited approximately 50% increased risk (adjusted odds ratio = 1.47; 95% confidence interval, 0.67-3.19). Associations for fetal mortality were less consistent. CONCLUSIONS: These results suggest that motor vehicle exhaust exposures may be a risk factor for perinatal mortality.

Disequilibrium Coefficient: A Bayesian Perspective

Hardy-Weinberg Equilibrium (HWE) is an important genetic property that populations should have whenever they are not observing adverse situations as complete lack of panmixia, excess of mutations, excess of selection pressure, etc. HWE for decades has been evaluated; both frequentist and Bayesian methods are in use today. While historically the HWE formula was developed to examine the transmission of alleles in a population from one generation to the next, use of HWE concepts has expanded in human diseases studies to detect genotyping error and disease susceptibility (association); Ryckman and Williams (2008). Most analyses focus on trying to answer the question of whether a population is in HWE. They do not try to quantify how far from the equilibrium the population is. In this paper, we propose the use of a simple disequilibrium coefficient to a locus with two alleles. Based on the posterior density of this disequilibrium coefficient, we show how one can conduct a Bayesian analysis to verify how far from HWE a population is. There are other coefficients introduced in the literature and the advantage of the one introduced in this paper is the fact that, just like the standard correlation coefficients, its range is bounded and it is symmetric around zero (equilibrium) when comparing the positive and the negative values. To test the hypothesis of equilibrium, we use a simple Bayesian significance test, the Full Bayesian Significance Test (FBST); see Pereira, Stern andWechsler (2008) for a complete review. The disequilibrium coefficient proposed provides an easy and efficient way to make the analyses, especially if one uses Bayesian statistics. A routine in R programs (R Development Core Team, 2009) that implements the calculations is provided for the readers.

Relationship between breastfeeding duration and prevalence of posterior crossbite in the deciduous dentition

Introduction: This cross-sectional retrospective epidemiologic study assessed the relationship between exclusive breastfeeding duration and the prevalence of posterior crossbite in the deciduous dentition. Methods: Clinical examinations were performed in 1377 Brazilian children (690 boys, 687 girls), 3 to 6 years old, from 11 public schools in Sao Paulo, Brazil. Based on questionnaires answered by the parents, the children were classified into 4 groups according to the duration of exclusive breastfeeding: G1, never (119 subjects); G2, less than 6 months (720 subjects); G3, 6 to 12 months (312 subjects); and G4, more than 12 months (226 subjects). The statistical analyses included the chi-square test (P < 0.05) and the odds ratio. Results: The posterior crossbite was observed in 31.1%, 22.4%, 8.3%, and 2.2% of the children, in groups G1, G2, G3, and G4, respectively. The results showed a statistically significant relationship between exclusive breastfeeding duration and the prevalence of posterior crossbite. Conclusions: Children who were breastfed for more than 12 months had a 20-fold lower risk for the development of posterior crossbite compared with children who were never breastfed and a 5-fold lower risk compared with those breastfed between 6 and 12 months. (Am J Orthod Dentofacial Orthop 2010; 137:54-8)

Bayesian nonlinear regression models with scale mixtures of skew-normal distributions: Estimation and case influence diagnostics

The purpose of this paper is to develop a Bayesian analysis for nonlinear regression models under scale mixtures of skew-normal distributions. This novel class of models provides a useful generalization of the symmetrical nonlinear regression models since the error distributions cover both skewness and heavy-tailed distributions such as the skew-t, skew-slash and the skew-contaminated normal distributions. The main advantage of these class of distributions is that they have a nice hierarchical representation that allows the implementation of Markov chain Monte Carlo (MCMC) methods to simulate samples from the joint posterior distribution. In order to examine the robust aspects of this flexible class, against outlying and influential observations, we present a Bayesian case deletion influence diagnostics based on the Kullback-Leibler divergence. Further, some discussions on the model selection criteria are given. The newly developed procedures are illustrated considering two simulations study, and a real data previously analyzed under normal and skew-normal nonlinear regression models. (C) 2010 Elsevier B.V. All rights reserved.

On estimation and influence diagnostics for log-Birnbaum-Saunders Student-t regression models: Full Bayesian analysis

The purpose of this paper is to develop a Bayesian approach for log-Birnbaum-Saunders Student-t regression models under right-censored survival data. Markov chain Monte Carlo (MCMC) methods are used to develop a Bayesian procedure for the considered model. In order to attenuate the influence of the outlying observations on the parameter estimates, we present in this paper Birnbaum-Saunders models in which a Student-t distribution is assumed to explain the cumulative damage. Also, some discussions on the model selection to compare the fitted models are given and case deletion influence diagnostics are developed for the joint posterior distribution based on the Kullback-Leibler divergence. The developed procedures are illustrated with a real data set. (C) 2010 Elsevier B.V. All rights reserved.

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.

Bayesian Analysis for the Generalized Lognormal Distribution Applied to Failure Time Analysis

There are several versions of the lognormal distribution in the statistical literature, one is based in the exponential transformation of generalized normal distribution (GN). This paper presents the Bayesian analysis for the generalized lognormal distribution (logGN) considering independent non-informative Jeffreys distributions for the parameters as well as the procedure for implementing the Gibbs sampler to obtain the posterior distributions of parameters. The results are used to analyze failure time models with right-censored and uncensored data. The proposed method is illustrated using actual failure time data of computers.

Factors associated with stages of change for red meat and vegetable intake by Japanese-Brazilians

Stages of change assess individual motivation for lifestyle changes, contributing to the development of more effective intervention strategies. The objective of the present study was to identify factors associated with stages of change for lower intake of red meat and higher intake of vegetables in a cross-sectional analysis of 578 Japanese-Brazilians aged 30-90 years. In adjusted logistic regression models, the odds ratios for women (OR = 1.89; 95%CI: 1.154; 3.103) and physically active individuals (OR = 1.00; 95%CI: 1.000; 1.001) were positively associated with stage of "action" for the higher intake of vegetables. Inverse associations were observed between central obesity (OR = 0.5; 95%CI: 0.351; 0.887) and highest tertile of red meat intake (OR = 0.50; 95%CI: 0.302; 0.817), as well as a positive association between age (OR = 1.04; 95%CI: 1.020; 1.070) and the stage of "action" to the lower intake of meat were verified. Motivation for Japanese-Brazilians to change their food intake was linked to lifestyle. Stage of change is an important factor in mediating food intake behavior change.

Body Mass Index, Cigarette Smoking, and Alcohol Consumption and Cancers of the Oral Cavity, Pharynx, and Larynx: Modeling Odds Ratios in Pooled Case-Control Data

study-specific results, their findings should be interpreted with caution

Assessing influence in survival data with a cure fraction and covariates

Diagnostic methods have been an important tool in regression analysis to detect anomalies, such as departures from error assumptions and the presence of outliers and influential observations with the fitted models. Assuming censored data, we considered a classical analysis and Bayesian analysis assuming no informative priors for the parameters of the model with a cure fraction. A Bayesian approach was considered by using Markov Chain Monte Carlo Methods with Metropolis-Hasting algorithms steps to obtain the posterior summaries of interest. Some influence methods, such as the local influence, total local influence of an individual, local influence on predictions and generalized leverage were derived, analyzed and discussed in survival data with a cure fraction and covariates. The relevance of the approach was illustrated with a real data set, where it is shown that, by removing the most influential observations, the decision about which model best fits the data is changed.