OpenMx 2.0 Discussion
carey
Joined: 10/19/2009
NLOPT fatal error -1
Am helping a grad student who got the error message:
Error in runHelper(model, frontendStart, intervals, silent, suppressWarnings, :
NLOPT fatal error -1
When I search on "NLOPT fatal error -1" on the site, the only thing I get is:
Rev 4444 and /trunk/ Rev 4445
which just gives code for writing ""NLOPT fatal error -1."
Anyway you can inform us mortals about what "NLOPT fatal error -1" means?
PS It would also be nice to give English outputs for "fatal error -2", -3, -4, etc.
- Read more about NLOPT fatal error -1
- 5 comments
- Log in or register to post comments
AdminRobK
Joined: 01/24/2014
Optimizer status code 6 (Mx status RED) in OpenMx v2.2.x
Some OpenMx users may have noticed that version 2.2.x of the package throws warnings about code 6 ("Mx status RED") more frequently than did prior versions. In particular, users may have noticed that OpenMx v2.2.x throws code-6 warnings in instances where earlier versions did NOT. A warning of code 6 indicates that the gradient (i.e., the vector of first partial derivatives of the fitfunction with respect to the free parameters) is not sufficiently close to the origin for the software to "be sure" it has found a minimum of the fitfunction.
rabil
Joined: 01/14/2010
SLSQP Confidence Intervals?
I ran a very simple one-factor with 3 indicators measurement error model. I used the default optimizer (SLSQP) and then the NPSOL optimizer that I've been using for a long time. Either NPSOL has been outputting bad intervals for a long time or SLSQP is not working very well for confidence intervals. Please see the attached text file that compares SLSQP, NPSOL, and lavaan results for the same model and same data. SLSQP does not compute some bounds, even though NPSOL has no problem computing bounds.
- Read more about SLSQP Confidence Intervals?
- 4 comments
- Log in or register to post comments
Cindy.s
Joined: 07/04/2015
mxRefModels output
Hi OpenMx Community!
I just ran a single-trait ACE model, with one observed variable, and no definition variables. I used mxRefModels to get the Saturated and the Independence models. I wanted to see if I get similar values defining the Saturated model and the Independence model by hand.
My results differ basically in every aspect. Is my Saturated and Independence model definition wrong (see attached script), or is this a limitation of mxRefModels:
-Df is 6 vs 3: I guess mxRefModels includes A,C,E parameters, that is why it is 3 more. Should they be included in the Saturated model?
- Read more about mxRefModels output
- 4 comments
- Log in or register to post comments
bwiernik
Joined: 01/30/2014
Information Criteria Bug?
I'm fitting a pseudo-bifactor model to a covariance matrix. When I run summary(), I get the following fit statistics:
observed statistics: 45
estimated parameters: 22
degrees of freedom: 23
fit value ( -2lnL units ): -3093.63
saturated fit value ( -2lnL units ): -3226.707
number of observations: 1131
chi-square: X2 ( df=23 ) = 133.0765, p = 1.748705e-17
Information Criteria:
| df Penalty | Parameters Penalty | Sample-Size Adjusted
AIC: 87.07654 177.0765 NA
- Read more about Information Criteria Bug?
- 2 comments
- Log in or register to post comments
joz
Joined: 06/30/2015
Scripts for Power analysis of Bivariate ACE in Twin Sample?
Hello,
This is my first time using this forum...any help is greatly appreciated!
Are there any OpenMx scripts for determining power in a twin sample using bivariate ACE analyses?
In a twin sample, I'm conducting bivariate ACE analysis of a personality trait and depressive symptoms. I found 4.6% in shared additive genetic variance, which is smaller than expected (compared to previous reports in the literature)
My univariate ACE results: A^2 for personality trait = 62%; A^2 for depressive symptoms = 32%
Salutfate
Joined: 11/22/2013
Bivariate Cholesky model adjusted for age-Error message for ACE
I am using a bivariate cholesky model for combined continuous and bivariate variables with Age effects on thresholds.
I am stuck with the error message in running model 'CholAce':A definition variable has been declared in model 'ACE' that does not contain a data set.
Could you tell me what the error means? Any help appreciated! Thanks.
Samaneh2013
Joined: 06/19/2013
Nested models in Multivariate Cholesky script - Error message for AE
Hi,
I have been running longitudinal multivariate cholesky decomposition. This seems o run ok but the problem is comparing nested models to ACE model and the Saturated model.
The error message is:
Error: The following error occurred while evaluating the subexpression 'solve(sqrt(CholAE.I * CholAE.C))' during the evaluation of 'iSDc' in model 'CholAE' : Lapack routine dgesv: system is exactly singular: U[1,1] = 0
The code:
# -----------------------------------------------------------------------------
# Cholesky Decomposition AE Model
Mike Cheung
Joined: 10/08/2009
OpenMx 2.2.3 with SLSQP optimizer
I am testing the SLSQP optimizer in the GitHub version (OpenMx 2.2.3). It seems that it will return an error code 6 when the estimates hit the lower bounds. The NPSOL optimizer seems to work fine in this case. The -2logLikelihoods of both solutions are equal (up to 4 decimal places). Please see the attached examples.
- Read more about OpenMx 2.2.3 with SLSQP optimizer
- 2 comments
- Log in or register to post comments
bee
Joined: 05/19/2015
Installation error
Hi, I installed OpenMX following instructions here http://openmx.psyc.virginia.edu/2015/01/version-201-now-available
Installation went through but I noticed an error
mkdir -p ../inst/debug
cp *.gcno ../inst/debug/
cp: cannot stat ‘*.gcno’: No such file or directory
Makevars:40: recipe for target 'save-gcno' failed
make: [save-gcno] Error 1 (ignored)
Not sure if that is important. I am using R 3.2 on Ubuntu 15.04 64 bit.
Thanks
- Read more about Installation error
- 1 comment
- Log in or register to post comments