# NOT RUN {
Y <- BGGM::bfi[, 1:20]
# analytic approach (sample by setting analytic = FALSE)
fit_analytic <- estimate(Y, analytic = TRUE)
# select the graph (edge set E)
E <- select(fit_analytic, ci_width = 0.95)
# plot
plt <- plot(E,
node_labels = letters[1:20],
node_labels_color = "white",
node_groups = rep(c("1", "2", "3", "4"), each = 5),
edge_colors = "classic",
alpha = 0.5, palette = "Set2")
# }
Run the code above in your browser using DataLab