OpenMx Help

Posted on
No user picture. twoertwein Joined: 12/11/2019

lavaan multi-level SEM in OpenMX

I'm trying to replicate this example two-level SEM from lavaan in OpenMx http://lavaan.ugent.be/tutorial/multilevel.html

Lavaan's log-likelihood is -23309.87 but with the following OpenMx code I get only -26495.56. I have also tried to use the estimated parameters from lavaan as fixed parameters in the OpenMx model - the log-likelihood gets even worse then. I assume that the model structure in OpenMx is not the same as the structure model in lavaan.

Posted on
No user picture. khusmann Joined: 08/11/2019

trouble installing in conda env

Hi all,

I'm trying to get openmx to compile in a conda env (using getopenmx.R), and am getting the following compile error.

It appears it's because the conda version of zstd is later than openmx is expecting, and is using deprecated features. I poked around a little, and it looks like openmx might bundle its own copy of zstd, which it should be compiling against instead, but is instead trying to use the conda one. Any help on fixing this would be appreciated! (error messages below)

Thanks,
Kyle

Posted on
No user picture. kaje Joined: 06/05/2019

factor scores curiosity -- comparing OpenMx and lavaan

Hi,
when estimating a very simple latent variable autoregressive panel (AR(1)) model with an integrated trend component, I noticed that the factor scores based on the regression method do not include the trend (or maybe they are a kind of detrended). I'm wondering whether there maybe is a reason or a rationale behind this that I don't see. If you add the model implied means to the factor scores, you get the same factor scores as when using lavaan.

Here is a minimal working example:

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

OpenMx does not work after updating R 3.6

Hi everyone,

After updating R3.6, OpenMx doesn't work. My OS is macOS Mojave version10.14.5. I am pasting the error message I got. Any advice about it? Thanks in advance.

Veronica

Posted on
No user picture. Kenji_UIUC Joined: 05/03/2019

NPSOL Optimizer error

I'm experiencing some issues with the twin code from the 2018 IBG workshop. It seems like the NPSOL optimizer isn't working.

Code: https://ibg.colorado.edu/cdrom2018/deleeuw/sexLimitation/oneACE5cSOLVED.R

source: https://ibg.colorado.edu/cdrom2018/deleeuw/sexLimitation/miFunctions.R

data: https://ibg.colorado.edu/cdrom2018/deleeuw/sexLimitation/sexLimACE.dat

The following happens :
###what I type
source("miFunctions.R")

Posted on
Picture of user. Mike Cheung Joined: 10/08/2009

CIs on mxAlgebra

Hi, all.

I have a mxAlgebra of a parameter multiplied by a constant, say new_x = 2*x. When I construct an LBCI on both x and new_x, I expect that the CIs on new_x are the same as the CIs on x multiplied by 2. It turns out that they are not exactly the same (see the "diff" in the following output). When the constants are larger, say 5 or 10, the CIs on the new_x even become NA.

Any ideas why this happens? Thanks.


## Multiplied by 2

Posted on
Picture of user. Mike Cheung Joined: 10/08/2009

CFA with categorical variables

Hi all,

I am fitting some CFA models with all categorical variables with FIML. The information matrix is often not positive definite when the number of variables is larger than 4 or above. Attached are some examples.

Any suggestions to improve its stability? Thanks in advance.

Best,
Mike

Posted on
No user picture. Aivil Joined: 03/14/2019

Error Message - LCGA 1 Class

Hello,
Could you help me understanding the following error message when fitting a 1class LCGA (Script and summary attached)?
Eventually, I obtain the results, but I am a bit skeptical due to the error message.

Fit attempt 10, fit=364650.77991012, new current best! (was 366100.551846931)Error in array(x, c(length(x), 1L), if (!is.null(names(x))) list(names(x), :
'data' must be of a vector type, was 'NULL'

Retry limit reached; solution not found. Best fit=364650.78 (started at 1881335.6) (11 attempt(s): 9 valid, 2 errors)

Posted on
No user picture. Aivil Joined: 03/14/2019

Error Message - Quadratic Growth Model (Even after mxTryHard)

Hello,

-> Can you help me understanding why I get this error message when fitting this unconditional quadratic growth model?
-> Once running the script attached, I get this error:

Running Quadratic Growth Model, Path Specification with 10 parameters
Error: The job for model 'Quadratic Growth Model, Path Specification' exited abnormally with the error message: fit is not finite (The continuous part of the model implied covariance (loc2) is not positive definite in data 'Quadratic Growth Model, Path Specification.data' row 8664.