# Load data
wmt <- wmt2[,7:24]
# Obtain correlation matrix
S <- auto.correlate(wmt)
# EBICglasso (default for EGA functions)
glasso_network <- network.estimation(
data = wmt, model = "glasso"
)
# Obtain fit (expects continuous variables!)
network.fit(network = glasso_network, n = nrow(wmt), S = S)
# Scaled metrics are not yet available for
# dichotomous or polytomous data!
Run the code above in your browser using DataLab