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

2 replies [Last post]
shaoyong.su's picture
Offline
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:
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
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.

Best,
Shaoyong

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

Thanks so much. This is realy helpful!