mxGetExpected {OpenMx} | R Documentation |
This function extracts the expected means, covariance, or thresholds from a model.
mxGetExpected(model, component, defvar.row=1, subname=model$name) imxGetExpectationComponent(model, component, defvar.row=1, subname=model$name)
model |
MxModel object from which to extract the expectation component. |
component |
A character. The name of the component to extract. |
defvar.row |
A row index. Which row to load for definition variables. |
subname |
Name of the submodel to evaluate. |
The expected means, covariance, or thresholds can be extracted from Normal (mxExpectationNormal), RAM (mxExpectationRAM), and LISREL (mxExpectationLISREL) models. When more than one component is requested, the components will be returned as a list.
If component 'vector' is requested then the non-redundent coefficients of the expected manifest distribution will be returned as a vector.
If component 'standVector' is requested then the same parameter structure as 'vector' is returned, but it is standardized. For Normal expectations the covariances are returned as correlations, the means are returned as zeros, and the thresholds are returned as z-scores. For the thresholds the z-scores are computed by using the model-implied means and variances.
Note that capitalization is ignored for the 'standVector' option, so 'standvector' is also acceptable.
See details.
The OpenMx User's guide can be found at http://openmx.ssri.psu.edu/documentation.
#---------- require(OpenMx) manifests <- paste("x", 1:5, sep="") latents <- c("G") factorModel <- 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), mxPath(from = 'one', to = manifests)) mxGetExpected(factorModel, "covariance") # oops. Starting values indicate a zero covariance matrix. # Probably should adjust them.