Copyright © 2007-2024 The OpenMx Project
require(OpenMx) data(demoOneFactor) latents = c("G") manifests = names(demoOneFactor) # 1. Make two models, and nest them in a multigroup m1 <- mxModel("model1", type = "RAM", manifestVars = manifests, latentVars = latents, mxPath(from = latents, to = manifests), mxPath(from = manifests, arrows = 2), mxPath(from = latents, arrows = 2, free = F, values = 1.0), mxData(cov(demoOneFactor), type = "cov", numObs = 500) ) m2 <- m1 m2 <- mxRename(m2, 'model2') # 2. Nest them in a multigroup supermodel and run m3 = mxModel("bob", m1, m2, mxFitFunctionMultigroup(c("model1","model2")) ) m3Run = mxRun(m3); # 3 request saturated model x = omxSaturatedModel(m3Run) Error in obsdata[selVars, selVars] : subscript out of bounds traceback() # 2: omxSaturatedModel(x[[grpnames[i]]], run = FALSE) # 1: omxSaturatedModel(m3)