showClass("pcAlgo")
## generate a pcAlgo object
p <- 8
set.seed(45)
myDAG <- randomDAG(p, prob = 0.3)
n <- 10000
d.mat <- rmvDAG(n, myDAG, errDist = "normal")
indepTest <- gaussCItest
suffStat <- list(C = cor(d.mat), n = n)
alpha <- 0.01
pc.fit <- pc(suffStat, indepTest, p, alpha)
## use methods of class pcAlgo
show(pc.fit)
if(require(Rgraphviz))
plot(pc.fit)
summary(pc.fit)
## access slots of this object
(g <- pc.fit@graph)
str(ss <- pc.fit@sepset, max=1)
Run the code above in your browser using DataLab