Description Usage Arguments Details References Examples

This function creates and, optionally, runs saturated and independence (null) models of a base model or data set for use with mxSummary to enable fit indices that depend on these models.

1 | ```
mxRefModels(x, run=FALSE, ..., distribution="default")
``` |

`x` |
A MxModel object, data frame, or matrix. |

`run` |
logical. If TRUE, runs the models before returning; otherwise returns built models without running. |

`...` |
Not used. Forces remaining arguments to be specified by name. |

`distribution` |
character. Which distribution to assume. |

For typical structural equation models the saturated model is the free-est possible model. Not only all variances (and, when possible, all means) are estimated, but also all covariances. In the case of ordinal data, the ordinal means are fixed to zero and the thresholds are estimated. For binary variables, those variances are also constrained to one. This is the free-est possible model that is identified. The saturated model is used in calculating fit statistics such as the RMSEA, and Chi-squared fit indices.

The independence model, sometimes called the null model, is a model in which each variable is treated as being completely independent of every other variable. As such, all the variances and, when possible, all means are estimated. However, covariances are set to zero. Ordinal variables are handled the same for the independence and saturated models. The independence model is used, along with the saturated model, to create CFI and TLI fit indices.

The saturated and independence models could be used to create further fit indices. However, OpenMx does not recommend using GFI, AGFI, NFI (aka Bentler-Bonett), or SRMR. The page for `mxSummary`

has information about why.

When the mxFitFunctionMultigroup fit function is used, `mxRefModels`

creates the appropriate multi-group saturated and independence models. Saturated and independence models are created separately for each group. Each group has its own saturated and independence model. The multi-group saturated model is a multi-group model where each group has its own saturated model, and similarly for the independence model.

One potentially important limitation of the `mxRefModels`

function is for behavior-genetic models. If variables 'x', 'y', and 'z' are measured on twins 1 and 2 creating the modeled variables 'x1', 'y1', 'z1', 'x2', 'y2', 'z2', then this function may not create the intended saturated or independence models. In particular, the means of 'x1' and 'x2' are estimated separately. Similarly, the covariance of 'x1' with 'y1' and 'x2' with 'y2' are allowed be be distinct: *cov(x1, y1) != cov{x2, y2}*. Moreover, the cross-twin covariances are estimated: e.g. *cov(x1, y2) != 0*.

Another potential misuse of this function is for models with definition variables. If definition variables are used, the saturated and independence model may not be correct because they do not account for the definition variables.

When an MxModel has been run, some effort is made to make the reference models for only the variables used in the model. For covariance data, all variables are modeled by default. For raw data when the model has been run, only the modeled variables are used in the reference models. This matches the behavior of mxModel.

In general, it is best practice to give `mxRefModels`

a model that has already been run.

The are a few considerations specific to IFA models
(mxExpectationBA81).
The independence model preserves equality constraints among item
parameters from the original model.
The saturated model is a multinomial distribution with the proportions
equal to the proportions in your data. For example, if you have 2
dichotomous items then there are 4 possible response patterns: 00, 01, 10, 11.
A multinomial distribution for these 2 items is fully specified by 3
proportions or 3 parameters: a, b, c, *1.0-(a+b+c)*.
Hence, there is no need to optimize the saturated model.
When there is no missing data,
the deviance is immediately known as *-2 * sum(log proportions)*.
Typical Bayesian priors involve latent
factors (various densities on the pseudo-guessing lower bound, log norm
on loading, and uniqueness prior). These priors cannot be included
in the independence model because there are no latent factors.
Therefore, exercise caution when comparing the independence model
to a model that includes Bayesian priors.

Multivariate normal models with all ordinal data and no missing values
can use the saturated multinomial distribution. This is much faster than
estimation of the saturated multivariate normal model. Use
`distribution='multinomial'`

to avail this option.

`mxRefModels()`

is not compatible with
GREML expectation, as there is no
sensible general definition for a saturated GREML-type model.

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

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 | ```
require(OpenMx)
data(demoOneFactor)
manifests <- names(demoOneFactor)
latents <- c("G")
factorModel <- mxModel("OneFactor", type = "RAM",
manifestVars = manifests, latentVars = latents,
mxPath(from = latents, to=manifests, values = diag(var(demoOneFactor))*.2),
mxPath(from = manifests, arrows = 2, values = diag(var(demoOneFactor))*.8),
mxPath(from = latents, arrows = 2, free = FALSE, values = 1),
mxData(cov(demoOneFactor), type = "cov", numObs = 500)
)
factorRun <- mxRun(factorModel)
factorSat <- mxRefModels(factorRun, run=TRUE)
summary(factorRun)
summary(factorRun, refModels=factorSat)
# A raw-data example where using mxRefModels adds fit indices
## Not run:
m1 <- mxModel("OneFactor", type = "RAM",
manifestVars = manifests, latentVars = latents,
mxPath(latents, to=manifests, values = diag(var(demoOneFactor))*.2),
mxPath(manifests, arrows = 2, values = diag(var(demoOneFactor))*.8),
mxPath(latents, arrows = 2, free = FALSE, values = 1),
mxPath("one", to = latents, free = FALSE, values = 0),
mxPath("one", to = manifests, values = 0),
mxData(demoOneFactor, type = "raw")
)
m1 <- mxRun(m1)
summary(m1) # CFI, TLI, RMSEA missing
summary(m1, refModels=mxRefModels(m1, run = TRUE))
## End(Not run)
``` |

OpenMx documentation built on Oct. 5, 2018, 5:05 p.m.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.