OpenMx General Help

Genetic correlation equals to 1 or -1 in Bivariate ACE model
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?

excessive memory usage when definition variables vary substantially

Specifying a Mediation Model in OpenMx using metaSEM
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.

twinData zygosity
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?
- Read more about twinData zygosity
- 2 comments
- Log in or register to post comments

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

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

how to deal with twin data with inexact zygosity diagnosis in OpenMx
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

Estimate individuals' values or scores on a latent variable
Thanks in advance!

how to test is a model is a RAM model
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
- Read more about how to test is a model is a RAM model
- 5 comments
- Log in or register to post comments

Source codes
- Read more about Source codes
- 2 comments
- Log in or register to post comments
Pagination
- Previous page
- Page 11
- Next page