# Load data
dep <- depression[,24:44]
# Estimate EGA
## plot.type = "qqraph" used for CRAN checks
## plot.type = "GGally" is the default
ega.dep <- EGA(data = dep, model = "glasso", plot.type = "qgraph")
# Compute entropy indices
vn.entropy(data = ega.dep$correlation, structure = ega.dep$wc)
Run the code above in your browser using DataLab