Hypothesis Testing With Rank Conditions in Phylogenetics.
Ontology highlight
ABSTRACT:
SUBMITTER: Long C
PROVIDER: S-EPMC8283673 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Functional data arise frequently in biomedical studies, where it is often of interest to investigate the association between functional predictors and a scalar response variable. While functional linear models (FLM) are widely used to address these questions, hypothesis testing for the functional association in the FLM framework remains challenging. A popular approach to testing the functional effects is through dimension reduction by functional principal component (PC) analysis. However, its power performance depends on the choice of the number of PCs, and is not systematically studied. In this article, we first investigate the power performance of the Wald-type test with varying thresholds in selecting the number of PCs for the functional covariates, and show that the power is sensitive to the choice of thresholds. To circumvent the issue, we propose a new method of ordering and selecting principal components to construct test statistics. The proposed method takes into account both the association with the response and the variation along each eigenfunction. We establish its theoretical properties and assess the finite sample properties through simulations. Our simulation results show that the proposed test is more robust against the choice of threshold while being as powerful as, and often more powerful than, the existing method. We then apply the proposed method to the cerebral white matter tracts data obtained from a diffusion tensor imaging tractography study.
| S-EPMC5518697 | biostudies-literature
Project description:In a multivariate setting, we consider the task of identifying features whose correlations with the other features differ across conditions. Such correlation shifts may occur independently of mean shifts, or differences in the means of the individual features across conditions. Previous approaches for detecting correlation shifts consider features simultaneously, by computing a correlation-based test statistic for each feature. However, since correlations involve two features, such approaches do not lend themselves to identifying which feature is the culprit. In this article, we instead consider a serial testing approach, by comparing columns of the sample correlation matrix across two conditions, and removing one feature at a time. Our method provides a novel perspective and favorable empirical results compared with competing approaches.
| S-EPMC5031944 | biostudies-literature
Project description:Human cytosolic sulfotransferases (SULTs) regulate the activities of hundreds of signaling metabolites via transfer of the sulfuryl moiety (-SO3) from activated sulfate (3'-phosphoadenosine 5'-phosphosulfate) to the hydroxyls and primary amines of xeno- and endobiotics. How SULTs select substrates from the scores of competing ligands present in a cytosolic milieu is an important issue in the field. Selectivity appears to be sterically controlled by a molecular pore that opens and closes in response to nucleotide binding. This point of view is fostered by structures showing nucleotide-dependent pore closure and the fact that nucleotide binding induces an isomerization that restricts access to the acceptor-binding pocket. Molecular dynamics models underscore the importance of pore isomerization in selectivity and predict that specific molecular linkages stabilize the closed pore in response to nucleotide binding. To test the pore model, these linkages were disrupted in SULT2A1 via mutagenesis, and the effects on selectivity were determined. The mutations uncoupled nucleotide binding from selectivity and produced enzymes that no longer discriminated between large and small substrates. The mutations did not affect the affinity or turnover of small substrates but resulted in a 183-fold gain in catalytic efficiently toward large substrates. Models predict that an 11-residue "flap" covering the acceptor-binding pocket can open and admit large substrates when nucleotide is bound; a mutant structure demonstrated that this is so. In summary, the model was shown to be a robust, accurate predictor of SULT structure and selectivity whose general features will likely apply to other members of the SULT family.
| S-EPMC3605680 | biostudies-literature
Project description:We consider large-scale studies in which thousands of significance tests are performed simultaneously. In some of these studies, the multiple testing procedure can be severely biased by latent confounding factors such as batch effects and unmeasured covariates that correlate with both primary variable(s) of interest (e.g., treatment variable, phenotype) and the outcome. Over the past decade, many statistical methods have been proposed to adjust for the confounders in hypothesis testing. We unify these methods in the same framework, generalize them to include multiple primary variables and multiple nuisance variables, and analyze their statistical properties. In particular, we provide theoretical guarantees for RUV-4 [Gagnon-Bartsch, Jacob and Speed (2013)] and LEAPP [Ann. Appl. Stat. 6 (2012) 1664-1688], which correspond to two different identification conditions in the framework: the first requires a set of "negative controls" that are known a priori to follow the null distribution; the second requires the true nonnulls to be sparse. Two different estimators which are based on RUV-4 and LEAPP are then applied to these two scenarios. We show that if the confounding factors are strong, the resulting estimators can be asymptotically as powerful as the oracle estimator which observes the latent confounding factors. For hypothesis testing, we show the asymptotic z-tests based on the estimators can control the type I error. Numerical experiments show that the false discovery rate is also controlled by the Benjamini-Hochberg procedure when the sample size is reasonably large.
| S-EPMC6706069 | biostudies-literature
Project description:The ratites have stimulated much debate as to how such large flightless birds came to be distributed across the southern continents, and whether they are a monophyletic group or are composed of unrelated lineages that independently lost the power of flight. Hypotheses regarding the relationships among taxa differ for morphological and molecular data sets, thus hindering attempts to test whether plate tectonic events can explain ratite biogeography. Here, we present the complete mitochondrial DNA genomes of two extinct moas from New Zealand, along with those of five extant ratites (the lesser rhea, the ostrich, the great spotted kiwi, the emu and the southern cassowary and two tinamous from different genera. The non-stationary base composition in these sequences violates the assumptions of most tree-building methods. When this bias is corrected using neighbour-joining with log-determinant distances and non-homogeneous maximum likelihood, the ratites are found to be monophlyletic, with moas basal, as in morphological trees. The avian sequences also violate a molecular clock, so we applied a non-parametric rate smoothing algorithm, which minimizes ancestor-descendant local rate changes, to date nodes in the tree. Using this method, most of the major ratite lineages fit the vicariance biogeography hypothesis, the exceptions being the ostrich and the kiwi, which require dispersal to explain their present distribution.
| S-EPMC1088691 | biostudies-other
Project description:Summary Causal Bayesian Networks are a special class of Bayesian networks in which the hierarchy directly encodes the causal relationships between the variables. This allows to compute the effect of interventions, which are external changes to the system, caused by e.g. gene knockouts or an administered drug. Whereas numerous packages for constructing causal Bayesian networks are available, hardly any program targeted at downstream analysis exists. In this paper we present CausalTrail, a tool for performing reasoning on causal Bayesian networks using the do-calculus. CausalTrail's features include multiple data import methods, a flexible query language for formulating hypotheses, as well as an intuitive graphical user interface. The program is able to account for missing data and thus can be readily applied in multi-omics settings where it is common that not all measurements are performed for all samples. Availability and Implementation CausalTrail is implemented in C++ using the Boost and Qt5 libraries. It can be obtained from https://github.com/dstoeckel/causaltrail.
| S-EPMC4743151 | biostudies-other
Project description:A number of hypotheses about compensatory mechanisms that allow ectothermic animals to cope with the latitudinal decrease in ambient temperature ( TA) have been proposed during the last century. One of these hypotheses, the 'metabolic homeostasis' hypothesis (MHH), states that species should show the highest thermal sensitivity of the metabolic rate ( Q10-SMR) at the colder end of the range of TAs they usually experience in nature. This way, species should be able to minimize maintenance costs during the colder hours of the day, but quickly take advantage of increases in TA during the warmer parts of the day. Here, we created a dataset that includes Q10-SMR values for 58 amphibian species, assessed at four thermal ranges, to evaluate three predictions derived from the MHH. In line with this hypothesis, we found that: (i) Q10-SMR values tended to be positively correlated with latitude when measured at lower TAs, but negative correlated with latitude when measured at higher TAs, (ii) Q10-SMR measured at lower TAs were higher in temperate species, whereas Q10-SMR measured at higher TAs were higher in tropical species, and (iii) the experimental TA at which Q10-SMR was maximal for each species decreased with latitude. This is the first study to our knowledge showing that the relationship between Q10-SMR and latitude in ectotherms changes with the TA at which Q10-SMR is assessed, as predicted from an adaptive hypothesis. This article is part of the theme issue 'Physiological diversity, biodiversity patterns and global climate change: testing key hypotheses involving temperature and oxygen'.
| S-EPMC6606461 | biostudies-literature
Project description:Although evolution by natural selection is widely regarded as the most important principle of biology, it is unknown whether phenotypic variations within and between species are mostly adaptive or neutral due to the lack of relevant studies of large, unbiased samples of phenotypic traits. Here, we examine 210 yeast morphological traits chosen because of experimental feasibility irrespective of their potential adaptive values. Our analysis is based on the premise that, under neutrality, the rate of phenotypic evolution measured in the unit of mutational size declines as the trait becomes more important to fitness, analogous to the neutral paradigm that functional genes evolve more slowly than functionless pseudogenes. However, we find faster evolution of more important morphological traits within and between species, rejecting the neutral hypothesis. By contrast, an analysis of 3,466 gene expression traits fails to refute neutrality. Thus, at least in yeast, morphological evolution appears largely adaptive, but the same may not apply to other classes of phenotypes. Our neutrality test is applicable to other species, especially genetic model organisms, for which estimations of mutational size and trait importance are relatively straightforward.
| S-EPMC5699056 | biostudies-literature
Project description:BackgroundSocial capital research has consistently shown positive associations between generalised trust and health outcomes over 2 decades. Longitudinal studies attempting to test causal relationships further support the theory that trust is an independent predictor of health. However, as the reverse causality hypothesis has yet to be empirically tested, a knowledge gap remains. The aim of this study, therefore, was to investigate if health status predicts trust.MethodsData employed in this study came from 4 waves of the British Household Panel Survey between years 2000 and 2007 (N=8114). The sample was stratified by baseline trust to investigate temporal relationships between prior self-rated health (SRH) and changes in trust. We used logistic regression models with random effects, as trust was expected to be more similar within the same individuals over time.ResultsFrom the 'Can trust at baseline' cohort, poor SRH at time (t-1) predicted low trust at time (t) (OR=1.38). Likewise, good health predicted high trust within the 'Cannot' trust cohort (OR=1.30). These patterns of positive association remained after robustness checks, which adjusted for misclassification of outcome (trust) status and the existence of other temporal pathways.ConclusionsThis study offers empirical evidence to support the circular nature of trust/health relationship. The stability of association between prior health status and changes in trust over time differed between cohorts, hinting at the existence of complex pathways rather than a simple positive feedback loop.
| S-EPMC4717376 | biostudies-literature
Project description:In this paper, we study the detection boundary for minimax hypothesis testing in the context of high-dimensional, sparse binary regression models. Motivated by genetic sequencing association studies for rare variant effects, we investigate the complexity of the hypothesis testing problem when the design matrix is sparse. We observe a new phenomenon in the behavior of detection boundary which does not occur in the case of Gaussian linear regression. We derive the detection boundary as a function of two components: a design matrix sparsity index and signal strength, each of which is a function of the sparsity of the alternative. For any alternative, if the design matrix sparsity index is too high, any test is asymptotically powerless irrespective of the magnitude of signal strength. For binary design matrices with the sparsity index that is not too high, our results are parallel to those in the Gaussian case. In this context, we derive detection boundaries for both dense and sparse regimes. For the dense regime, we show that the generalized likelihood ratio is rate optimal; for the sparse regime, we propose an extended Higher Criticism Test and show it is rate optimal and sharp. We illustrate the finite sample properties of the theoretical results using simulation studies.
| S-EPMC4522432 | biostudies-literature