Check Markov equivalence of two DAGs:
Check Markov equivalence of two DAGs
Description
Check Markov equivalence of two DAGs.
Usage
equivdags(g1, g2)
Arguments
g1
The matrix of a DAG or a partially directed graph as produced from pc.or or any other algorithm.
g2
The matrix of a DAG or a partially directed graph as produced from pc.or or any other algorithm.
Value
A list including:
A list including:
Details
Two DAGs are Markov equivalent if a) they have the same adjancencies (regardlsee of the mark, arrowhead, tail or nothing) and b) they have the same unshielded colliders.
References
Tsamardinos, Brown and Aliferis (2006). The max-min hill-climbing Bayesian network structure learning algorithm. Machine learning, 65(1), 31-78.
y <- rdag(1000, 10, 0.3)
tru <- y$G
x <- y$x
mod <- pc.con(x)
eg <- dag2eg(y$G) ## make it essential graph first est <- pc.or(mod)$G
equivdags(est, tru)