Hypothesis testing for differentially correlated features.
Ontology highlight
ABSTRACT:
SUBMITTER: Sheng E
PROVIDER: S-EPMC5031944 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Given data obtained under two sampling conditions, it is often of interest to identify variables that behave differently in one condition than in the other. We introduce a method for differential analysis of second-order behavior called Differential Correlation Mining (DCM). The DCM method identifies differentially correlated sets of variables, with the property that the average pairwise correlation between variables in a set is higher under one sample condition than the other. DCM is based on an iterative search procedure that adaptively updates the size and elements of a candidate variable set. Updates are performed via hypothesis testing of individual variables, based on the asymptotic distribution of their average differential correlation. We investigate the performance of DCM by applying it to simulated data as well as to recent experimental datasets in genomics and brain imaging.
| S-EPMC6927674 | biostudies-literature
Project description:Male mate choice might be based on both absolute and relative strategies. Cues of female attractiveness are thus likely to reflect both fitness and reproductive potential, as well as compatibility with particular male phenotypes. In humans, absolute clues of fertility and indices of favorable developmental stability are generally associated with increased women's attractiveness. However, why men exhibit variable preferences remains less studied. Male mate choice might be influenced by uncertainty of paternity, a selective factor in species where the survival of the offspring depends on postnatal paternal care. For instance, in humans, a man might prefer a woman with recessive traits, thereby increasing the probability that his paternal traits will be visible in the child and ensuring paternity. Alternatively, attractiveness is hypothesized to be driven by self-resembling features (homogamy), which would reduce outbreeding depression. These hypotheses have been simultaneously evaluated for various facial traits using both real and artificial facial stimuli. The predicted preferences were then compared to realized mate choices using facial pictures from couples with at least 1 child. No evidence was found to support the paternity uncertainty hypothesis, as recessive features were not preferred by male raters. Conversely, preferences for self-resembling mates were found for several facial traits (hair and eye color, chin dimple, and thickness of lips and eyebrows). Moreover, realized homogamy for facial traits was also found in a sample of long-term mates. The advantages of homogamy in evolutionary terms are discussed.
| S-EPMC3504097 | biostudies-literature
Project description:BackgroundFor many practical hypothesis testing (H-T) applications, the data are correlated and/or with heterogeneous variance structure. The regression t-test for weighted linear mixed-effects regression (LMER) is a legitimate choice because it accounts for complex covariance structure; however, high computational costs and occasional convergence issues make it impractical for analyzing high-throughput data. In this paper, we propose computationally efficient parametric and semiparametric tests based on a set of specialized matrix techniques dubbed as the PB-transformation. The PB-transformation has two advantages: 1. The PB-transformed data will have a scalar variance-covariance matrix. 2. The original H-T problem will be reduced to an equivalent one-sample H-T problem. The transformed problem can then be approached by either the one-sample Student's t-test or Wilcoxon signed rank test.ResultsIn simulation studies, the proposed methods outperform commonly used alternative methods under both normal and double exponential distributions. In particular, the PB-transformed t-test produces notably better results than the weighted LMER test, especially in the high correlation case, using only a small fraction of computational cost (3 versus 933 s). We apply these two methods to a set of RNA-seq gene expression data collected in a breast cancer study. Pathway analyses show that the PB-transformed t-test reveals more biologically relevant findings in relation to breast cancer than the weighted LMER test.ConclusionsAs fast and numerically stable replacements for the weighted LMER test, the PB-transformed tests are especially suitable for "messy" high-throughput data that include both independent and matched/repeated samples. By using our method, the practitioners no longer have to choose between using partial data (applying paired tests to only the matched samples) or ignoring the correlation in the data (applying two sample tests to data with some correlated samples). Our method is implemented as an R package 'PBtest' and is available at https://github.com/yunzhang813/PBtest-R-Package .
| S-EPMC6466736 | biostudies-literature
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: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: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:A phylogenetic model of sequence evolution for a set of n taxa is a collection of probability distributions on the 4 n possible site patterns that may be observed in their aligned DNA sequences. For a four-taxon model, one can arrange the entries of these probability distributions into three flattening matrices that correspond to the three different unrooted leaf-labeled four-leaf trees, or quartet trees. The flattening matrix corresponding to the tree parameter of the model is known to satisfy certain rank conditions. Methods such as ErikSVD and SVDQuartets take advantage of this observation by applying singular value decomposition to flattening matrices consisting of empirical data. Each possible quartet is assigned an "SVD score" based on how close the flattening is to the set of matrices of the predicted rank. When choosing among possible quartets, the one with the lowest score is inferred to be the phylogeny of the four taxa under consideration. Since an n-leaf phylogenetic tree is determined by its quartets, this approach can be generalized to infer larger phylogenies. In this article, we explore using the SVD score as a test statistic to test whether phylogenetic data were generated by a particular quartet tree. To do so, we use several results to approximate the distribution of the SVD score and to give upper bounds on the p-value of the associated hypothesis tests. We also apply these hypothesis tests to simulated phylogenetic data and discuss the implications for interpreting SVD scores in rank-based inference methods.
| S-EPMC8283673 | biostudies-literature