skeleton
, pc
, fci
.gaussCItest(x, y, S, suffStat)
zStat
for details.dsepTest
, disCItest
and
binCItest
for similar functions for a d-separation
oracle, conditional independence test on discrete variables and
conditional independence test on binary variables, respectively.## simulate data: x -> y -> z
set.seed(29)
x <- rnorm(100)
y <- 3*x + rnorm(100)
z <- 2*y + rnorm(100)
dat <- cbind(x,y,z)
## analyze data
suffStat <- list(C = cor(dat), n = nrow(dat))
gaussCItest(1,3,NULL,suffStat) ## dependent
gaussCItest(1,3,2,suffStat) ## independent
Run the code above in your browser using DataLab