Causal Discovery from Discrete Data using Hidden Compact Representation

This code provides a method to fit the hidden compact representation model as well as to identify the causal direction on discrete data. We implement an effective solution to recover the above hidden compact representation under the likelihood framework. Please see the Causal Discovery from Discrete Data using Hidden Compact Representation from NIPS 2018 by Ruichu Cai, Jie Qiao, Kun Zhang, Zhenjie Zhang and Zhifeng Hao (2018) < https://nips.cc/Conferences/2018/Schedule?showEvent=11274> for a description of some of our methods.


Overview

This code provides a method to fit the hidden compact representation model as well as to discover the causal direction on discrete data.

Please cite "Ruichu Cai, Jie Qiao, Kun Zhang, Zhenjie Zhang, Zhifeng Hao. Causal Discovery from Discrete Data using Hidden Compact Representation. NIPS 2018."

Installation

install.packages("HCR")

Quick Start

This package contains the data synthetic process for HCR model. Here are some examples to make a quick start:

# HCR data
library(data.table)
set.seed(1)
data=simuXY(sample_size=2000)
r1<-HCR(data$X,data$Y,score_type = "aic")
r2<-HCR(data$Y,data$X,score_type = "aic")
if(r1$score>r2$score){
  print("X cause Y.")
}else{
  print("Y cause X.")
}
unique(r1$data[,c("X","Yp")]) # The fitted hidden representation
unique(data.frame(data$X,data$Yp)) # The true hidden representation

# A fast implementation of HCR
r1<-HCR.fast(data$X,data$Y)
r2<-HCR.fast(data$Y,data$X)
if(r1$score>r2$score){
  print("X cause Y.")
}else{
  print("Y cause X.")
}
unique(r1$data[,c("X","Yp")]) # The fitted hidden representation
unique(data.frame(data$X,data$Yp)) # The true hidden representation

# ANM data
# Note that the complete generate of the ANM data can be found in http://webdav.tuebingen.mpg.de/causality/
set.seed(0)
X=sample(1:5,2000,replace = T)
N=sample(-1:1,2000,replace = T)
Y=X+N
r1=HCR(X,Y,score_type = "bic",is_anm = T)
r2=HCR(Y,X,score_type = "bic",is_anm = T)
if(r1$score>r2$score){
  print("X cause Y.")
}else{
  print("Y cause X.")
}


News

HCR 0.1.1

  • Fix an error in the score function.

Reference manual

It appears you don't have a PDF plugin for this browser. You can click here to download the reference manual.

install.packages("HCR")

0.1.1 by Jie Qiao, 8 months ago


Browse source code at https://github.com/cran/HCR


Authors: Jie Qiao [aut, cre] , Ruichu Cai [ths, aut] , Kun Zhang [ths, aut] , Zhenjie Zhang [ths, aut] , Zhifeng Hao [ths, aut]


Documentation:   PDF Manual  


GPL (>= 2) license


Imports data.table, methods


See at CRAN