Tue, 06/17/2014 - 09:46

Hi,

I am new to OpenMx and trying to use it for twin study (monozygotic vs dizygotic twins). I have some questions and hope someone may help me with the following questions.

I tried the following two scripts:

http://openmx.psyc.virginia.edu/svn/tags/stable-1.2/demo/UnivariateTwinA...

http://openmx.psyc.virginia.edu/svn/tags/stable-1.2/demo/UnivariateTwinA...

I knew their input data were different, but is there any difference for data processing steps between these two?

I see in order to decide between ACE and AE models, the aforementioned scripts do the log-likelihood ratio test at the end of the script. They got the difference between the log-likelihood of AE model and ACE model. But how could I decide which model I should choose based on this value? Should I do chi-square test using this value and degrees of freedom?

Also, is there any script available online which could also calculate the p-value for the parameters a2, c2 and e2?

Thank you very much!

These scripts appear to me to be doing the same model.

You could compare the AE model with the ACE model by using the

mxCompare()function which does a likelihood ratio (change Chi-Squared / change df) test and shows AIC.To evaluate the significance of parameters in any model, you could include confidence intervals wtih

mxCI()objects in the model that refer to any named matrix, algebra or parameter in the model. You could also construct models where parameters were set to zero and usemxCompare()as above to do likelihood ratio tests to see if estimating that parameter is worthwhile.