umxSummary.MxModel: Shows a compact, publication-style, summary of a RAM model

Description

Report the fit of a model in a compact form suitable for a journal. Emits a "warning" when model fit is worse than accepted criterion (TLI >= .95 and RMSEA <= .06; (Hu & Bentler, 1999; Yu, 2002).

Usage

# S3 method for MxModel
umxSummary(model, refModels = NULL,
  showEstimates = c("raw", "std", "none", "both"), digits = 2,
  report = c("1", "table", "html"), filter = c("ALL", "NS", "SIG"),
  SE = TRUE, RMSEA_CI = FALSE, matrixAddresses = FALSE, std = NULL, ...)

Arguments

model

The mxModel whose fit will be reported

refModels

Saturated models if needed for fit indices (see example below: If NULL will be competed on demand. If FALSE will not be computed. Only needed for raw data.

showEstimates

What estimates to show. By default, the raw estimates are shown (Options = c("raw", "std", "none", "both").

digits

How many decimal places to report (default = 2)

report

If "html", then show results in browser ("1", "table", "html")

filter

whether to show significant paths (SIG) or NS paths (NS) or all paths (ALL)

Whether to compute SEs... defaults to TRUE. In rare cases, you might need to turn off to avoid errors.

RMSEA_CI

Whether to compute the CI on RMSEA (Defaults to FALSE)

matrixAddresses

Whether to show "matrix address" columns (Default = FALSE)

std

deprecated: use show = "std" instead!

...

Other parameters to control model summary

Value

- parameterTable returned invisibly, if estimates requested

Details

Note: For some (multi-group) models, you will need to fall back on summary

CIs and Identification This function uses the standard errors reported by OpenMx to produce the CIs you see in umxSummary These are used to derive confidence intervals based on the formula 95

Sometimes they appear NA. This often indicates a model which is not identified (seehttp://davidakenny.net/cm/identify.htm). This can include empirical under-identification - for instance two factors that are essentially identical in structure. use mxCheckIdentification to check identification.

One or more paths estimated at or close to zero suggests that fixing one or two of these to zero may fix the standard error calculation, and alleviate the need to estimate likelihood-based or bootstrap CIs

If factor loadings can flip sign and provide identical fit, this creates another form of under-identification and can break confidence interval estimation. Fixing a factor loading to 1 and estimating factor variances can help here.

References

- Hu, L., & Bentler, P. M. (1999). Cutoff criteria for fit indexes in covariance structure analysis: Conventional criteria versus new alternatives. Structural Equation Modeling, 6, 1-55.

- Yu, C.Y. (2002). Evaluating cutoff criteria of model fit indices for latent variable models with binary and continuous outcomes. University of California, Los Angeles, Los Angeles. Retrieved from http://www.statmodel.com/download/Yudissertation.pdf

http://tbates.github.io

Examples

Run this code

# NOT RUN {
require(umx)
data(demoOneFactor)
latents  = c("G")
manifests = names(demoOneFactor)
m1 <- umxRAM("One Factor",
	data = mxData(cov(demoOneFactor), type = "cov", numObs = 500),
	umxPath(latents, to = manifests),
	umxPath(var = manifests),
	umxPath(var = latents, fixedAt = 1)
)
umxSummary(m1, showEstimates = "std")
# output as latex
umx_set_table_format("latex")
umxSummary(m1, showEstimates = "std")
umx_set_table_format("markdown")
# output as raw
umxSummary(m1, show = "raw")
m1 <- mxModel(m1,
  mxData(demoOneFactor[1:100,], type = "raw"),
  umxPath(mean = manifests),
  umxPath(mean = latents, fixedAt = 0)
)
m1 <- mxRun(m1)
umxSummary(m1, showEstimates = "std", filter = "NS")
# }

Run the code above in your browser using DataLab

Last chance! 50% off unlimited learning