ggm (version 2.3)

MarkEqRcg: Markov equivalence for regression chain graphs.

Description

MarkEqMag determines whether two RCGs (or subclasses of RCGs) are Markov equivalent.

Usage

MarkEqRcg(amat, bmat)

Arguments

amat

An adjacency matrix of an RCG or a graph that can be a graphNEL or an igraph object or a vector of length \(3e\), where \(e\) is the number of edges of the graph, that is a sequence of triples (type, node1label, node2label). The type of edge can be "a" (arrows from node1 to node2), "b" (arcs), and "l" (lines).

bmat

The same as amat.

Value

"Markov Equivalent" or "NOT Markov Equivalent".

Details

The function checks whether the two graphs have the same skeleton and unshielded colliders.

References

Wermuth, N. and Sadeghi, K. (2011). Sequences of regressions and their independences. TEST, To appear. http://arxiv.org/abs/1103.2523.

See Also

MarkEqMag, msep

Examples

# NOT RUN {
H1<-matrix(c(0,100,0,0,0,100,0,100,0,0,0,100,0,0,0,1,0,0,0,100,0,0,1,100,0),5,5)
H2<-matrix(c(0,0,0,0,0,1,0,1,0,0,0,0,0,0,0,1,0,0,0,100,0,0,1,100,0),5,5)
H3<-matrix(c(0,0,0,0,0,1,0,1,0,0,0,0,0,0,0,1,0,0,0,0,0,0,1,1,0),5,5)
#MarkEqRcg(H1,H2)
#MarkEqRcg(H1,H3)
#MarkEqRcg(H2,H3)
# }