Testing allele homogeneity: the problem of nested hypotheses

Background: The evaluation of associations between genotypes and diseases in a case-control framework plays an important role in genetic epidemiology. This paper focuses on the evaluation of the homogeneity of both genotypic and allelic frequencies. The traditional test that is used to check allelic homogeneity is known to be valid only under Hardy-Weinberg equilibrium, a property that may not hold in practice. Results: We first describe the flaws of the traditional (chi-squared) tests for both allelic and genotypic homogeneity. Besides the known problem of the allelic procedure, we show that whenever these tests are used, an incoherence may arise: sometimes the genotypic homogeneity hypothesis is not rejected, but the allelic hypothesis is. As we argue, this is logically impossible. Some methods that were recently proposed implicitly rely on the idea that this does not happen. In an attempt to correct this incoherence, we describe an alternative frequentist approach that is appropriate even when Hardy-Weinberg equilibrium does not hold. It is then shown that the problem remains and is intrinsic of frequentist procedures. Finally, we introduce the Full Bayesian Significance Test to test both hypotheses and prove that the incoherence cannot happen with these new tests. To illustrate this, all five tests are applied to real and simulated datasets. Using the celebrated power analysis, we show that the Bayesian method is comparable to the frequentist one and has the advantage of being coherent. Conclusions: Contrary to more traditional approaches, the Full Bayesian Significance Test for association studies provides a simple, coherent and powerful tool for detecting associations.

Litterfall and leaf decomposition in forest fragments under different successional phases on the Atlantic Plateau of the state of Sao Paulo, Brazil

Litterfall and litter decomposition are vital processes in tropical forests because they regulate nutrient cycling. Nutrient cycling can be altered by forest fragmentation. The Atlantic Forest is one of the most threatened biomes in the world due to human occupation over the last 500 years. This scenario has resulted in fragments of different size, age and regeneration phase. To investigate differences in litterfall and leaf decomposition between forest successional phases, we compared six forest fragments at three different successional phases and an area of mature forest on the Atlantic Plateau of Sao Paulo, Brazil. We sampled litter monthly from November 2008 to October 2009. We used litterbags to calculate leaf decomposition rate of an exotic species, Tipuana tipu (Fabaceae), over the same period litter sampling was performed. Litterfall was higher in the earliest successional area. This pattern may be related to the structural properties of the forest fragments, especially the higher abundance of pioneer species, which have higher productivity and are typical of early successional areas. However, we have not found significant differences in the decomposition rates between the studied areas, which may be caused by rapid stabilization of the decomposition environment (combined effect of microclimatic conditions and the decomposers activities). This result indicates that the leaf decomposition process have already been restored to levels observed in mature forests after a few decades of regeneration, although litterfall has not been entirely restored. This study emphasizes the importance of secondary forests for restoration of ecosystem processes on a regional scale.

Parameter recovery for a skew-normal IRT model under a Bayesian approach: hierarchical framework, prior and kernel sensitivity and sample size

Item response theory (IRT) comprises a set of statistical models which are useful in many fields, especially when there is an interest in studying latent variables (or latent traits). Usually such latent traits are assumed to be random variables and a convenient distribution is assigned to them. A very common choice for such a distribution has been the standard normal. Recently, Azevedo et al. [Bayesian inference for a skew-normal IRT model under the centred parameterization, Comput. Stat. Data Anal. 55 (2011), pp. 353-365] proposed a skew-normal distribution under the centred parameterization (SNCP) as had been studied in [R. B. Arellano-Valle and A. Azzalini, The centred parametrization for the multivariate skew-normal distribution, J. Multivariate Anal. 99(7) (2008), pp. 1362-1382], to model the latent trait distribution. This approach allows one to represent any asymmetric behaviour concerning the latent trait distribution. Also, they developed a Metropolis-Hastings within the Gibbs sampling (MHWGS) algorithm based on the density of the SNCP. They showed that the algorithm recovers all parameters properly. Their results indicated that, in the presence of asymmetry, the proposed model and the estimation algorithm perform better than the usual model and estimation methods. Our main goal in this paper is to propose another type of MHWGS algorithm based on a stochastic representation (hierarchical structure) of the SNCP studied in [N. Henze, A probabilistic representation of the skew-normal distribution, Scand. J. Statist. 13 (1986), pp. 271-275]. Our algorithm has only one Metropolis-Hastings step, in opposition to the algorithm developed by Azevedo et al., which has two such steps. This not only makes the implementation easier but also reduces the number of proposal densities to be used, which can be a problem in the implementation of MHWGS algorithms, as can be seen in [R.J. Patz and B.W. Junker, A straightforward approach to Markov Chain Monte Carlo methods for item response models, J. Educ. Behav. Stat. 24(2) (1999), pp. 146-178; R. J. Patz and B. W. Junker, The applications and extensions of MCMC in IRT: Multiple item types, missing data, and rated responses, J. Educ. Behav. Stat. 24(4) (1999), pp. 342-366; A. Gelman, G.O. Roberts, and W.R. Gilks, Efficient Metropolis jumping rules, Bayesian Stat. 5 (1996), pp. 599-607]. Moreover, we consider a modified beta prior (which generalizes the one considered in [3]) and a Jeffreys prior for the asymmetry parameter. Furthermore, we study the sensitivity of such priors as well as the use of different kernel densities for this parameter. Finally, we assess the impact of the number of examinees, number of items and the asymmetry level on the parameter recovery. Results of the simulation study indicated that our approach performed equally as well as that in [3], in terms of parameter recovery, mainly using the Jeffreys prior. Also, they indicated that the asymmetry level has the highest impact on parameter recovery, even though it is relatively small. A real data analysis is considered jointly with the development of model fitting assessment tools. The results are compared with the ones obtained by Azevedo et al. The results indicate that using the hierarchical approach allows us to implement MCMC algorithms more easily, it facilitates diagnosis of the convergence and also it can be very useful to fit more complex skew IRT models.

Ecological aspects of the sandfly fauna (Diptera, Psychodidae) in an American cutaneous leishmaniasis endemic area under the influence of hydroelectric plants in Paranapanema river, State of Parana, Brazil

Introduction: An epidemiological study was undertaken to identify determinant factors in the occurrence of American cutaneous leishmaniasis in areas under the influence of hydroelectric plants in Paranapanema river, State of Parana, Brazil. The ecological aspects of the phlebotomine fauna were investigated. Methods: Sandflies were sampled with automatic light traps from February 2004 to June 2006 at 25 sites in the urban and rural areas of Itambaraca, and in Porto Almeida and Sao Joaquim do Pontal. Results: A total of 3,187 sandflies of 15 species were captured. Nyssomyia neivai predominated (34.4%), followed by Pintomyia pessoai (32.6%), Migonemyia migonei (11.6%), Nyssomyia whitmani (8.8%), and Pintomyia fischeri (2.7%), all implicated in the transmission of Leishmania. Males predominated for Ny. neivai, and females for the other vector species, with significant statistical differences (p < 0.001). Nyssomyia neivai, Pi. pessoai, Ny. whitmani, Brumptomyia brumpti, Mg. migonei, and Pi. fischeri presented the highest values for the Standardized Species Abundance Index (SSAI). The highest frequencies and diversities were found in the preserved forest in Porto Almeida, followed by forests with degradation in Sao Joaquim do Pontal and Vila Rural. Conclusions: Sandflies were captured in all localities, with the five vectors predominating. Ny. neivai had its highest frequencies in nearby peridomestic environments and Pi. pessoai in areas of preserved forests. The highest SSAI values of Ny. neivai and Pi. pessoai reflect their wider dispersion and higher frequencies compared with other species, which seems to indicate that these two species may be transmitting leishmaniasis in the area.

An Extreme Case of Plant-Insect Codiversification: Figs and Fig-Pollinating Wasps

It is thought that speciation in phytophagous insects is often due to colonization of novel host plants, because radiations of plant and insect lineages are typically asynchronous. Recent phylogenetic comparisons have supported this model of diversification for both insect herbivores and specialized pollinators. An exceptional case where contemporaneous plant-insect diversification might be expected is the obligate mutualism between fig trees (Ficus species, Moraceae) and their pollinating wasps (Agaonidae, Hymenoptera). The ubiquity and ecological significance of this mutualism in tropical and subtropical ecosystems has long intrigued biologists, but the systematic challenge posed by >750 interacting species pairs has hindered progress toward understanding its evolutionary history. In particular, taxon sampling and analytical tools have been insufficient for large-scale cophylogenetic analyses. Here, we sampled nearly 200 interacting pairs of fig and wasp species from across the globe. Two supermatrices were assembled: on an average, wasps had sequences from 77% of 6 genes (5.6 kb), figs had sequences from 60% of 5 genes (5.5 kb), and overall 850 new DNA sequences were generated for this study. We also developed a new analytical tool, Jane 2, for event-based phylogenetic reconciliation analysis of very large data sets. Separate Bayesian phylogenetic analyses for figs and fig wasps under relaxed molecular clock assumptions indicate Cretaceous diversification of crown groups and contemporaneous divergence for nearly half of all fig and pollinator lineages. Event-based cophylogenetic analyses further support the codiversification hypothesis. Biogeographic analyses indicate that the present-day distribution of fig and pollinator lineages is consistent with a Eurasian origin and subsequent dispersal, rather than with Gondwanan vicariance. Overall, our findings indicate that the fig-pollinator mutualism represents an extreme case among plant-insect interactions of coordinated dispersal and long-term codiversification.

Comparing assemblages of Asteraceae and their insect herbivores under different land-use regimens

Accelerated tropical landscape changes occurring over recent decades have produced environmental mosaics comprising remaining isolated green areas and mixed land-use types. Our objective was to study the effects of alterations in the natural landscape on the species composition and structure of assemblages of Asteraceae and their endophagous insects through comparisons between cerrado (savanna), pastures and Eucalyptus stands. We first investigated whether similarities between assemblages of Asteraceae and their insects varied among land uses or localities. Secondly, we asked whether assemblages of Eucalyptus stands and pastures are subsets of those within the cerrado. We sampled within randomly deployed transects in 15 areas. Land use was found to be an important factor in determining plant composition similarity; however, locality did not exert any significant influence. Pastures were less similar to one another, suggesting high beta diversity. Similarities among insect assemblages were correlated with plant assemblage composition, but not with land use or locality. Species of Tephritidae were distributed along localities independently of land use. High beta diversity in Asteraceae assemblages among cerrados and pastures was supported by nestedness analysis. Plant assemblages in Eucalyptus stands were subsets of cerrado, but pasture assemblages were only partial subsets. A higher degree of nestedness in insect assemblages than in plant assemblages indicated lower beta diversity within these herbivores. Our data indicate that many herbivores are specialized on widely distributed plant genera. Conservation of Asteraceae species and their flower head insects depends not only on maintenance of landscape fragments but also on the correct matching of management form and land use. Such management may contribute to reducing isolation of plant and insect species by increasing the connectivity of remaining cerrado tracts, allowing population maintenance even at low abundances.

A comparison of deterministic, reliability-based and risk-based structural optimization under uncertainty

In this paper, the effects of uncertainty and expected costs of failure on optimum structural design are investigated, by comparing three distinct formulations of structural optimization problems. Deterministic Design Optimization (DDO) allows one the find the shape or configuration of a structure that is optimum in terms of mechanics, but the formulation grossly neglects parameter uncertainty and its effects on structural safety. Reliability-based Design Optimization (RBDO) has emerged as an alternative to properly model the safety-under-uncertainty part of the problem. With RBDO, one can ensure that a minimum (and measurable) level of safety is achieved by the optimum structure. However, results are dependent on the failure probabilities used as constraints in the analysis. Risk optimization (RO) increases the scope of the problem by addressing the compromising goals of economy and safety. This is accomplished by quantifying the monetary consequences of failure, as well as the costs associated with construction, operation and maintenance. RO yields the optimum topology and the optimum point of balance between economy and safety. Results are compared for some example problems. The broader RO solution is found first, and optimum results are used as constraints in DDO and RBDO. Results show that even when optimum safety coefficients are used as constraints in DDO, the formulation leads to configurations which respect these design constraints, reduce manufacturing costs but increase total expected costs (including expected costs of failure). When (optimum) system failure probability is used as a constraint in RBDO, this solution also reduces manufacturing costs but by increasing total expected costs. This happens when the costs associated with different failure modes are distinct. Hence, a general equivalence between the formulations cannot be established. Optimum structural design considering expected costs of failure cannot be controlled solely by safety factors nor by failure probability constraints, but will depend on actual structural configuration. (c) 2011 Elsevier Ltd. All rights reserved.

A transmission problem for beams on nonlinear supports

A transmission problem involving two Euler-Bernoulli equations modeling the vibrations of a composite beam is studied. Assuming that the beam is clamped at one extremity, and resting on an elastic bearing at the other extremity, the existence of a unique global solution and decay rates of the energy are obtained by adding just one damping device at the end containing the bearing mechanism.

Different approaches for dealing with rare variants in family-based genetic studies: application of a Genetic Analysis Workshop 17 problem

Rare variants are becoming the new candidates in the search for genetic variants that predispose individuals to a phenotype of interest. Their low prevalence in a population requires the development of dedicated detection and analytical methods. A family-based approach could greatly enhance their detection and interpretation because rare variants are nearly family specific. In this report, we test several distinct approaches for analyzing the information provided by rare and common variants and how they can be effectively used to pinpoint putative candidate genes for follow-up studies. The analyses were performed on the mini-exome data set provided by Genetic Analysis Workshop 17. Eight approaches were tested, four using the trait’s heritability estimates and four using QTDT models. These methods had their sensitivity, specificity, and positive and negative predictive values compared in light of the simulation parameters. Our results highlight important limitations of current methods to deal with rare and common variants, all methods presented a reduced specificity and, consequently, prone to false positive associations. Methods analyzing common variants information showed an enhanced sensibility when compared to rare variants methods. Furthermore, our limited knowledge of the use of biological databases for gene annotations, possibly for use as covariates in regression models, imposes a barrier to further research.

Ecological aspects of the sandfly fauna (Diptera, Psychodidae) in an American cutaneous leishmaniasis endemic area under the influence of hydroelectric plants in Paranapanema river, State of Paraná, Brazil

INTRODUCTION: An epidemiological study was undertaken to identify determinant factors in the occurrence of American cutaneous leishmaniasis in areas under the influence of hydroelectric plants in Paranapanema river, State of Paraná, Brazil. The ecological aspects of the phlebotomine fauna were investigated. METHODS: Sandflies were sampled with automatic light traps from February 2004 to June 2006 at 25 sites in the urban and rural areas of Itambaracá, and in Porto Almeida and São Joaquim do Pontal. RESULTS: A total of 3,187 sandflies of 15 species were captured. Nyssomyia neivai predominated (34.4%), followed by Pintomyia pessoai (32.6%), Migonemyia migonei (11.6%), Nyssomyia whitmani (8.8%), and Pintomyia fischeri (2.7%), all implicated in the transmission of Leishmania. Males predominated for Ny. neivai, and females for the other vector species, with significant statistical differences (p < 0.001). Nyssomyia neivai, Pi. pessoai, Ny. whitmani, Brumptomyia brumpti, Mg. migonei, and Pi. fischeri presented the highest values for the Standardized Species Abundance Index (SSAI). The highest frequencies and diversities were found in the preserved forest in Porto Almeida, followed by forests with degradation in São Joaquim do Pontal and Vila Rural. CONCLUSIONS: Sandflies were captured in all localities, with the five vectors predominating. Ny. neivai had its highest frequencies in nearby peridomestic environments and Pi. pessoai in areas of preserved forests. The highest SSAI values of Ny. neivai and Pi. pessoai reflect their wider dispersion and higher frequencies compared with other species, which seems to indicate that these two species may be transmitting leishmaniasis in the area.

Models for capacitated lot-sizing problem with backlogging, setup carryover and crossover

Setup operations are significant in some production environments. It is mandatory that their production plans consider some features, as setup state conservation across periods through setup carryover and crossover. The modelling of setup crossover allows more flexible decisions and is essential for problems with long setup times. This paper proposes two models for the capacitated lot-sizing problem with backlogging and setup carryover and crossover. The first is in line with other models from the literature, whereas the second considers a disaggregated setup variable, which tracks the starting and completion times of the setup operation. This innovative approach permits a more compact formulation. Computational results show that the proposed models have outperformed other state-of-the-art formulation.