I downloaded from the OpenMX website the correlated factors sex limitation model, and amended it for our data (removing the opposite sex components as our twins are all same sex).
The script runs and converges but there is no output (seeming that it is not fitting). There are no warning messages and no error messages. The output we get is:
observed statistics: 4500
estimated parameters: 0
degrees of freedom: 4500
-2 log likelihood: NA
saturated -2 log likelihood: NA
number of observations: 750
AIC (Mx): NA
BIC (Mx): NA
timestamp: 2013-11-15 09:54:40
frontend time: 9.456 secs
backend time: 0.125 secs
independent submodels time: 0 secs
wall clock time: 9.581 secs
cpu time: 9.581 secs
openmx version number: 1.0.7-1706
A number of us have looked at this, and just can't work out where the glitch is.
The data and script are attached.
Thank you so much for any assistance in advance!
|Sex limitation correlated factors.R ||19.81 KB|
|DASS_Data_Dummy.csv ||47.9 KB|