Revision of OpenMx for mx 1.x users from Mon, 08/09/2010 - 05:46

The revisions let you track differences between multiple versions of a post.

Wiki home page

Please add material here as you learn...

If you have questions not answers, then add those here: That's how a wiki works.

There are a number of changes between Mx and OpenMx, mostly easy to learn.

You will still need to read in data, propose algebra, possibly create more than one group, estimate the model, and post-process and display the results. Much of the structure of an mx script translates easily into OpenMx, just with more structure. COnstraints becomes mxConstraint(), data statements such as SElect are replaced with calls to mxData(), and DEfinition-variables are implemented by . COvariance and MEans statements are replaced with calls to mxAlgebra(), mxMLObjective(), mxRAMObjective(), or mxFIMLObjective.. etc.

Data

R has much more powerful access to data than did mx, including reading a range of formats, not barfing on names, and accessing data from a range of locations, including over the web etc.
You will need to replace Mx1.x Data, Missing, file, and SELECT statements with some R code to read in the data and a call of mxData()

This, for instance:

 Data NI= 13 ! nr of inputvars per family
 Missing=-1.000 ! missing values must be exactly this string "-1.000"
 RE file=example.dat ! read in raw datafile, must be in same directory and not too long a path

would become:

	<em># Read in data (could be a URL. Names are in the file (Header =T))</em>
	data &lt;- read.table("example.dat", header=TRUE, sep="\t")
	<strong>mxData</strong>(data, type="raw") # apply the data to a model

DEfinition variables (covariates)

In mx 1.x, to model a variable at the individual level, you DEfined these variables from among the SElected variables.

In OpenMx, you no longer need to explicitly DEfine variables. Instead prefixing the name of a path or matrix label with "data." causes OpenMx to evaluate the matrix fit on subject-level data.

OpenMx automatically removes any such data.variables from the variables being models, and reserves it as a definition variable.

PS: this explains why "." style names (though legal in R) are illegal in OpenMx.

From Groups to Models

Mx used groups to structure its analyses. The comparable structure in OpenMx is the model. Because R provides a lot of functionality, the mapping is not complete. This is clearest in mx groups used for standardising results rather than defining a SEM model. These will often be replaced entirely with manipulation of the data in regular R code. In multi-group models, mx groups will often be mapped directly into OpenMx (sub-)models. Thus a two-group twin design may well be implemented with two models in OpenMx, contained in a parent model. mx scripts often begin with a "Define" group specifying a set of matrices, and this can usually be translated directly to an OpenMx model.

Algebra

Operators

OpenMx uses R-style operators, which entails some very straightforward find-and-replace type changes.

Two very common functions: matrix multiplication and transposition, signified by "*" and "~" in old Mx, need to be changed to %*% and t(x) respectively.

A full list of function mappings is provided here operators and functions.

The other major change in OpenMx algebra is the addressing of symbols. For instance, in Mx 1.x, to include a matrix from another group in an algebra, you would import the matrix into the second group, then refer to it by its name. In OpenMx, you simply refer to objects in other models by their full name, i.e., to refer to matrix "A" in model "mz", you simply say "mz.A".

So if in mx you said
Begin Matrices = Group 1
Means M ;

In OpenMx, you would say
"mxFIMLObjective("expMean", "name.M"),

where "name" is the name of what was group 1.

To compare fits between models use mxCompare().

To set the saturated 2*LL, you say:

	 <strong>summary</strong>(modelOut, SaturatedLikelihood=300)

Things that are no longer needed

save and get

The result of an MxRun is itself a model, populated with the final estimated values of the matrices and algebras you proposed.

So you can say:

  fit&lt;- mxRun(model)
  fit@matrices$a@values=0
  fit@matrices$a@free=FALSE
  fit2&lt;- mxRun(fit) 

  • Option Sat

    Options like NDecimals have not been re-implemented in OpenMx, as these are handled flexibly in R using round() etc.

    \part(A,B)

    Where you would previously have called \part(A,B) to extract the rectangle described by B out of matrix A, now you can use algebraic matrix labels and Rs matrix addressing to meet these needs.

    Confidence intervals

    Numerical estimates of Standard Errors are part of the standard Summary().

    Maximum likelihood confidence intervals are requested by including mxCI() in your model, then mxRun(model, intervals=T) to compute the CI's (which are then shown in the summary() )

  • Comments

    deadpixel's picture

    Great insight