udag2pdagRelaxed {pcalg}R Documentation

Extend a pcAlgo-object containing a skeleton to a PDAG

Description

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. There is no check whether the result is extendible to a DAG

Usage

udag2pdagRelaxed(gInput,verbose)

Arguments

gInput pcAlgo-object containing skeleton and cond. ind. information
verbose 0: No output; 1: Details

Details

The skeleton is extended to a PDAG using rules by Spirtes and Pearl (see References). There is no test whether the result is really extendible.

Value

pcObj Oriented pc-Object

Author(s)

Markus Kalisch (kalisch@stat.math.ethz.ch)

References

P. Spirtes, C. Glymour and R. Scheines (2000) Causation, Prediction, and Search, 2nd edition, The MIT Press.

J. Pearl (2000), Causality, Cambridge University Press.

See Also

pdag2dag,dag2cpdag,udag2pdag,udag2pdagRelaxed,udag2pdagSpecial

Examples

## 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)


[Package pcalg version 0.1-9 Index]