udag2pdag {pcalg} | R Documentation |
This function extends a pcAlgo-object containing a skeleton and corresponding conditional independence information to a Partially Directed Acyclic Graph (PDAG). The result is a pcAlgo-object as well.
udag2pdag(gInput,verbose)
gInput |
pcAlgo-object containing skeleton and cond. ind. information |
verbose |
0: No output; 1: Details |
The skeleton is extended to a PDAG using rules by Spirtes and Pearl (see References).
pcObj |
Oriented pc-Object |
Markus Kalisch (kalisch@stat.math.ethz.ch)
P. Spirtes, C. Glymour and R. Scheines (2000) Causation, Prediction, and Search, 2nd edition, The MIT Press.
J. Pearl (2000), Causality, Cambridge University Press.
pdag2dag
,dag2cpdag
,udag2pdag
,udag2pdagRelaxed
,udag2pdagSpecial
## simulate data set.seed(123) p <- 10 myDAG <- randomDAG(p, prob = 0.2) trueCPDAG <- dag2cpdag(myDAG) n <- 1000 d.mat <- rmvDAG(n, myDAG, errDist = "normal") ## estimate skeleton resU <- pcAlgo(d.mat, alpha = 0.05, corMethod = "standard",directed=FALSE) ## orient edges using three different methods resD1 <- udag2pdagRelaxed(resU,verbose=0) resD2 <- udag2pdagSpecial(resU,verbose=0,n.max=100) resD3 <- udag2pdag(resU,verbose=0)