# Load data
wmt <- wmt2[,7:24]
if (FALSE) # Estimate EGA network
## plot.type = "qqraph" used for CRAN checks
## plot.type = "GGally" is the default
ega.wmt <- EGA(data = wmt, model = "glasso", plot.type = "qgraph")
# Estimate dimension stability
boot.wmt <- bootEGA(data = wmt, iter = 100, typicalStructure = TRUE,
plot.typicalStructure = TRUE, model = "glasso", plot.type = "qgraph",
type = "parametric", ncores = 2)
# Estimate item stability statistics
res <- itemStability(boot.wmt)
# Changing plot features (ggplot2)
## Changing colors (ignore warnings)
### qgraph Defaults
res$plot +
ggplot2::scale_color_manual(values = rainbow(max(res$membership$unique)))
### Pastel
res$plot +
ggplot2::scale_color_brewer(palette = "Pastel1")
## Changing Legend (ignore warnings)
res$plot +
ggplot2::scale_color_discrete(labels = "Intelligence")
Run the code above in your browser using DataLab