## 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