Behavioral Genetics Models
lf-araujo
Joined: 11/25/2020
Cannot invert?!
Hi all,
I am getting a cannot invert error in attempting to run a CLPM in twins:
Error in runHelper(model, frontendStart, intervals, silent, suppressWarnings, : MZ.data: weight matrix is rank deficient; cannot invert
I suppose there is something wrong with my identity and beta matrices. Can someone please advise? The minimal example that gets you there:
- Read more about Cannot invert?!
- 1 comment
- Log in or register to post comments
JuanJMV
Joined: 07/20/2016
Bivariate model — non-linear association
Hi All,
I am working on a bivariate model for two continuous variables (normally distributed).
I have checked the association between these variables and looking at the plots the relationship seems to be non-linear. I have run a regression model and the association is significant but it is probably better explained by a quadratic or U-shaped association.
- Read more about Bivariate model — non-linear association
- 4 comments
- Log in or register to post comments
Ilaria
Joined: 07/22/2024
Multivariate code ACE models (path and variance based models)
Hi all!
I have two questions regarding the script I am using to fit a multivariate model exploring the association between body dissatisfaction (exposure X1) and depressive symptoms at ages 21 and 26 (outcomes Y1 and Y2) whilst accounting for twins' sex and age (residuals being regressed).
I have tried to adapt Hermine's code but I am not sure this is correct. I do get estimates and correlations for all the parameters (X1, Y1, Y2).
However, I am not sure about the following:
Zhiyang
Joined: 12/11/2024
Very small estimate with a very large standard error in univariate modelling.
Dear all,
This is the first time to have a post in this forum, and thanks for the platform.
I am a researcher working on the Finnish Twin Cohorts, and I am doing a univariate twin modeling on depression. I have 113 DZ pairs and 69 MZ pairs. The outcome of depression is basically normal-distributed.
Then in the ADE model, I get the data summary and results like this:
Stefano5
Joined: 06/30/2023
JuanJMV
Joined: 07/20/2016
How to add object to a model without re-run
I have fitted a very computationally demanding model (around one week) and now I need to add an mxAlgebra object into this model:
threG2 <- mxAlgebra( expression= threG[,1:7], name="threG2" )
I know I can do it in this way:
modelIPAE2 <- mxModel( fitIPAE,threG2, name="mulIPAE2c" )
fitIPAE2 <- mxTryHardOrdinal( modelIPAE2, intervals=F, extraTries=11 )
But I wonder if there is any way to modify the model without re-running it
- Read more about How to add object to a model without re-run
- 7 comments
- Log in or register to post comments
tvo015
Joined: 01/18/2024
Van der Sluis in OpenMX?
- Read more about Van der Sluis in OpenMX?
- 1 comment
- Log in or register to post comments
JuanJMV
Joined: 07/20/2016
Fix mean latent variable
I am working on a CP model (using a reference script from Boulder). I would like to know if it is possible to fix (or freely estimate) the mean of the latent variables (here the latent factor) using a model specification like this one (see below). Would it be possible by fixing one of the manifest means? I know it is possible (and must be done) to fix the variance to one so that the model is identified but I am also curious if the mean of latent variables can be fixed (or freely estimated). If so how could I do it?
- Read more about Fix mean latent variable
- 4 comments
- Log in or register to post comments
tvo015
Joined: 01/18/2024
vjvuorio
Joined: 01/01/2024
Multivariate moderation
I'm an undergrad working on my dissertation, and I've managed to build my first models thanks to these forums, the Boulder Workshops, as well as [Hermine's wonderful stash of scripts](https://hermine-maes.squarespace.com/). I have several questions, however, just to make sure I'm doing things the right way, and if what I want to do is possible.
- Read more about Multivariate moderation
- 3 comments
- Log in or register to post comments