Attachment | Size |
---|---|
ACE script (KB).txt [6] | 12.75 KB |
Univariate ACE model BMI output.txt [7] | 38.17 KB |
Hi there,
I'm very new to OpenMx and working with twin data (OpenMx version: 2.5.2 [GIT v2.5.2] R version: R version 3.2.4 Revised (2016-03-16 r70336) Platform: x86_64-w64-mingw32 Default optimiser: SLSQP).
I created a script based on online resources from this site, where I want to run saturated, ACE and submodels, controlling for age and gender.
As far as I can tell, the script is running, with no warnings/ errors. However, I am getting p values of 0 for the Chi square difference test for the ACE model compared to saturated. I am wondering if this is alright/ presentable?
When comparing submodels for the best fit, do we compare to saturated or to the ACE model? when comparing nested submodels to ACE, I am getting values of 1 and 0.
I am working with small numbers (45 MZ and 24 DZ pairs), and raw data (not normalised), would this be an issue?
I would just like someone to check over the script, so I can confidently present my results!
I've attached a script and an output for my BMI data.
Appreciate any help/ tips you can give! :)
Rebecca