library(OpenMx) library(lme4) (fm1 <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy)) summary(fm1) data(demoOneFactor) manifests <- names(demoOneFactor) latents <- c("G") factorModel <- mxModel( "One Factor", mxMatrix("Full", 5, 1, values=0.8, free=TRUE, name="A"), mxMatrix("Symm", 1, 1, values=1, free=FALSE, name="L"), mxMatrix("Diag", 5, 5, values=1, free=TRUE, name="U"), mxAlgebra(A %*% L %*% t(A) + U, name="R"), mxExpectationNormal(covariance = "R", dimnames = names(demoOneFactor)), mxFitFunctionML(), mxData(cov(demoOneFactor), type="cov", numObs=500)) factorModelFit <- mxRun(factorModel) summary(factorModelFit) save.image("loadtest181129b.RData")