umx (version 1.4.0)

confint.MxModel: Get confidence intervals from an MxModel

Description

Implements confidence interval function for OpenMx models. Note: Currently requested CIs are added to existing CIs, and all are run, even if they alrady exist in the output. This might change in the future.

Usage

"confint"(object, parm = list("existing", c("vector", "of", "names"), "default = add all"), level = 0.95, run = FALSE, showErrorCodes = FALSE, ...)

Arguments

object

An mxModel, possibly already containing mxCIs that have been mxRun with intervals = TRUE))

parm

A specification of which parameters are to be given confidence intervals. Can be "existing", "all", or a vector of names.

level

The confidence level required (default = .95)

run

Whether to run the model (defaults to FALSE)

showErrorCodes

(default = FALSE)

...

Additional argument(s) for umxConfint.

Value

- mxModel

Details

Unlike confint, if parm is missing, all CIs requested will be added to the model, but (because these can take time to run) by default only CIs already computed will be reported.

CIs will be run only if run is TRUE, allowing this function to be used to add CIs without automatically having to run them. If parm is empty, and run = FALSE, a message will alert you to add run = TRUE. Even a few CIs can take too long to make running the default.

References

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

See Also

- confint

Other Reporting functions: RMSEA.MxModel, RMSEA.summary.mxmodel, RMSEA, extractAIC.MxModel, loadings, logLik.MxModel, plot.MxModel, residuals.MxModel, umxCI_boot, umxCI, umxCompare, umxExpCov, umxExpMeans, umxFitIndices, umxPlotACEcov, umxPlotACE, umxPlotCP, umxPlotGxE, umxPlotIP, umxSummary.MxModel, umxSummaryACE, umx_drop_ok, umx_standardize_RAM

Examples

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)
m2 = confint(m1) # default: CIs added, but user prompted to set run = TRUE
m2 = confint(m2, run = TRUE) # CIs run and reported
# Add CIs for asymmetric paths in RAM model, report them, save m1 with this CI added
m1 = confint(m1, parm = "G_to_x1", run = TRUE) 
# Add CIs for asymmetric paths in RAM model, report them, save m1 with mxCIs added
m1 = confint(m1, parm = "A", run = TRUE)
confint(m1, parm = "existing") # request existing CIs (none added yet...)