ggm (version 2.5.1)

MarkEqRcg: Markov equivalence for regression chain graphs.

Description

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

Usage

MarkEqRcg(amat, bmat)

Value

"Markov Equivalent" or "NOT Markov Equivalent".

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.

Author

Kayvan Sadeghi

Details

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

References

Wermuth, N. and Sadeghi, K. (2012). Sequences of regressions and their independences. Test 21:215–252.

See Also

MarkEqMag, msep

Examples

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)