OpenMx Help

Posted on
No user picture. Veronica_echo Joined: 02/23/2018

bootstrap coverage probability and the bootstrap replication

Hi everyone,

I am conducting a simulation study of a growth model and would like to evaluate the bootstrap CP of it. I kept the simulation replication as 1000 and set bootstrap replication as 1000 and 2000, respectively. The results seemed wiered, since the CPs of bootstrap 1000 (all of them were located between (0.93, 0.97)) were much better than those of bootstrap 2000 (some CPs were quite low, say 0.86). Any advice about this issue? Should I increase the bootstrap replication to a larger number, say 5000? Thank you in advance.

Posted on
Picture of user. dtofighi Joined: 10/12/2009

OpenMx on R 3.5.0 (2018-04-23)

Hello,

I used the following code to install the latest version of OpenMx on macOS 10.13.4. I got the following error.

Posted on
No user picture. felix.fischer Joined: 04/26/2018

compilation of OpenMx fails in ubuntu 14.04

Hi,

i am trying to install OpenMx on a ubuntu server, R version 3.4.4, 64 bit, but compilation fails and i cannot really make sense of the error message. I ran sudo apt-get upgrade and dist-upgrade, but it did not help.


ffischer@common-metrics:~$ sudo su -c "R -e \"install.packages('OpenMx', repos = 'http://cran.rstudio.com' )\""
[sudo] password for ffischer:

Posted on
No user picture. Veronica_echo Joined: 02/23/2018

parallel computation in linux cluster

Hi everyone,

I want to conduct several 10,000-simulations for an OpenMx model in Linux cluster. I have over 20 parameters in the model and for each parameter, I need 1000-bootstrap confidence interval to test the coverage. I've add lines in below screenshot, and there are 16 cores available for OpenMx. Even though, the simulation may take about one month at my best guess. Could someone kindly advise some other ways to speed up my simulations? Thanks in advance.

Posted on
No user picture. Veronica_echo Joined: 02/23/2018

Create a function to generate mxMatrix and mxAlgebra

Hi everyone,

I am writing an R function which can create and run an OpenMx model. I can do it in a simple way. Partial codes are shown below. However, I have over 20 definition variables and for each I need corresponding mxAlgebra. I may want to write a function/a loop to generate mxMatrix and mxAlgebra. I've tried a couple of methods, yet none of them worked. The mainly problem is that in the expression of mxAlgebra, I need call the name of corresponding mxMatrix. Thank you in advance!

Posted on
No user picture. amr438 Joined: 03/09/2016

Converting ClassicMx script to OpenMx

Hello,

I am currently trying to recreate ClassicMx script for an ECOT model into OpenMx and I am coming across an error that I am not sure how to solve.
When I run the model I get the following error:

Error: The following error occurred while evaluating the subexpression 'solve(ecot.matF %*% solve(ecot.matI - ecot.matA) %*% ecot.matM %*% ecot.matU)' during the evaluation of 'MZtoss.expMeanMZtoss' in model 'modelECOT' : 'a' (4 x 1) must be square

I have attached the ClassicMX code and OpenMx code.

Posted on
No user picture. jkarch Joined: 03/15/2011

I fail to install OpenMx from source on MacOs 10.12.6

The short story is that R seems to ignore my intention to use the gcc* compilers.

This is what I do (following https://openmx.ssri.psu.edu/wiki/howto-build-openmx-source-repository):

Posted on
No user picture. jkarch Joined: 03/15/2011

Performance issues many algebras

Hey,

I am currently implementing a new model specification language for longitudinal panel models. As "backend", I use OpenMx. Everything works but it is painfully slow.

Posted on
Picture of user. paulT Joined: 11/03/2017

Empirical Underidentification with a bifactor type model

First post by new user: Can anyone give some advice on an empirical under-identification issue please? I am trying to fit the following model:

resVars <- mxPath( from=mylabels, arrows=2,
free=TRUE, values=rep(1,12),
labelatVars <- mxPath( from=c("X1","X2"), arrows=2, connect="unique.pairs",
free=c(TRUE,FALSE,TRUE), values=c(1,0,1), labels=c("varX1","cov","varX2") )