MarkEqRcg {ggm}R Documentation

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 3e3e, where ee 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.

Details

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

Value

"Markov Equivalent" or "NOT Markov Equivalent".

Author(s)

Kayvan Sadeghi

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)

[Package ggm version 2.5.1 Index]