library(easybgm)
library(bgms)
data <- na.omit(Wenchuan)
fit <- easybgm(data, type = "continuous",
iter = 1000 # for demonstration only (> 5e4 recommended)
)
plot_network(fit)
# Shows all edges with an inclusion probability larger than 0.1
plot_network(fit, exc_prob = 0.1)
# Indicate which edges have insufficient evidence for inclusion through a dashed line
plot_network(fit, dashed = TRUE, evidence_thresh = 10)
Run the code above in your browser using DataLab