Attachment | Size |
---|---|
Multi_ACE_Twin Nov 2016 By Gender Nov 10.R [6] | 16.6 KB |
Table 1.pdf [7] | 43.16 KB |
I am running a multivariate model (3 dependent vars) with sex (female MZ, male MZ, female DZ, male DZ, F&M DZ, M&F DZ) the model runs fine with the SLSQP optimizer.
No SEs are suspect and the return code is zero.
However when I try the command:
FitCholMod <- mxRun(CholMod, intervals =TRUE)
It runs as before but I receive this error message.
Error in runHelper(model, frontendStart, intervals, silent, suppressWarnings, :
NLOPT fatal error -1
So the confidence intervals are not produced.
I have tried the other optimizes. With CSOLNP a number of SEs are suspect. With NPSOL MX returns a code of 1 and more SEs are suspect.
I can obtain CIs with a worse fitting (non-preferred) model with the SLSQP optimizer. The non-preferred model has a significantly worse fit than the full-Cholesky, so I don't want to use that model.
With other models, just MZ vs DZ and a cohort model (young MZ, old MZ, young DZ, old DZ) I can obtain CIs.
I have several suspicions:
- I don't have enough male dizygotic twins !!!!
- I should rescale the dependent variables (years of education is left as is 8-18, income have been divided by $10,000, but occupational status still ranges from 0-100). I could divide occupational status by 10.
*It did seem to work before I downloaded the new MX version but I am at all not certain that I was analyzing the same model.
Any suggestions?