OpenMx Help

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

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

mxFactorScores - Error: The KA matrix associated with ... is not an N x 1 matrix

Not sure why I'm getting an error when I'm trying to compute factor scores using mxFactorScores with type="ML". The model has two correlated latent factors.

> fit.lisrel <- mxRun(model.lisrel)
Running Model LISREL
> summary(fit.lisrel)
Summary of Model LISREL

free parameters:
name matrix row col Estimate Std.Error A lbound ubound
1 b_e_nc LX c_nc_el mu_e 1.13324087 0.16178139
2 b_o_nc LX c_nc_or mu_o 1.11548157 0.13534684
3 v_e PH mu_e mu_e 0.11987041 0.03554072 0.01

Posted on
No user picture. Rachel Joined: 02/08/2012

mxTryHard status codes

Dear all,

I have been running a model within a loop over several variables, which does not converge well for all variables. Therefore, I applied mxTryHard. There are two things I noticed:

1) for models that run well without calling mxTryHard the best fit does not seem to change using mxTryHard, but the status code changes from 0 to -1

> fit1 <- mxRun(model,intervals =F)
> fit1$output$Minus2LogLikelihood
[1] 4595.342
> fit1$output$status
$code
[1] 0
$status
[1] 0

fit2<-mxTryHard(model,intervals=F)
> fit2$output$Minus2LogLikelihood
[1] 4595.342
fit2$output$status
$code