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.041687e-44 .
# > testing F -> E with conditioning set ' '.
# > p-value is 4.598853e-08 .
# @ removing F -> E .
plot(res, highlight = c("E", "F"))
Run the code above in your browser using DataLab