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

2 replies [Last post]'s picture
Joined: 04/23/2013

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:

The key commands are below:

G2: All twin pairs
Data NInput_vars=3 Nmodel=2
Labels t1 t2 pmz
Select t1 t2 pmz ;
Definition pmz ;

Matrices= Group 1
P Full 1 1
End Matrices;

Specify P pmz

Means N;
Covariances A+C+D+E | A+C+D _
A+C+D | A+C+D+E _
A+C+D+E | H@A+C+Q@D _
H@A+C+Q@D | A+C+D+E ;
Weights P_I-P;

My question is how we can do the same analysis in OpenMx? Any help will be greatly appreciated.


neale's picture
Joined: 07/31/2009
Yes, in the example code's picture
Joined: 04/23/2013
Thanks so much. This is realy

Thanks so much. This is realy helpful!