A Comparative Study of Statistical Methods Used to Identify

Total Page:16

File Type:pdf, Size:1020Kb

A Comparative Study of Statistical Methods Used to Identify Briefings in Bioinformatics Advance Access published August 20, 2013 BRIEFINGS IN BIOINFORMATICS. page1of13 doi:10.1093/bib/bbt051 A comparative study of statistical methods used to identify dependencies between gene expression signals Suzana de Siqueira Santos, Daniel YasumasaTakahashi, Asuka Nakata and Andre¤Fujita Submitted: 18th April 2013; Received (in revised form): 10th June 2013 Abstract One major task in molecular biology is to understand the dependency among genes to model gene regulatory net- works. Pearson’s correlation is the most common method used to measure dependence between gene expression signals, but it works well only when data are linearly associated. For other types of association, such as non-linear Downloaded from or non-functional relationships, methods based on the concepts of rank correlation and information theory-based measures are more adequate than the Pearson’s correlation, but are less used in applications, most probably because of a lack of clear guidelines for their use. This work seeks to summarize the main methods (Pearson’s, Spearman’s and Kendall’s correlations; distance correlation; Hoeffding’s D measure; Heller^Heller^Gorfine measure; mutual in- http://bib.oxfordjournals.org/ formation and maximal information coefficient) used to identify dependency between random variables, especially gene expression data, and also to evaluate the strengths and limitations of each method. Systematic Monte Carlo simulation analyses ranging from sample size, local dependence and linear/non-linear and also non-functional relationships are shown. Moreover, comparisons in actual gene expression data are carried out. Finally, we provide a suggestive list of methods that can be used for each type of data set. Keywords: correlation; dependence; network; data analysis at Princeton University on August 25, 2013 INTRODUCTION particular measure of stochastic dependence among In bioinformatics, the notion of dependence is cen- many. It is the canonical measure in the world of tral to model gene regulatory networks. The func- multivariate normal distributions, and more generally tional relationships and interactions among genes for spherical and elliptical distributions. However, and their products are usually inferred by using stat- empirical research in molecular biology shows that istical methods that identify dependence among the distributions of gene expression signals may signals [1–4]. not belong in this class [4]. To identify associations One well-known method to identify dependence not limited to linear associations, but dependent in between gene expression signals is the Pearson’s cor- a broad sense, several methods have been de- relation. But Pearson’s correlation, although it is one veloped, most of them based on ranks or information of the most ubiquitous concepts in modern molecu- theory. lar biology, is also one of the most misunderstood The main aim of this article is to clarify the essen- concepts. Some of the confusion may arise from the tial mathematical ideas behind several methods literary use of the word to cover any notion of de- [Pearson’s correlation coefficient [5,6], Spearman’s pendence. To a statistician, correlation is only one correlation coefficient [7], Kendall’s correlation Corresponding author. Andre´ Fujita, Department of Computer Science, Institute of Mathematics and Statistics, University of Sa˜o Paulo, Rua do Mata˜o, 1010–Cidade Universita´ria–SP, 05508-090, Brazil. Tel.: þ 55 11 3091 5177; Fax.: þ 55 11 3091 6134; E-mail: [email protected] Suzana de Siqueira Santos obtained her BS in Computer Science in 2012 at the University of Sa˜o Paulo. She is currently a master course student in Computer Science at the same university. DanielYasumasaTakahashi received his Ph.D. in Bioinformatics in 2009 at the University of Sa˜o Paulo. He is currently a post-doc at Princeton University. Asuka Nakata received her Ph.D. in Biosciences in 2009. She is currently a post-doc at the University of Tokyo. Andre¤Fujita received his Ph.D. in Bioinformatics in 2007. He is currently Assistant Professor at the University of Sa˜o Paulo. ß The Author 2013. Published by Oxford University Press. For Permissions, please email: [email protected] page 2 of 13 de Siqueira Santos et al. coefficient [8], distance correlation [9], Hoeffding’s MEASURES OF DEPENDENCE D measure [10], Heller–Heller–Gorfine (HHG) BETWEEN RANDOM VARIABLES measure [11], mutual information (MI) [12] and Let ðÞx1, y1 , ðÞx2, y2 , ..., ðÞxn, yn be a set of joint n maximal information coefficient (MIC) [13]] used observations from two univariate random variables X to identify dependence between random variables and Y. that anyone wishing to model dependent phenom- The test of dependence between X and Y is ena should know. Furthermore, we illustrate by described as a hypothesis test as follows: Monte Carlo simulations where we varied sample size, local dependence and linear/non-linear and H0: X and Y are not dependent (null hypothesis). also non-functional relationships, the strengths and H1: X and Yare dependent (alternative hypothesis). limitations concerning the different measures used to identify dependent signals. Finally, we illustrate In the next section, we will describe frequently used the application of the methods in actual gene ex- methods to identify dependent data and also show by pression data with known dependence structure be- simulations that some methods such as Pearson’s, tween the genes. Thus we hope to provide the Spearman’s and Kendall’s correlations can only Downloaded from necessary elements for a better comprehension of detect linear or non-linear monotonic (strictly the methods and also the choice of a suitable increasing or strictly decreasing function, i.e. a func- dependence test method, based on practical con- tion that preserves the given order) relationships, straints and goals. whereas others such as distance correlation, HHG measure, Hoeffding’s D measure and MI are able http://bib.oxfordjournals.org/ to identify non-monotonic and non-functional rela- tionships also. Furthermore, we will see that STATISTICAL INDEPENDENCE although MIC is not mathematically proven to be BETWEEN TWO RANDOM consistent against all general alternatives, it can detect VA R I A BLES some non-monotonic relationships. Statistical independence indicates that there is no re- lation between two random variables. If the variables at Princeton University on August 25, 2013 are statistically independent, then the distribution of Pearson’s product-moment correlation one of them is the same no matter at which fixed coefficient levels the other variable is considered, and observa- The Pearson’s product-moment correlation or tions for such variables will lead correspondingly to simply Pearson’s correlation [5,6] is a measure of nearly equal frequency distributions. On the other linear dependence, as the slope obtained by the hand, if there is dependence, then the levels of one linear regression of Y by X is Pearson’s correlation of the variables vary with changing levels of the multiplied byP that ratio ofP standard deviations. n n x y other. In other words, under independence, know- i¼1 i i¼1 i Let x ¼ n and y ¼ n be the means of X ledge about one feature remains unaffected by infor- and Y, respectively. Then, the Pearson’s correlation mation provided about the other, whereas under coefficient rPearson is defined as follows: dependence, it follows which level of one variable P n ðÞÀ ðÞÀ occurs as soon as the level of the other variable is qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffii¼1 xi x yi y rPearsonðÞ¼X, Y P P n 2 n 2 known. i¼1 ðÞxi À x i¼1 ðÞyi À y Formally, two random variables X and Y with cumulative distribution functions FXðÞx and FY ðÞy , For joint normal distributions, Pearson’s correl- and probability densities fX ðÞx and fY ðÞy , are inde- ation coefficient under H0 follows a Student’s pendent if and only if the combined random variable t-distribution with n À 2 degrees of freedom. The ðÞX, Y has a joint cumulative distribution function t statistic is as follows: pffiffiffiffiffiffiffiffiffiffiffi FX, Y ðÞ¼x, y FXðÞx FY ðÞy , or, equivalently, a joint r ðÞX, Y n À 2 density f ðÞ¼x, y f ðÞx f ðÞy . We say that two t ¼ pPearsonffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi X, Y X Y 1 À r2 ðÞX, Y random variables X and Y are dependent if they Pearson are not independent. The problem then is how to When the random variables are not jointly nor- measure and detect dependence from the observa- mally distributed, the Fisher’s transformation is used tion of the two random variables. to get an asymptotic normal distribution. Correlation and dependence in biological data page 3 of 13 In the case of perfect linear dependence, we have monotonically increasing relationship (for all x1 and rPearsonðÞ¼ÆX, Y 1. The Pearson correlation is þ1 x2 such that x1 < x2,wehavey1 < y2), and À1in in the case of a perfect positive (increasing) linear the case of a perfect monotonically decreasing rela- relationship and À1 in the case of a perfect negative tionship (for all x1 and x2 such that x1 < x2,wehave (decreasing) linear relationship. In the case of linearly y1 > y2). In the case of monotonically independent independent random variables, rPearsonðÞ¼X, Y 0, random variables, rSpearmanðÞ¼X, Y 0, and in the and in the case of imperfect linear dependence, case of imperfect monotonically dependence, À1 < rPearsonðÞX, Y < 1. These last two cases are À1 < rSpearmanðÞX, Y < 1. Again, similarly to the ones for which misinterpretations of correlation Pearson’s correlation coefficient, rSpearmanðÞ¼X, Y 0 are possible because it is usually assumed that non- does not mean that random variables X and Y are correlated X and Y means independent variables, independent, but only that they are monotonically whereas in fact, they may be associated in a non- independent. linear fashion that Pearson’s correlation coefficient The R function for Spearman’s test is cor.test with is not able to identify.
Recommended publications
  • Testing the Causal Direction of Mediation Effects in Randomized Intervention Studies
    Prevention Science (2019) 20:419–430 https://doi.org/10.1007/s11121-018-0900-y Testing the Causal Direction of Mediation Effects in Randomized Intervention Studies Wolfgang Wiedermann1 & Xintong Li1 & Alexander von Eye2 Published online: 21 May 2018 # Society for Prevention Research 2018 Abstract In a recent update of the standards for evidence in research on prevention interventions, the Society of Prevention Research emphasizes the importance of evaluating and testing the causal mechanism through which an intervention is expected to have an effect on an outcome. Mediation analysis is commonly applied to study such causal processes. However, these analytic tools are limited in their potential to fully understand the role of theorized mediators. For example, in a design where the treatment x is randomized and the mediator (m) and the outcome (y) are measured cross-sectionally, the causal direction of the hypothesized mediator-outcome relation is not uniquely identified. That is, both mediation models, x → m → y or x → y → m, may be plausible candidates to describe the underlying intervention theory. As a third explanation, unobserved confounders can still be responsible for the mediator-outcome association. The present study introduces principles of direction dependence which can be used to empirically evaluate these competing explanatory theories. We show that, under certain conditions, third higher moments of variables (i.e., skewness and co-skewness) can be used to uniquely identify the direction of a mediator-outcome relation. Significance procedures compatible with direction dependence are introduced and results of a simulation study are reported that demonstrate the performance of the tests. An empirical example is given for illustrative purposes and a software implementation of the proposed method is provided in SPSS.
    [Show full text]
  • Local White Matter Architecture Defines Functional Brain Dynamics
    2018 IEEE International Conference on Systems, Man, and Cybernetics Local White Matter Architecture Defines Functional Brain Dynamics Yo Joong Choe* Sivaraman Balakrishnan Aarti Singh Kakao Dept. of Statistics and Data Science Machine Learning Dept. Seongnam-si, South Korea Carnegie Mellon University Carnegie Mellon University [email protected] Pittsburgh, PA, USA Pittsburgh, PA, USA [email protected] [email protected] Jean M. Vettel Timothy Verstynen U.S. Army Research Laboratory Dept. of Psychology and CNBC Aberdeen Proving Ground, MD, USA Carnegie Mellon University [email protected] Pittsburgh, PA, USA [email protected] Abstract—Large bundles of myelinated axons, called white the integrity of the myelin sheath is critical for synchroniz- matter, anatomically connect disparate brain regions together ing information transmission between distal brain areas [2], and compose the structural core of the human connectome. We fostering the ability of these networks to adapt over time recently proposed a method of measuring the local integrity along the length of each white matter fascicle, termed the local [4]. Thus, variability in the myelin sheath, as well as other connectome [1]. If communication efficiency is fundamentally cellular support mechanisms, would contribute to variability constrained by the integrity along the entire length of a white in functional coherence across the circuit. matter bundle [2], then variability in the functional dynamics To study the integrity of structural connectivity, we recently of brain networks should be associated with variability in the introduced the concept of the local connectome. This is local connectome. We test this prediction using two statistical ap- proaches that are capable of handling the high dimensionality of defined as the pattern of fiber systems (i.e., number of fibers, data.
    [Show full text]
  • From Distance Correlation to Multiscale Graph Correlation Arxiv
    From Distance Correlation to Multiscale Graph Correlation Cencheng Shen∗1, Carey E. Priebey2, and Joshua T. Vogelsteinz3 1Department of Applied Economics and Statistics, University of Delaware 2Department of Applied Mathematics and Statistics, Johns Hopkins University 3Department of Biomedical Engineering and Institute of Computational Medicine, Johns Hopkins University October 2, 2018 Abstract Understanding and developing a correlation measure that can detect general de- pendencies is not only imperative to statistics and machine learning, but also crucial to general scientific discovery in the big data age. In this paper, we establish a new framework that generalizes distance correlation — a correlation measure that was recently proposed and shown to be universally consistent for dependence testing against all joint distributions of finite moments — to the Multiscale Graph Correla- tion (MGC). By utilizing the characteristic functions and incorporating the nearest neighbor machinery, we formalize the population version of local distance correla- tions, define the optimal scale in a given dependency, and name the optimal local correlation as MGC. The new theoretical framework motivates a theoretically sound ∗[email protected] arXiv:1710.09768v3 [stat.ML] 30 Sep 2018 [email protected] [email protected] 1 Sample MGC and allows a number of desirable properties to be proved, includ- ing the universal consistency, convergence and almost unbiasedness of the sample version. The advantages of MGC are illustrated via a comprehensive set of simula- tions with linear, nonlinear, univariate, multivariate, and noisy dependencies, where it loses almost no power in monotone dependencies while achieving better perfor- mance in general dependencies, compared to distance correlation and other popular methods.
    [Show full text]
  • From Distance Correlation to Multiscale Graph Correlation
    From Distance Correlation to Multiscale Graph Correlation Cencheng Shen∗1, Carey E. Priebey2, and Joshua T. Vogelsteinz3 1Department of Applied Economics and Statistics, University of Delaware 2Department of Applied Mathematics and Statistics, Johns Hopkins University 3Department of Biomedical Engineering and Institute of Computational Medicine, Johns Hopkins University January 19, 2019 Abstract Understanding and developing a correlation measure that can detect general de- pendencies is not only imperative to statistics and machine learning, but also crucial to general scientific discovery in the big data age. In this paper, we establish a new framework that generalizes distance correlation — a correlation measure that was recently proposed and shown to be universally consistent for dependence testing against all joint distributions of finite moments — to the Multiscale Graph Correla- tion (MGC). By utilizing the characteristic functions and incorporating the nearest neighbor machinery, we formalize the population version of local distance correla- tions, define the optimal scale in a given dependency, and name the optimal local correlation as MGC. The new theoretical framework motivates a theoretically sound ∗[email protected] [email protected] [email protected] 1 Sample MGC and allows a number of desirable properties to be proved, includ- ing the universal consistency, convergence and almost unbiasedness of the sample version. The advantages of MGC are illustrated via a comprehensive set of simula- tions with linear, nonlinear, univariate, multivariate, and noisy dependencies, where it loses almost no power in monotone dependencies while achieving better perfor- mance in general dependencies, compared to distance correlation and other popular methods. Keywords: testing independence, generalized distance correlation, nearest neighbor graph 2 1 Introduction p q Given pairs of observations (xi; yi) 2 R ×R for i = 1; : : : ; n, assume they are generated by independently identically distributed (iid) FXY .
    [Show full text]
  • University of Aberdeen the Distance Standard Deviation Edelmann, Dominic; Richards, Donald ; Vogel, Daniel
    University of Aberdeen The Distance Standard Deviation Edelmann, Dominic; Richards, Donald ; Vogel, Daniel Publication date: 2017 Document Version Other version Link to publication Citation for published version (APA): Edelmann, D., Richards, D., & Vogel, D. (2017). The Distance Standard Deviation. (pp. 1-27). ArXiv. https://arxiv.org/abs/1705.05777 General rights Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. ? Users may download and print one copy of any publication from the public portal for the purpose of private study or research. ? You may not further distribute the material or use it for any profit-making activity or commercial gain ? You may freely distribute the URL identifying the publication in the public portal ? Take down policy If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim. Download date: 02. Oct. 2021 The Distance Standard Deviation Dominic Edelmann,∗ Donald Richards,y and Daniel Vogelz May 17, 2017 Abstract The distance standard deviation, which arises in distance correlation analy- sis of multivariate data, is studied as a measure of spread. New representations for the distance standard deviation are obtained in terms of Gini's mean differ- ence and in terms of the moments of spacings of order statistics. Inequalities for the distance variance are derived, proving that the distance standard deviation is bounded above by the classical standard deviation and by Gini's mean difference.
    [Show full text]
  • Comparative Evaluation of Statistical Dependence Measures
    University of Arkansas, Fayetteville ScholarWorks@UARK Theses and Dissertations 12-2020 Comparative Evaluation of Statistical Dependence Measures Eman Abdel Rahman Ibrahim University of Arkansas, Fayetteville Follow this and additional works at: https://scholarworks.uark.edu/etd Part of the Applied Statistics Commons Citation Ibrahim, E. A. (2020). Comparative Evaluation of Statistical Dependence Measures. Theses and Dissertations Retrieved from https://scholarworks.uark.edu/etd/3903 This Thesis is brought to you for free and open access by ScholarWorks@UARK. It has been accepted for inclusion in Theses and Dissertations by an authorized administrator of ScholarWorks@UARK. For more information, please contact [email protected]. Comparative Evaluation of Statistical Dependence Measures A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Statistics and Analytics by Eman Abdel Rahman Ibrahim University of Arkansas Bachelor of Arts in Mathematics, 2015 December 2020 University of Arkansas This thesis is approved for recommendation to the Graduate Council. ______________________________ Qingyang Zhang, PhD. Thesis Director _______________________________ _____________________________ Mark E. Arnold, PhD. Giovanni Petris, PhD. Committee Member Committee Member Abstract Measuring and testing dependence between random variables is of great importance in many scientific fields. In the case of linearly correlated variables, Pearson’s correlation coefficient is a commonly used measure of the correlation
    [Show full text]
  • Measuring and Testing Dependence by Correlation of Distances
    The Annals of Statistics 2007, Vol. 35, No. 6, 2769–2794 DOI: 10.1214/009053607000000505 c Institute of Mathematical Statistics, 2007 MEASURING AND TESTING DEPENDENCE BY CORRELATION OF DISTANCES By Gabor´ J. Szekely,´ 1 Maria L. Rizzo and Nail K. Bakirov Bowling Green State University, Bowling Green State University and USC Russian Academy of Sciences Distance correlation is a new measure of dependence between ran- dom vectors. Distance covariance and distance correlation are anal- ogous to product-moment covariance and correlation, but unlike the classical definition of correlation, distance correlation is zero only if the random vectors are independent. The empirical distance depen- dence measures are based on certain Euclidean distances between sample elements rather than sample moments, yet have a compact representation analogous to the classical covariance and correlation. Asymptotic properties and applications in testing independence are discussed. Implementation of the test and Monte Carlo results are also presented. 1. Introduction. Distance correlation provides a new approach to the problem of testing the joint independence of random vectors. For all dis- tributions with finite first moments, distance correlation generalizes the idea of correlation in two fundamental ways: R (i) (X,Y ) is defined for X and Y in arbitrary dimensions; (ii) R(X,Y ) = 0 characterizes independence of X and Y . R Distance correlation has properties of a true dependence measure, analogous to product-moment correlation ρ. Distance correlation satisfies 0 1, and = 0 only if X and Y are independent. In the bivariate normal≤R≤ case, isR a function of ρ, and (X,Y ) ρ(X,Y ) with equality when ρ = 1.
    [Show full text]
  • The Perfect Marriage and Much More: Combining Dimension Reduction, Distance Measures and Covariance
    The Perfect Marriage and Much More: Combining Dimension Reduction, Distance Measures and Covariance Ravi Kashyap SolBridge International School of Business / City University of Hong Kong July 15, 2019 Keywords: Dimension Reduction; Johnson-Lindenstrauss; Bhattacharyya; Distance Measure; Covariance; Distribution; Uncertainty JEL Codes: C44 Statistical Decision Theory; C43 Index Numbers and Aggregation; G12 Asset Pricing Mathematics Subject Codes: 60E05 Distributions: general theory; 62-07 Data analysis; 62P20 Applications to economics Edited Version: Kashyap, R. (2019). The perfect marriage and much more: Combining dimension reduction, distance measures and covariance. Physica A: Statistical Mechanics and its Applications, XXX, XXX-XXX. Contents 1 Abstract 4 2 Introduction 4 3 Literature Review of Methodological Fundamentals 5 3.1 Bhattacharyya Distance . .5 3.2 Dimension Reduction . .9 4 Intuition for Dimension Reduction 10 arXiv:1603.09060v6 [q-fin.TR] 12 Jul 2019 4.1 Game of Darts . 10 4.2 The Merits and Limits of Four Physical Dimensions . 11 5 Methodological Innovations 11 5.1 Normal Log-Normal Mixture . 11 5.2 Normal Normal Product . 13 5.3 Truncated Normal / Multivariate Normal . 15 5.4 Discrete Multivariate Distribution . 19 5.5 Covariance and Distance . 20 1 6 Market-structure, Microstructure and Other Applications 24 6.1 Asset Pricing Application . 25 6.2 Biological Application . 26 7 Empirical Illustrations 27 7.1 Comparison of Security Prices across Markets . 27 7.2 Comparison of Security Trading Volumes, High-Low-Open-Close Prices and Volume-Price Volatilities . 30 7.3 Speaking Volumes Of: Comparison of Trading Volumes . 30 7.4 A Pricey Prescription: Comparison of Prices (Open, Close, High and Low) .
    [Show full text]
  • Distance Correlation Coefficient
    Journal of Modern Applied Statistical Methods Volume 13 | Issue 1 Article 23 5-1-2014 Distance Correlation Coefficient: An Application with Bayesian Approach in Clinical Data Analysis Atanu Bhattacharjee Malabar Cancer Centre, Kerala, India, [email protected] Follow this and additional works at: http://digitalcommons.wayne.edu/jmasm Part of the Applied Statistics Commons, Social and Behavioral Sciences Commons, and the Statistical Theory Commons Recommended Citation Bhattacharjee, Atanu (2014) "Distance Correlation Coefficient: An Application with Bayesian Approach in Clinical Data Analysis," Journal of Modern Applied Statistical Methods: Vol. 13 : Iss. 1 , Article 23. DOI: 10.22237/jmasm/1398918120 Available at: http://digitalcommons.wayne.edu/jmasm/vol13/iss1/23 This Regular Article is brought to you for free and open access by the Open Access Journals at DigitalCommons@WayneState. It has been accepted for inclusion in Journal of Modern Applied Statistical Methods by an authorized editor of DigitalCommons@WayneState. Journal of Modern Applied Statistical Methods Copyright © 2014 JMASM, Inc. May 2014, Vol. 13, No. 1, 354-366. ISSN 1538 − 9472 Distance Correlation Coefficient: An Application with Bayesian Approach in Clinical Data Analysis Atanu Bhattacharjee Malabar Cancer Centre Kerala, India The distance correlation coefficient – based on the product-moment approach – is one method by which to explore the relationship between variables. The Bayesian approach is a powerful tool to determine statistical inferences with credible intervals. Prior information about the relationship between BP and Serum cholesterol was applied to formulate the distance correlation between the two variables. The conjugate prior is considered to formulate the posterior estimates of the distance correlations. The illustrated method is simple and is suitable for other experimental studies.
    [Show full text]
  • Package 'Multivariance'
    Package ‘multivariance’ December 10, 2020 Title Measuring Multivariate Dependence Using Distance Multivariance Version 2.4.0 Date 2020-12-10 Description Distance multivariance is a measure of dependence which can be used to detect and quantify dependence of arbitrarily many random vectors. The necessary functions are implemented in this packages and examples are given. It includes: distance multivariance, distance multicorrelation, dependence structure detection, tests of independence and copula versions of distance multivariance based on the Monte Carlo empirical transform. Detailed references are given in the package description, as starting point for the theoretic background we refer to: B. Böttcher, Dependence and Dependence Structures: Estimation and Visualization Using the Unifying Concept of Distance Multivariance. Open Statistics, Vol. 1, No. 1 (2020), <doi:10.1515/stat-2020-0001>. Depends R (>= 3.3.0) License GPL-3 Encoding UTF-8 LazyData true Imports igraph, graphics, stats, Rcpp, microbenchmark RoxygenNote 7.1.1 Suggests testthat LinkingTo Rcpp NeedsCompilation yes Author Björn Böttcher [aut, cre], Martin Keller-Ressel [ctb] Maintainer Björn Böttcher <[email protected]> Repository CRAN Date/Publication 2020-12-10 15:50:03 UTC 1 2 R topics documented: R topics documented: multivariance-package . .3 anscombe.extended . .5 cdm .............................................6 cdms.............................................7 clean.graph . .8 coins . .9 copula.multicorrelation . 10 copula.multicorrelation.test . 11 copula.multivariance . 11 dependence.structure . 12 dep_struct_iterated_13_100 . 15 dep_struct_ring_15_100 . 16 dep_struct_several_26_100 . 17 dep_struct_star_9_100 . 17 emp.transf . 18 fastdist . 18 fastEuclideanCdm . 19 find.cluster . 19 independence.test . 21 layout_on_circles . 22 m.multivariance . 23 multicorrelation . 24 multivariance . 27 multivariance.pvalue . 28 multivariance.test . 29 multivariance.timing . 31 multivariances.all .
    [Show full text]
  • Causal Inference on Discrete Data Via Estimating Distance Correlations
    1 Causal Inference on Discrete Data via Estimating Distance Correlations Furui Liu, Laiwan Chan Department of Computer Science and Engineering, The Chinese University of Hong Kong Abstract In this paper, we deal with the problem of inferring causal directions when the data is on discrete domain. By considering the distribution of the cause P (X) and the condi- tional distribution mapping cause to effect P (Y jX) as independent random variables, we propose to infer the causal direction via comparing the distance correlation between P (X) and P (Y jX) with the distance correlation between P (Y ) and P (XjY ). We in- fer “X causes Y ” if the dependence coefficient between P (X) and P (Y jX) is smaller. Experiments are performed to show the performance of the proposed method. Editor: Aapo Hyvarinen¨ 1 Introduction Inferring the causal direction between two variables from observational data becomes arXiv:1803.07712v3 [stat.ML] 7 Aug 2018 a hot research topic. Additive noise models (ANMs) (Zhang and Hyvarinen,¨ 2009a; Shimizu et al., 2006; Hoyer et al., 2009; Peters et al., 2011; Hoyer and Hyttinen, 2009; Zhang and Hyvarinen,¨ 2009b; Shimizu et al., 2011; Hyvarinen¨ and Smith, 2013; Hyvarinen¨ et al., 2010; Hyttinen et al., 2012; Mooij et al., 2011) are preliminary trials to solve this problem. They assume that the effect is governed by the cause and an ad- ditive noise, and the causal inference is done by finding the direction that admits such a model. Recently, under another view of exploiting the asymmetry between cause and effect, the linear trace method (LTr) (Zscheischler et al., 2011; Janzing et al., 2010) and information geometric causal inference (IGCI) (Janzing et al., 2012) are proposed.
    [Show full text]
  • Lecture 2: Measures of Correlation and Dependence 4Mm Foundations of Data Science
    1 Lecture 2: Measures of Correlation and Dependence Foundations of Data Science: Algorithms and Mathematical Foundations Mihai Cucuringu [email protected] CDT in Mathematics of Random System University of Oxford 2 Measures of Correlation/Dependence I Pearson I Spearman I Kendall I Hoeffding’s I Maximal Correlation I Distance Correlation I Mutual Information I Maximal Information Coefficient (MIC) 3 Review of correlation measures Pearson correlation ρ is a measure of linear dependence between variables. I In the population: given random variables X; Y 2 R Cov(X; Y ) Cor(X; Y ) = pVar(X)pVar(Y ) I In the sample: given vectors x; y 2 R cov(x; y) (x − x¯1)T (y − y¯1) cor(x; y) = p p = var(x) var(y) jjx − x¯1jj2jjy − y¯1jj2 If x; y are have been centered xT y cor(x; y) = jjxjj2jjyjj2 4 Properties of the population correlation Cov(X; Y ) ρ =def Cor(X; Y ) = pVar(X)pVar(Y ) Properties of ρ I Cor(X; X) = 1 I Cor(X; Y ) = Cor(Y ; X) I Cor(aX + b; Y ) = sign(a)Cor(X; Y ) for any a; b 2 R I −1 ≤ Cor(X; Y ) ≤ 1 I jCor(X; Y )j = 1 if and only if Y = aX + b for some a; b 2 R, with a 6= 0 I If X; Y are independent then Cor(X; Y ) = 0 I If Cor(X; Y ) = 0 then X; Y need not be independent!!! I If (X; Y ) is bivariate normal and Cor(X; Y ) = 0,then X; Y are independent 5 Bivariate normal distribution Two-dimensional Gaussian distribution The random vector Z = (X; Y ) 2 R2 has a bivariate normal dist.
    [Show full text]