OpenMx Help

Posted on
Picture of user. rabil Joined: 01/14/2010

Handling Exogenous Latent Variables in OpenMx

Suppose you have two exogenous latent variable. (There will be directed paths from them to a response variable.) Is the covariance between the two exogenous latent variables part of the model? That is, do I have to include a path specification for the covariance or otherwise the covariance will be constrained to zero? I understand that the covariances among exogenous manifest variables are not part of the model.

Posted on
No user picture. AshleeM Joined: 01/20/2016

Error in if (class(fit) == "try-error" || is.na(fit$output$minimum) || : missing value where TRUE/FALSE needed

I'm running a measurement model (two time point common factor model) and I am getting an error message I've never seen before when I drop a parameter. I talked to AdminRobK about this, and he suggested I post it to the forums. The error reads:

Error in if (class(fit) == "try-error" || is.na(fit$output$minimum) || : missing value where TRUE/FALSE needed
In addition: Warning message:
In is.na(fit$output$minimum) : is.na() applied to non-(list or vector) of type 'NULL'

Posted on
Picture of user. rabil Joined: 01/14/2010

How can you tell if OpenMx is using multiple cores?

I'm running OpenMx on 4 different computers running Ubuntu 15.04 and the latest version of OpenMx and current versions of R. On one, when look at the output for htop, it seems clear that all 4 cores are being used by OpenMx. Yet on two of the others, htop seems to indicate that not all 4 cores are being used simultaneously and the same model takes 2 to 3 times longer to run.

Posted on
No user picture. Jorien Joined: 02/08/2013

Joint Ordinal-Continuous Model

Hi everyone,

I am trying to run a joint ordinal-continuous model.
In my model I have a continuous variable (varA) and an ordinal variable with 1 threshold (varB).
There are 5 zygosity groups (MZM, DZM, MZF, DZF, DOS) and a covariate age which
I modeled separately for men and women.

The model is actually running perfectly fine (no errors) and the estimates for the thresholds for varB
look good and are what I would have expected.
However, the estimates for the means are all way too low when I compare them with the raw data.

Posted on
Picture of user. tbates Joined: 07/31/2009

access model estimated parameters (EP) without calling summary?

Hi OpenMxers,
I want to determine automatically which order to put models into mxCompare so that the model with more degrees of freedom comes first.

Both models have been run, and I can access a models number of estimated parameters from summary(model).

summary() however, can be costly in terms of execution time: Is there a way to access the EP (or df) property from the model directly without calling summary?

Besten Wünsche,
tim

FYI, the situation I am trying to automate is avoiding this:

|Model        | EP|Delta-2LL |Delta df|p  |      AIC|Compare with Model |
Posted on
No user picture. Cindy.s Joined: 07/04/2015

mxRefModels error

Hi all,

I'm running a Qualitative Sex Differences ACE model and I'm using the mxRefModels to get the reference models, but for some variables I get a error:

>mxRefModels(QualAceFit_rg, run = T)
Error in chol.default(sampcov) :
the leading minor of order 1 is not positive definite
> mxVersion()
OpenMx version: 2.3.1 [GIT v2.3.1]
R version: R version 3.2.1 (2015-06-18)
Platform: x86_64-w64-mingw32
Default optimiser: NPSOL

Posted on
No user picture. Barati Joined: 10/20/2015

executing error, band dispersion problem

Hello everyone,

I am struggling with a “band dispersion” problem.
I am trying to execute my input file with the following command:

% ./openmx input_file_name.dat

and I get the following error:

Error: unexpected input in "% ./openmx input_file_name.dat"

any idea how to fix it?

PS:
I'm working in the directory where my input file is.
I am using the latest version of OpenMX under R 3.2.2 on Win10 (64-bit).

Posted on
No user picture. a.ortega Joined: 02/05/2011

Error message when Installing in server

Hello everyone.
I have tried to install OpenMx for R in a computer cluster running with Linux and I have encountered the following error message:

> cc1plus: error: unrecognized command line option "-Wvla"
> cc1plus: error: unrecognized command line option "-Wno-unused-local-typedefs"
> make: *** [Compute.o] Error 1
> ERROR: compilation failed for package 'OpenMx'

It seems that there is an issue with "gcc"(?), but I'm not sure.

The R version is is 3.2.1 while the Cluster is running on Linux CentOS 5.1.

Any help you may provide will be highly appreciated.

Posted on
Picture of user. rabil Joined: 01/14/2010

Convert RAM Model to LISREL

I need to compute factor scores which only works in OpenMx if the model type is LISREL. I've attached the model code that includes creating a path diagram for the RAM version. The RAM model code works fine in estimating the parameters. When I convert it to "type="LISREL"", I get various error messages.

This shows the specification for the endogenous and exogenous variables in the LISREL version:

manifestVars=list(
endo=c("c_el","c_nc_el","c_or","c_nc_or",

Posted on
No user picture. Lisa M Joined: 08/19/2013

P value of 1 in mxCompare

Hello,
I am working on an analysis of a univariate sex limitation model for depression and I have run into a potential problem. When comparing the homogeneity ACE model to the AE model, I get a p-value of 1 in mxCompare. The fit of the models are very similar:


> mxCompare(HomACEModelFit,HomAEModelFit)
base comparison ep minus2LL df AIC diffLL diffdf p
1 HomACE 4 16217.9 21258 -26298.1 NA NA NA