OpenMx Help

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.

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.
- Read more about Joint Ordinal-Continuous Model
- 4 comments
- Log in or register to post comments

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 |

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
- Read more about mxRefModels error
- 7 comments
- Log in or register to post comments

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).
- Read more about executing error, band dispersion problem
- 1 comment
- Log in or register to post comments

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.
- Read more about Error message when Installing in server
- 1 comment
- Log in or register to post comments

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",
- Read more about Convert RAM Model to LISREL
- 9 comments
- Log in or register to post comments

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
- Read more about P value of 1 in mxCompare
- 2 comments
- Log in or register to post comments

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

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
- Read more about mxTryHard status codes
- 9 comments
- Log in or register to post comments
Pagination
- Previous page
- Page 18
- Next page