Runs PC algorithm given a dataset dat.
pc_search(dat, alpha = 0.05)Square binary matrix (or NA in case of error).
The CPDAG estimated from the data.
The function pcalg:pc() is called with the following
arguments:
suffStat = list(C = cor(dat), n = NROW(dat)),
indepTest = gaussCItest,
u2pd = "retry" — this ensures that the produced CPDAG
is extendible to a DAG,
skel.method = "stable".