You are here

Home › Forums › OpenMx 2.0 Discussion

The OpenMx website will be down for maintenance from 9 AM EDT on Tuesday, September 17th, and is expected to return by the end of the day on Wednesday, September 18th. During this period, the backend will be updated and the website will get a refreshed look.

R crashes

10 posts / 0 new
Log in or register to post comments
Last post
Thu, 09/04/2014 - 05:46
#1
karobro's picture
karobro
Offline
Joined: 01/16/2013 - 10:41
R crashes

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:

Faulting application name: Rgui.exe, version: 3.11.579.0, time stamp: 0x53be708a
Faulting module name: OpenMx.dll, version: 0.0.0.0, time stamp: 0x53ee0f11
Exception code: 0xc0000005
Fault offset: 0x0001a2be
Faulting process id: 0x1820
Faulting application start time: 0x01cfc778c91b0525
Faulting application path: C:\Program Files\R\R-3.1.1\bin\i386\Rgui.exe
Faulting module path: C:\Program Files\R\R-3.1.1\library\OpenMx\libs\i386\OpenMx.dll
Report Id: 3be81308-336c-11e4-87e3-005056ad3f12

Does anyone have any idea what could be causing this? I've tried different computers and R Versions.

Top
Thu, 09/04/2014 - 18:13
RobK's picture
RobK
Offline
Joined: 04/19/2011 - 21:00
NPSOL?

Does it crash when NPSOL is the default optimizer? You can switch from CSOLNP to NPSOL with mxOption(NULL,"Default optimizer","NPSOL").

Also, out of curiosity, what do you see when you enter options()$mxOption$"Number of Threads" into R?

Edit:
Two other things. First, on which other systems have you observed this crash? Second, could you attach your script or one similar to it?

Top
Fri, 09/05/2014 - 09:33
(Reply to #2)
karobro's picture
karobro
Offline
Joined: 01/16/2013 - 10:41
mxOption$Number of threads = 0

Thanks for the reply!

Yes, it also crashes with NPSOL.

This is what comes out when entering
options()$mxOption$"Number of Threads"
[1] 0

It also crashed on R 3.0.1, and I've only tried Windows. I should ask the author of the script for permission to post it, so I will come back to this.

Top
Fri, 09/05/2014 - 09:52
(Reply to #3)
tbates's picture
tbates
Offline
Joined: 07/31/2009 - 14:25
try 1 core?

try running this first to stop it trying to go multi-threaded?

mxOption(NULL, "Number of Threads", 1)
Top
Mon, 09/08/2014 - 05:00
(Reply to #4)
karobro's picture
karobro
Offline
Joined: 01/16/2013 - 10:41
still crashing

Thanks, but this unfortunately did not work either.

Attached is a version of the script I've used.

Top
Mon, 09/08/2014 - 11:47
(Reply to #5)
RobK's picture
RobK
Offline
Joined: 04/19/2011 - 21:00
What do nv, nth, and Vars

What do nv, nth, and Vars happen to be in your case?

Top
Mon, 09/08/2014 - 13:58
(Reply to #6)
karobro's picture
karobro
Offline
Joined: 01/16/2013 - 10:41
nv, nth, Vars

nv <- 3
nth <- 6
Vars <- c('M07nnT','F07nnT','T07nnT')

Top
Mon, 09/08/2014 - 16:37
(Reply to #7)
RobK's picture
RobK
Offline
Joined: 04/19/2011 - 21:00
Bug with mxFitFunctionMultigroup() ?

I think there's a bug with mxFitFunctionMultigroup(). Try redefining your function to

makeModel  <- function(name) {
  parsZf    <- list( pathAf, pathCf, pathEf, pathRaf, pathRcf, pathRef, covAf, covCf, covEf, covPf, estVarsZf )
  parsZm    <- list( pathAm, pathCm, pathEm, pathRam, pathRcm, pathRem, covAm, covCm, covEm, covPm, estVarsZm )
  parsZfm   <- list( corRao, corRco, covAmf, covCmf, covAfm, covCfm)
  modelMZf  <- mxModel( parsZf, mtGf, covMZf, dataMZf, expFunctionMZf, fitfunction, name="MZf" )
  modelDZf  <- mxModel( parsZf, mtGf, covDZf, dataDZf, expFunctionDZf, fitfunction, name="DZf" )
  modelMZm  <- mxModel( parsZm, mtGm, covMZm, dataMZm, expFunctionMZm, fitfunction, name="MZm" )
  modelDZm  <- mxModel( parsZm, mtGm, covDZm, dataDZm, expFunctionDZm, fitfunction, name="DZm" )
  modelDZo  <- mxModel( parsZf, parsZm, parsZfm, mtGo, covDZo, dataDZo, expFunctionDZo, fitfunction, name="DZo" )
 
  mgalg <- mxAlgebra(MZf.fitfunction + DZf.fitfunction + MZm.fitfunction + DZm.fitfunction, name="mgalg")
  mgfitfunction <- mxFitFunctionAlgebra("mgalg")
 
 
  name      <- mxModel( name, parsZf, parsZm, parsZfm, modelMZf, modelDZf, modelMZm, modelDZm, modelDZo, mgalg,
                        mgfitfunction)
}

and see if that prevents the crash.

Top
Tue, 09/09/2014 - 06:05
(Reply to #8)
karobro's picture
karobro
Offline
Joined: 01/16/2013 - 10:41
it worked

Yes, this prevented the crash! Thank you so much for the help!

Top
Tue, 09/09/2014 - 12:05
neale's picture
neale
Offline
Joined: 07/31/2009 - 15:14
The multi group fit function bug has been fixed, yay!

At present the fix is only available for those building from source. Beta 3 should include the fix when it is released. Sorry for the inconvenience.

Top

Copyright © 2007-2024 The OpenMx Project