# NOT RUN {
# BFI Extraversion data from psychTools package:
library("psychTools")
data(bfi)
# Subset of data:
bfiSub <- bfi[1:250,1:25]
# Estimate ggmModSelect networks (not stepwise to increase speed):
Network <- estimateNetwork(bfiSub], default = "ggmModSelect", corMethod = "cor",
stepwise = FALSE)
# Bootstrap 100 values, using 8 cores (100 to incease speed, preferably 1000+):
boots <- bootnet(Network, nBoots = 100, nCores = 8)
# Threshold network:
Network_inclusion <- bootInclude(boots)
# Plot:
plot(Network_inclusion)
# }
Run the code above in your browser using DataLab