Markov equivalence for regression chain graphs.
MarkEqMag
determines whether two RCGs (or subclasses of RCGs) are
Markov equivalent.
MarkEqRcg(amat, bmat)
amat |
An adjacency matrix of an RCG or a graph that can be a |
bmat |
The same as |
The function checks whether the two graphs have the same skeleton and unshielded colliders.
"Markov Equivalent" or "NOT Markov Equivalent".
Kayvan Sadeghi
Wermuth, N. and Sadeghi, K. (2011). Sequences of regressions and their independences. TEST, To appear. http://arxiv.org/abs/1103.2523.
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)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.