summary-MxModel {OpenMx}R Documentation

Model Summary

Description

This function returns summary statistics of a model after it has been run

Usage

summary(object, ..., verbose=FALSE)

Arguments

object

A MxModel object.

...

Any number of named arguments (see below).

verbose

logical. Changes the printing style for summary (see Details)

Details

mxSummary allows the user to set or override the following parameters of the model:

numObs

Numeric. Specify the total number of observations for the model.

numStats

Numeric. Specify the total number of observed statistics for the model.

refModels

List of MxModel objects. Specify a saturated and independence likelihoods in single argument for testing.

SaturatedLikelihood

Numeric or MxModel object. Specify a saturated likelihood for testing.

SaturatedDoF

Numeric. When SaturatedLikelihood is numeric, specify the degrees of freedom of the saturated likelihood for testing.

IndependenceLikelihood

Numeric or MxModel object. Specify an independence likelihood for testing.

IndependenceDoF

Numeric. When IndependenceLikelihood is numeric, specify the degrees of freedom of the independence likelihood for testing.

indep

Logical. Set to FALSE to ignore independent submodels in summary.

The verbose argument changes the printing style for the summary of a model. When verbose=FALSE, a relatively minimal amount of information is printed: the free parameters, the likelihood, and a few fit indices. When more information is available, more is printed. For example, when the model has a saturated likelihood, several additional fit indices are printed. On the other hand, when verbose=TRUE, the compute plan, the data summary, and additional timing information are always printed. Moreover, available fit indices are printed regarless of whether or not they are defined. The undefined fit indices are printed as NA. Running a saturated model and including it with the call to summary will define these fit indices and they will dislay meaningful values. It should be noted that the verbose argument only changes the printing style, all of the same information is calculated and exists in the output of summary. More information is displayed when verbose=TRUE, and less when verbose=FALSE.

The Information Criteria (AIC, BIC) are reported in a table. The table shows different versions of the information criteria. Each entry in the table is an AIC or BIC obtained using different penalties. In particular, the entries of the table do not show different penalties, but rather different versions of AIC and BIC. For example the AIC is reported with both a Parameters Penalty and a Degrees of Freedom Penalty. AIC generally takes the form -2LL + 2*k. With the Parameters Penalty k is the number of free parameters. With the Degrees of Freedom Penalty, k is the model degrees of freedom. BIC is defined similarly: -2LL + k*log(N) where k is either the number of free parameters or the degrees of freedom. The Sample-Size Adjusted BIC is only defined for the parameters penalty: -2LL + k*log((N+2)/24).

The refModels, SaturatedLikelihood, SaturatedDoF, IndependenceLikelihood, and IndependenceDoF arguments can be used to obtain further fit statistics (RMSEA, CFI, TLI, Chi-Squared). For covariance data, saturated and independence models are fitted automatically so all fit indices are reported. For raw data, these reference models are not estimated to save computational time. For an easy way to make reference models for most cases is provided by the mxRefModels function.

This function can report Error codes as follows:

References

The OpenMx User's guide can be found at http://openmx.psyc.virginia.edu/documentation.

Examples


library(OpenMx)
data(demoOneFactor)  # load the demoOneFactor dataframe
manifests <- names(demoOneFactor) # set the manifest to the 5 demo variables
latents <- c("G")  # define 1 latent variable
model <- mxModel(model="One Factor", type="RAM",
    manifestVars = manifests,
    latentVars = latents,
    mxPath(from = latents, to=manifests, labels = paste("b", 1:5, sep = "")),
    mxPath(from = manifests, arrows = 2, labels = paste("u", 1:5, sep = "")),
    mxPath(from = latents, arrows = 2, free = FALSE, values = 1.0),
    mxData(cov(demoOneFactor), type = "cov", numObs = 500)
)
model <- mxRun(model) # Run the model, returning the result into model

# Show summary of the fitted model
summary(model)

# Compute the summary and store in the variable "statistics"
statistics <- summary(model)

# Access components of the summary
statistics$parameters
statistics$SaturatedLikelihood

# Specify a saturated likelihood for testing
summary(model, SaturatedLikelihood = -3000)

# Add a CI and view it in the summary
model = mxRun(mxModel(model=model, mxCI("b5")), intervals = TRUE)
summary(model)


[Package OpenMx version 2.0.0-3948 Index]