OpenMx General Help
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**:
hansfredriksunde
Joined: 12/10/2021
Wrong estimates when using definition variables and WLS
Hi everyone!
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
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.
- Read more about Estimating a non-linear model in OpenMx
- 5 comments
- Log in or register to post comments
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
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):
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:
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.
- Read more about Automation of monophenotype analyses
- 1 comment
- Log in or register to post comments
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.
- Read more about Bivariate ACE model with moderator
- 7 comments
- Log in or register to post comments
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