OpenMx Help

Posted on
Picture of user. siti nur azizah Joined: 03/12/2014

modified model

hi, I am a beginner using Mx open packages, so I do not quite understand. I have run my program. output is given as attached.
The value of RSMEA, CFI and TLI its very ugly because its value doesn't meet either criteria . What can I do in advance to improve the model?
thank you.

Posted on
No user picture. mverdam Joined: 02/04/2014

Referring to an OpenMx matrix in the OpenMx environment

I'm trying to run a model in which I have several matrices that should contain the same parameter estimates (they overlap partly). I know I can use labels or mxConstraint to put equality constraints on parameter estimates, but I was wondering if it is also possible to refer directly to the parameter estimates of a matrix in other parts of the OpenMx environment. For example, I have tried to build an mxAlgebra that refers to the parameter estimates of an mxMatrix by using:

mxAlgebra(expression=cbind(Model.mxMatrix[,i],Model.mxMatrix[,j]), name="part_mxMatrix")

Posted on
No user picture. mverdam Joined: 02/04/2014

Avoid default abort when error occurs

I'm trying to run multiple SEM-models within a loop (i.e., for simulation purposes), but am experiencing some problems when a model does not run properly (i.e., objective function returns NaN). The default action is to abort and display an error message, which results in aborting the whole loop.

I am looking for a way to make it possible to make sure the loop will continue when one of the models is not returning a result (and of course, I would also like to save the non-result in some way).

Posted on
No user picture. fife Joined: 07/01/2010

Computing R squared with missing data

This may be an easy question, but I can't think of the answer. I've built a regression model (predicting Y) with two independent variables (X and Z). I want to compute R^2. Is that built into the model somehow? If not, any ideas on how to compute it? Here's the model I have:

multiRegModel <- mxModel("Multiple Regression, All Variables",
type="RAM",
manifestVars=c("x", "y", "z"),
# variance paths
mxPath(
from=c("x", "y", "z"),
arrows=2,
free=TRUE,
values = c(.5, .5, .5),
labels=c("varx", "residual", "varz")
),

Posted on
Picture of user. Julia Joined: 03/29/2012

Adding new algebras in submodels and constraining them

Hi.
I am running a bivariate moderation model and would like to set some constraints in order to test for some nonlinear effects. I am struggling though with adding new algebras into submodels (with subsequent equating them). Let's say that I specify


pathAm <- mxMatrix(name = "am", type = "Lower", nrow = nv, ncol = nv, free=T, labels = c("amM","amC","amU"), values=pathVal)
pathCm <- mxMatrix(name = "cm", type = "Lower", nrow = nv, ncol = nv, free=T, labels = c("cmM","cmC","cmU"), values=pathVal)

Posted on
No user picture. a9mike Joined: 06/28/2013

Confidence Intervals

I'm not sure the best way to get confidence intervals for my estimates. I'm doing a bifactor model with a large dataset (HRS) with lots of missing data, and the mxCI are taking impossibly long (days long).

Anyone have suggestions? Would bootstrapping be faster? If so, what would that script look like? (I've never bootstrapped before)

Posted on
Picture of user. neale Joined: 07/31/2009

International Workshop on Statistical Genetic Methods for Complex Human Traits. March 3-7, Boulder Colorado

This workshop will be taught at the Millennium Harvest Hotel in Boulder Colorado, March 3-7 2014, and will use OpenMx extensively. Those interested in a general education in OpenMx use would be well-served by the course, at the expense of learning about genetics and modeling data from twins and other relatives.

http://www.colorado.edu/ibg/content/march-3-march-7-2014-2014-international-workshop-statistical-genetic-methods-human-complex

Posted on
No user picture. cuttance Joined: 01/20/2014

OpenMx in the Cloud

I’m looking to integrate MSEM (one of: OpenMx, MPlus, LISREL, etc) into an application I am developing for schools on a .NET Microsoft Azure Cloud platform — multiple schools will access the application simultaneously, and the application will call the MSEM modelling unit to estimate a specific model and return the estimates to a database for benchmarking purposes.

Would I need to run OpenMx inside a VM, or is there a server version? Can multiple instances of OpenMx run inside the same VM? or Would we need to run multiple VMs, each for an instant of OpenMx?

Posted on
No user picture. c3094213 Joined: 12/30/2013

Genetic correlation equals to 1 or -1 in Bivariate ACE model

Hello All,

I had just done bivariate ACE model for two continuous variables (bivariate correlated factors model). The results turned out to be a little weird as follows:

A1 = 0.25, A2 = 0.09
C1 = 0.12, C2 = 0.08
E1 = 0.62, E2 = 0.82
Rph = -0.24
Ra = -1
Rc = 1
Re = -0.26

We observe Ra is -1, but we personally think it is not reasonable that the two variables are perfectly correlated. Is it possible that there are some problems with the program?