Comparison of strategies for scalable causal discovery of latent variable models from mixed data

• Published in: International journal of data science and analytics 2018, Vol.6 (1), p.33-45
• Main Authors: Raghu, Vineet K., Ramsey, Joseph D., Morris, Alison, Manatakis, Dimitrios V., Sprites, Peter, Chrysanthis, Panos K., Glymour, Clark, Benos, Panayiotis V.
• Format: Article
• Language: English
• Subjects: Artificial Intelligence; Business Information Systems; Computational Biology/Bioinformatics; Computer Science; Data Mining and Knowledge Discovery; Database Management; Regular Paper
• ISSN: 2364-415X; 2364-4168
• DOI: 10.1007/s41060-018-0104-3

Modern technologies allow large, complex biomedical datasets to be collected from patient cohorts. These datasets are comprised of both continuous and categorical data (“Mixed Data”), and essential variables may be unobserved in this data due to the complex nature of biomedical phenomena. Causal inference algorithms can identify important relationships from biomedical data; however, handling the challenges of causal inference over mixed data with unmeasured confounders in a scalable way is still an open problem. Despite recent advances into causal discovery strategies that could potentially handle these challenges; individually, no study currently exists that comprehensively compares these approaches in this setting. In this paper, we present a comparative study that addresses this problem by comparing the accuracy and efficiency of different strategies in large, mixed datasets with latent confounders. We experiment with two extensions of the Fast Causal Inference algorithm: a maximum probability search procedure we recently developed to identify causal orientations more accurately, and a strategy which quickly eliminates unlikely adjacencies in order to achieve scalability to high-dimensional data. We demonstrate that these methods significantly outperform the state of the art in the field by achieving both accurate edge orientations and tractable running time in simulation experiments on datasets with up to 500 variables. Finally, we demonstrate the usability of the best performing approach on real data by applying it to a biomedical dataset of HIV-infected individuals.