# NOT RUN {
loading <- matrix(0, 6, 1)
loading[1:6, 1] <- NA
LY <- bind(loading, 0.7)
RPS <- binds(diag(1))
RTE <- binds(diag(6))
CFA.Model <- model(LY = LY, RPS = RPS, RTE = RTE, modelType="CFA")
# Write a function to extract the modification index from lavaan object
outfun <- function(out) {
result <- inspect(out, "mi")
}
# We will use only 5 replications to save time.
# In reality, more replications are needed.
Output <- sim(5, n=200, model=CFA.Model, outfun=outfun)
# Get the modification index of each replication
getExtraOutput(Output)
# }
Run the code above in your browser using DataLab