data(learning.test)
cpdag = pc.stable(learning.test)
plot(cpdag)
## highlight node B and related arcs.
plot(cpdag, highlight = "B")
## highlight B and its Markov blanket.
plot(cpdag, highlight = c("B", mb(cpdag, "B")))
## a more compact plot.
par(oma = rep(0, 4), mar = rep(0, 4), mai = rep(0, 4),
plt = c(0.06, 0.94, 0.12, 0.88))
plot(cpdag)
Run the code above in your browser using DataLab