ACE model error: In model 'ACE' Optimizer returned a non-zero status code 6.
Posted on

Forums
Hello,
I am quite new to behavioral genetic modelling and I have a question about an error that I am receiving since the OpenMX update.
I have used the same scripts as before the update, but my results are very different.
in addition, I now get the error:
In model 'ACE' Optimizer returned a non-zero status code 6. The model does not satisfy the first-order optimality conditions to the required accuracy, and no improved point for the merit function could be found during the final linesearch (Mx status RED)
Does anyone else received this message since the update?
Is there a way to solve it?
My data and scripts are attached.
The "old" results were: A=0,55 C=0 E=0,45, AIC= 252,12
Many thanks in advance!
Best, Michelle
switch optimizers?
FTR, the message you're reporting is a warning, not an error. Anyhow, you might be able to avoid it if you switch optimizers. Put
mxOption(NULL,"Default optimizer","SLSQP")
in your script, after you load OpenMx.
Log in or register to post comments
In reply to switch optimizers? by AdminRobK
SLSQP works for me
I tried running your script with SLSQP instead of CSOLNP (which is now the on-load default). I recovered the standardized point estimates and AIC you report in your post.
Log in or register to post comments
In reply to switch optimizers? by AdminRobK
Thank you!
I added the mxOption to the script, and now it works fine!
Thank you so much :)
Best, Michelle
Log in or register to post comments
starting values
Mike Neale pointed out off line that this is a starting value problem, where some of our optimisers can cope beginning a very long distance from the solution
Starting a, c, and e at 1/3 * sd is best. In multivariate space you'd start with only the diagonals of Cholesky matrices initialized at these values (and zero elsewhere).
#Figure out reasonable starting values
svMean <- colMeans(AmygdalaHPC[AmygdalaHPC$zygosity==1, 2:3])[1]
svPaths <- sqrt(var(AmygdalaHPC[AmygdalaHPC$zygosity==1, 2:3], na.rm=T)[1,1])/3
Log in or register to post comments