data(learning.test)
res = gs(learning.test)
# the arc between E and F has no direction
plot(res)
res = choose.direction(res, c("E", "F"), learning.test, debug = TRUE)
# * testing E - F for direction.
# > testing E -> F with conditioning set ' '.
# > p-value is 2.174867e-197 .
# > testing F -> E with conditioning set ' B '.
# > p-value is 3.935648e-245 .
# @ removing E -> F .
plot(res, highlight = c("E", "F"))Run the code above in your browser using DataLab