OpenMx 2.0 Discussion

When will a version of OpenMx compatible with R3.2 be released?

> mxVersion()
OpenMx version: 2.0.0.3838
R version: R version 3.1.1 (2014-07-10)
Platform: x86_64-w64-mingw32
Default optimiser: CSOLNP

I'm using SaturatedLikelihood= and IndependenceLikelihood= commands to obtain global fit indices and noticed that the degrees of freedom associated with the chi-square statistic is incorrect. The fit of the saturated model is:

observed statistics: 2221
estimated parameters: 35
degrees of freedom: 2186
-2 log likelihood: 15732.61
number of observations: 933

and the fit of the fitted model is:

observed statistics: 2221

I am updating my metaSEM package so that it will work in OpenMx 2.0. In the Description file, it depends on OpenMx (>= 2.0.0). In the Namespace, it includes import(OpenMx).

When I check the package, it returns the following error message:
* checking whether the namespace can be loaded with stated dependencies ... WARNING
Error: object '.DollarNames' not found whilst loading namespace 'OpenMx'
Execution halted

A namespace must be able to be loaded with just the base namespace
loaded: otherwise if the namespace gets loaded by a saved object, the
session will be unable to start.

Hi,
I' getting the following errors when trying to run a submodel (i.e. putting on constraints). I have tried changing starting values, putting on bounds, but nothing seems to help. I'm running this in Windows.

Error in runHelper(model, frontendStart, intervals, silent, suppressWarnings, :
mxLog failed with errno=-1
or
Error: The job for model 'HomAceIP' exited abnormally with the error message: MxComputeGradientDescent: fitfunction HomAceIP.fitfunction evaluated to 1.#INF00 ()
or

Hi,
I am getting a warning that some dependencies are not available each time I download the beta. Is this sth to worry about? The analysis will still run.

> source('http://openmx.psyc.virginia.edu/getOpenMxBeta.R')
Installing package into ‘M:/pc/Dokumenter/R/win-library/3.1’
(as ‘lib’ is unspecified)
Warning: dependencies ‘digest’, ‘RcppEigen’ are not available
trying URL 'http://openmx.psyc.virginia.edu/testing/bin/windows/contrib/3.1/OpenMx_2.0.0-3838.zip'
Content type 'application/zip' length 3404044 bytes (3.2 Mb)
opened URL
downloaded 3.2 Mb

The following issues are mentioned in the second beta's announcement:

• IMPORTANT: Windows users should NOT try to run OpenMx with multiple threads. Doing so could crash R. Users encountering unexplained crashes under Windows may wish to set the option mxOption(NULL,"Number of Threads",1).
• The log of changes in the beta version's User Guide is not current.
  

Hi,
I'm having problems with R crashing with the new beta version. It does not seem to be a problem of excess memory usage, as I've seen in another thread here. What is strange is that after I've downloaded the beta things look fine (see mxVersion below), and it works with a simple twin multivariate saturated script. When I try a more complicated model, this is when R crashes.

> mxVersion()
OpenMx version: 2.0.0.3751
R version: R version 3.1.1 (2014-07-10)
Platform: i386-w64-mingw32
Default optimiser: CSOLNP

Here is information from Windows event Viewer when R crashes:

Should I expect confidence intervals to be working for csolnp at the moment?
I ran a number of fits, generating random start values, using the same data. The following plots of four of the model parameters and confidence intervals show quite some inconsistency / strangeness, I think. Although the likelihood varied from fit to fit, the differences were very slight. Estimates were all from good fits (no code reds / blues).

The x axis of the plot is both parameter (1-4) and within each parameter section is iteration 1-20. scaled likelihood values 1-20 are in the top section of plot.

Hi,
I am trying the Beta version to run a multivariate saturated model with ordinal data. However, I am getting no estimated parameters. Could I be missing something in the script (attached)? Any help is greatly appreciated.

Here is the output message:

> SatFit<- mxRun(SatModel, intervals=F)
Running Sat
> summary(SatFit)
Summary of Sat

observed statistics: 7099
estimated parameters: 0
degrees of freedom: 7099
-2 log likelihood: NA
number of observations: 1394
Information Criteria:
| df Penalty | Parameters Penalty | Sample-Size Adjusted

I am using the newly released beta version (2.0) to run a multivariate twin analysis and have some questions. Due to rather low prevalence (but high number of observations, so have good power anyway) I am having some numerical issues but theses seems to be handled by using the options mvnMaxPointsA/B/C and mvnAbsEps in mxOption. However, I am curious of what these actually do, and cannot find it in the documentation. I am assuming the mvnMaxPointsA/B/C adjust the number of integration points in the numerical integration of the multivariate normal density.