data(learning.test)
res = gs(learning.test)
# the arc between E and F has no direction
plot(res, highlight = c("E", "F"))
res2 = choose.direction(res, c("E", "F"), learning.test)
plot(res, highlight = c("E", "F"))
compare(res, res2, debug = TRUE)
# * arcs in r1 not present in r2:
# [1] "F -> E"
# * arcs in r2 not present in r1:
# character(0)
# [1] FALSE
Run the code above in your browser using DataLab