OpenMx General Help

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?

Posted on
No user picture. CharlesD Joined: 04/30/2013

excessive memory usage when definition variables vary substantially

Hi. I've been doing some work on a package to overlay openmx and manage continuous time SEM models. Everything is working great but I am quite stuck with what seems to be some sort of memory management issue with openmx. I've spoken with Tim Brick regarding this previously, but thought I'd put it up here for all to see :)

Posted on
No user picture. xander Joined: 11/09/2013

Specifying a Mediation Model in OpenMx using metaSEM

Hello All,

First of all thank you for your hard work over the years to provide us with unique ways to combine meta-analytic principles with SEM--and providing a conduit to do this in R. My research team and I really appreciate it.

I am doing a two stage MASEM with the metaSEM package, and am having some difficulty specifying a mediation analysis in the A/Smatrix. I have four variables: w,x,y,z. And my linear equations are as follows:

X ~ W + e
Z ~ W + e
Y ~ W + X + Z + e

I would like to see if X or Z (or X and Z) mediate the relationship between W and Y.

Posted on
No user picture. IvanVoronin Joined: 08/18/2013

twinData zygosity

Hello!

I use Australian twin sample biometric data to test my scripts, but there is something unclear with zygosity variable: it has 10 levels instead of 5 (as declared on help page). How should I actually interpret these levels?

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

Will someone check over my script and model (latent change model)?

So I just started learning SEM this past year, teaching myself with some guidance from my advisor. I learned on AMOS, but I'm trying to teach myself OpenMX (just started last week).

Posted on
Picture of user. rabil Joined: 01/14/2010

OpenMx Refuses to Compute Confidence Intervals Using Ubuntu and Windows - But Not Mac

I adapted OpenMx ordinal threshold model example code to estimate a model with 5 raters using a 5-point ordinal scale. (The code uses threshold deviations.) The model converges (after some tweaking of starting values) with no errors and no warning messages - this is on a computer running Ubuntu 13.04 and the latest version of OpenMx. The same model was estimated using MPlus - both yield very similar factor loading estimates and threshold estimates (although MPlus used weighted least squares - ML wouldn't work).
Posted on
No user picture. shaoyong.su Joined: 04/23/2013

how to deal with twin data with inexact zygosity diagnosis in OpenMx

Dear OpenMx experts,

I'm now working on a twin dataset with zygosity diagnosis is not exact, i.e. I only have one variable indicating the probability of one twin pair being MZ or DZ. In Mx we can "specify" this variable and "weight" it when estimating the covariance matrix. Dr. Neale has discussed this issue and provided the Mx script file here:
ftp://ftp.vcu.edu/pub/mx/examples/zygosity

The key commands are below:

G2: All twin pairs
Data NInput_vars=3 Nmodel=2
Labels t1 t2 pmz
Rectangular File=sim1.mz
Select t1 t2 pmz ;
Definition pmz ;

Matrices= Group 1
P Full 1 1

Posted on
No user picture. dadrivr Joined: 01/19/2010

Estimate individuals' values or scores on a latent variable

Is it possible with OpenMx to estimate individuals' values or scores on a latent variable? For example, using the one factor model (http://openmx.psyc.virginia.edu/svn/trunk/demo/OneFactorModel_PathRaw.R), is it possible to estimate what each case's value is on the latent factor? If so, how can this be done?

Thanks in advance!

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

how to test is a model is a RAM model

The new 2-part objectives in OpenMx mean an older method for testing if a model is a RAM model don't work

I used to say this:

(isS4(obj) && is(obj, "MxModel") && class(obj$objective)[1] == "MxRAMObjective")

What code do people recommend for backward (and preferably forward compatible testing of model type? Do we need to add a type= field to MXModel?

best, tim

Posted on
No user picture. Joper Joined: 02/23/2013

Source codes

Hello! Does anybody know how can we download the source codes of OpenMx? We didnt find links to the source codes from the OpenMx Website. Thank you!