Attachment | Size |
---|---|
ordinal_ACE_version_test.R | 7.18 KB |
Hello,
I am running an ACE analyses on ordinal variables and I am finding that my script is failing with an error about stack imbalance when I use OpenMx 2.6.9 but not when I use OpenMx 2.5.2 Below is the Error I received, and I have attached a script that recreates the error when using 2.6.9 and finishes without error when using 2.5.2 adapted from Hermine Maes example script (http://ibg.colorado.edu/cdrom2016/maes/UnivariateAnalysis/onea/oneACEoa.R)
Thanks in advance for all of your help!
Britt
Error
Calls: mxTryHardOrdinal -> mxTryHard
Error in runHelper(model, frontendStart, intervals, silent, suppressWarnings, :
Confidence intervals are not supported for units 0
Error: C stack usage 845910713468 is too close to the limit
Execution halted
Warning: stack imbalance in 'lazyLoadDBfetch', 76 then 71
Warning: stack imbalance in '<-', 30 then 29
Errors during final fit for Hessian/SEs/CIs
Warning: stack imbalance in '{', 27 then 26
Warning: stack imbalance in 'if', 25 then 24
Hi, Brittany. What operating system and version of R are you using?
I ask because I am unable to reproduce the error you report on one of my systems, which is 32-bit Windows (on Intel x86) with R version 3.1.1 (and running OpenMx v2.6.9, of course). However, your script causes a segfault(!) and crashes R under x86_64 Linux/GNU with R version 3.3.1.
Edit: Here's the text of the error message I get, under Linux:
Below is the operating system and R version that I am using.
R version 3.2.2 (2015-08-14)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 15.10
OK, I'll look into it. In the meantime, would you mind trying your script again in OpenMx v2.6.9, but using NPSOL instead of SLSQP? Use
mxOption(NULL,"Default optimizer","NPSOL")
to switch. I'm curious to know if you still get similar errors, and if the results (parameter estimates, etc.) look reasonable.Thanks so much!! I tried running the script again with the NPSOL optimizer, however I got the following error:
I then tried getting the optimizer with the function and received a fatal error. So I am not sure how to proceed to test the NPSOL optimizer. Thanks!
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
omxMLFitFunction.cpp:19:33: fatal error: stan/math/mix/mat.hpp: No such file or directory
compilation terminated.
/usr/lib/R/etc/Makeconf:143: recipe for target 'omxMLFitFunction.o' failed
make: *** [omxMLFitFunction.o] Error 1
ERROR: compilation failed for package ‘OpenMx’
* removing ‘/export/home/brittany/R/x86_64-pc-linux-gnu-library/3.2/OpenMx’
* restoring previous ‘/export/home/brittany/R/x86_64-pc-linux-gnu-library/3.2/OpenMx’
The downloaded source packages are in
‘/tmp/Rtmp5wFFFw/downloaded_packages’
Warning message:
In install.packages(pkgs = c("OpenMx"), contriburl = contrib.url(repos, :
installation of package ‘OpenMx’ had non-zero exit status
Oops forgot to specify that when I tried to get the NPSOL optimizer I used the omxGetNPSOL() function. Thanks!
You can install NPSOL-enabled OpenMx by following the instructions in our Download page. Basically, copy-paste
into your R prompt.
Thanks! And to clarify will this over write the OpenMx version 2.5.2 I have installed?
Most likely, yes, it will.
Looks like that new OpenMx with the NPSOL fixed it, thanks!! :)
That's good to hear!
I'd like to thank you for reporting the error you encountered, and especially, for including a script that reproduces it. While looking into the matter, I noticed that
mxTryHard()
calculates and reports standard errors when there are MxConstraints involved. It's not supposed to, though. So, I fixed that behavior in the source repository, and it will be fixed in the next binary release. Making that change to the source also appears to repair the issue you encountered with SLSQP. In the meantime, don't trust the standard errors you're getting if you have any MxConstraints in your MxModel; go by the confidence intervals instead.Thanks for the heads up on the MxContraints! :)
Looks like that new OpenMx with the NPSOL fixed it, thanks!! :)