# NOT RUN {
# Load dataset:
library("lavaan")
data(HolzingerSwineford1939)
Data <- HolzingerSwineford1939[,7:15]
# Measurement model:
Lambda <- matrix(0, 9, 3)
Lambda[1:3,1] <- NA
Lambda[4:6,2] <- NA
Lambda[7:9,3] <- NA
# Search best fitting omega_theta:
# }
# NOT RUN {
res <- lvnetLasso(Data, "omega_theta", lambda = Lambda)
res$best
summary(res)
# }
Run the code above in your browser using DataLab