## Not run:
# ### Using bfi dataset from psych ###
# library("psych")
# data(bfi)
#
# # Compute correlations:
# CorMat <- cor_auto(bfi[,1:25])
#
# # Compute graph with tuning = 0 (BIC):
# BICgraph <- EBICglasso(CorMat, nrow(bfi), 0)
#
# # Compute graph with tuning = 0.5 (EBIC)
# EBICgraph <- EBICglasso(CorMat, nrow(bfi), 0.5)
#
# # Plot both:
# layout(t(1:2))
# BICgraph <- qgraph(BICgraph, layout = "spring", title = "BIC", details = TRUE)
# EBICgraph <- qgraph(EBICgraph, layout = "spring", title = "EBIC")
#
# # Compare centrality and clustering:
# layout(1)
# centralityPlot(list(BIC = BICgraph, EBIC = EBICgraph))
# clusteringPlot(list(BIC = BICgraph, EBIC = EBICgraph))
# ## End(Not run)
Run the code above in your browser using DataLab