umxCI: umxCI

Description

umxCI adds mxCI() calls for all free parameters in a model, runs the CIs, and reports a neat summary.

Usage

umxCI(model = NULL, which = c("ALL", NA, "list of your making"), remove = FALSE, run = c("no", "yes", "if necessary"), showErrorCodes = TRUE)

Arguments

model

The mxModel you wish to report mxCIs on

which

What CIs to add: c("ALL", NA, "list of your making")

remove

= FALSE (if set, removes existing specified CIs from the model)

run

Whether or not to compute the CIs. Valid values = "no" (default), "yes", "if necessary".

showErrorCodes

Whether to show errors (default == TRUE)

Value

- mxModel

Details

This function also reports any problems computing a CI. The codes are standard OpenMx errors and warnings

1: The final iterate satisfies the optimality conditions to the accuracy requested, but the sequence of iterates has not yet converged. NPSOL was terminated because no further improvement could be made in the merit function (Mx status GREEN)
2: The linear constraints and bounds could not be satisfied. The problem has no feasible solution.
3: The nonlinear constraints and bounds could not be satisfied. The problem may have no feasible solution.
4: The major iteration limit was reached (Mx status BLUE).
6: The model does not satisfy the first-order optimality conditions to the required accuracy, and no improved point for the merit function could be found during the final linesearch (Mx status RED)
7: The function derivates returned by funcon or funobj appear to be incorrect.
9: An input parameter was invalid

If runCIs is FALSE, the function simply adds CIs to be computed and returns the model.

References

- http://www.github.com/tbates/umx/

Examples

Run this code

require(umx)
data(demoOneFactor)
latents  = c("G")
manifests = names(demoOneFactor)
m1 <- 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(cov(demoOneFactor), type = "cov", numObs = 500)
)
m1 = umxRun(m1, setLabels = TRUE, setValues = TRUE)
m1$intervals # none yet list()
m1 = umxCI(m1)
m1$intervals # $G_to_x1
m1 = umxCI(m1, remove = TRUE) # Add CIs for all free parameters, and return model
## Not run: 
# umxCI(model, run = "yes") # force update of CIs
# # Don't force update of CIs, but if they were just added, then calculate them
# umxCI(model, run = "if necessary")
# ## End(Not run)

Run the code above in your browser using DataLab

State of Data and AI Literacy Report 2025