OpenMx General Help

Posted on
No user picture. lf-araujo Joined: 11/25/2020

Cascade / stealth model script example for OpenMx

Hi all,

I hope everyone is healthy.

I have been looking into extended family design, starting with S Medland papers [1,2]. Also, M Keller have some basic scripts [here](http://www.matthewckeller.com/html/cascade.html) which does not specify Medland's model completely, nor the Stealth model. My question:

- Does anyone care to share a script example for training purposes on either/or models? Or point to where I can find these?

Many thanks,

Luis

**Edit**:

Posted on
No user picture. MaximilianStefan Joined: 12/04/2018

Performance regression of OpenMx for simple CFA model

I specified a simple CFA model (5 factors, 20 Items each), simulated data from it (25 observations for each parameter) and fitted the model in lavaan and OpenMx. Surprisingly, OpenMx takes about 100x longer to fit the model (same starting values in lavaan and OpenMx). I attached an MWE.
The results on my machine are
lavaan: 0.16 seconds/31 iterations
OpenMx: 13.2 seconds/135 iterations

Posted on
No user picture. pgok Joined: 11/13/2019

Estimating a non-linear model in OpenMx

Hi, I'm pretty new to OpenMx, having only used it off-and-on over the past few years. That said, I am pretty interested in using it more, especially given the multilevel modeling capabilities (e.g., as discussed in Pritkin et al., 2017). However, my question isn't about the multilevel modeling capabilities, my question is about fitting non-linear models.

Posted on
No user picture. cjvanlissa Joined: 04/10/2019

Get all log likelihoods from mxTryHard() extratries

Dear all,

is it possible to obtain the log likelihoods from all tries of mxTryHard()? It prints a lot to console, but I can't find this information in the returned object.

Thank you kindly,
Caspar

Posted on
No user picture. IvanVoronin Joined: 08/18/2013

CIs are NA but their status is OK - is this OK?

Hello OpenMx community,

I am trying to compute confidence intervals for the parameters that I define as matrix algebra
summary(model) gives the following (only CIs showing):

Posted on
No user picture. MaximilianStefan Joined: 12/04/2018

Parallelization of the objective function for definition variables

We are currently trying to benchmark our SEM software against OpenMx for a SEM with unique definition variables per person.
With ML estimation, because there is a unique model-implied covariance matrix per person, which has to be inverted per person, I assumed that parallelizing the objective function should improve performance drastically. However, changing the number of threads does not change the performance - so I assume I probably did something wrong. I followed this wiki pages instructions:
Posted on
No user picture. sar2259 Joined: 11/10/2020

Automation of monophenotype analyses

Hello,

please could you let me know whether it is possible to automate monophenotype analyses ( and also when adjusting for age and sex)?

I was searching whether this might be possible, and found one similar question (https://openmx.ssri.psu.edu/node/4673) with very similar code to what I am currently using. However, it is unclear whether it is possible to do this or not. I have 142 monophenotypes I need to calculate heritability for.

Thank-you.

Posted on
No user picture. EmilieRH Joined: 09/24/2020

Bivariate ACE model with moderator

Hi all,

I am currently working on a project on how school achievement moderates intelligence among young adults. My sample consists of 4,084 German twins aged 17 and 23-24, respectively, who have information on both GPA and intelligence. However, the twins' GPAs come from various types of secondary schools so to check whether I need to use school type as a moderator in my GxE analyses, I first want to see what happens if I run separate bivariate variance decomposition analyses for the three school types constraining their parameters equal.

Posted on
No user picture. Soldor880 Joined: 10/14/2020

Applying the univariate ACE model to more than one phenotype at a time

I have 75 metabolites. I am running a univariate ACE model analysis for each one to decompose the phenotypic variation into estimated genetic and environmental contributions.

Is it possible to adjust the code attached to run the model for each of these simultaneously to get the output (h, c and e), or is it only possible to run the model for each metabolite individually?

Thank-you