Check Markov equivalence of two DAGs {MXM}

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.

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.

Value

A list including:

apofasi

A boolean variable, TRUE of FALSE.

mes

A message specyfying the result, the dimensions of the adjacency matrices do not match for example, or the number of adjancencies is not the same, they do not share the same unshilded colliders, or they are Markov equivalent.

Author(s)

Michail Tsagris

R implementation and documentation: Giorgos Athineou <[e-mail]> and Michail Tsagris [e-mail]

References

Tsamardinos, Brown and Aliferis (2006). The max-min hill-climbing Bayesian network structure learning algorithm. Machine learning, 65(1), 31-78.

See Also

pc.or, mmhc.skel, pc.con

Examples

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)