## Not run:
# library(OpenMx)
# data(demoOneFactor)
# manifests <- names(demoOneFactor)
# latents <- c("G")
# factorModel <- mxModel("One Factor",
# type="RAM",
# manifestVars=manifests,
# latentVars=latents,
# mxPath(from=latents, to=manifests),
# mxPath(from=manifests, arrows=2),
# mxPath(from=latents, arrows=2, free=FALSE, values=1.0),
# mxData(observed=cov(demoOneFactor), type="cov", numObs=500)
# )
# factorFit <- mxRun(factorModel)
# round(fitMeasuresMx(factorFit), 3)
#
# # Compare with lavaan
# library(lavaan)
# script <- "f1 =~ x1 + x2 + x3 + x4 + x5"
# fitMeasures(cfa(script, sample.cov = cov(demoOneFactor), sample.nobs = 500, std.lv = TRUE))
# ## End(Not run)
Run the code above in your browser using DataLab