Behavioral Genetics Models

Posted on
No user picture. diana Joined: 04/17/2019

How to fix the genetic correlation larger than 1

Hi everyone, I am fitting a bivariate ACE model and I am very confused about the Ra (genetic correlation of two traits) that is larger than 1. Since Ra is a correlation coefficient, it should range from -1 to 1, am I right? I don't know how to explain this result.

In the ACE model, Ra=1.37 ,Rc=-0.04, Re=0.22.
Then I fitted ADE model, Ra=0.07, Rd=NA, Re=0.22.

The results of the goodness-of-fit test between saturated model and ACE or ADE are as follows.

Posted on
No user picture. qingwen Joined: 08/15/2019

CI is confusing

I fit a nuclear twin family model. AE model is the best. However, the CI of e makes me confused. Here are the output in free parameters and CI details section. The CI details looks strange, ranging from 0 to 2.053589e+06.

Posted on
No user picture. liagym1 Joined: 11/06/2019

ACE with moderator-help interpreting results

I am using regular Mx (unfortunatly i cannot find an active forum), and would like to clarify the meanings of the output received after running in an ACE with an ordinal moderator script.
1. Are the estimates of T,U,V the moderated components? (often denoted by betas)
2. If the estimate for a certain parameter is not zero, yet the CI contains zero- does this imply that this component is not significant?
Posted on
No user picture. liagym1 Joined: 11/06/2019

Error in output

I am trying to run an ACE with moderator script, and this is the error I received in the output. What changes should i make in my data?

Error found in rectangular data; record number 1 :

 1004 -999 -999 -999 -999 -999

We have a problem whose error code is 176

and which I ran across at line number 38

of your input script

RE FILE=AGE11DATA.DAT

~~~

An error has occurred while reading a rectangular file

Make sure that you *don't* have a FORMAT at the beginning

Posted on
No user picture. ReddyM Joined: 10/02/2019

Difference between optimizers

In my twin model i calculate Rmz and Rdz from the saturated model, however when i run it (i run it on several variables) for a few variables CSOLNP fails to calculate confidence intervals, while SLSQP does caculate all. Is there a particular reason for this? Should i be concerned by this, or just use slsqp for the cases csolnp fails?

Here is the sat part, if it helps:


#sat model, rmz, rdz:

PathMZ <- mxMatrix("Lower",nrow = 2,ncol = 2,free = TRUE,values=starting_value_pathcoeff,name = "PathMZ")#cholesky for pos definitness

Posted on
No user picture. lior abramson Joined: 07/21/2017

direct comparison between traits' heritability with no raw data

Dear forum
I am trying to conduct a meta-analysis on the heritability of two traits that are different, yet related (I found six studies in which the twins were measured in both traits).

From what I know, in order to directly compare between the two traits heritability, I need to enter them in the same model, and for that, I need either the raw data or the full correlation matrix (and I have neither).

Posted on
No user picture. George Richardson Joined: 10/18/2018

Multivariate ACE Vs. Multilevel Model

Hi All,
I'm trying to understand the relationship between multivariate ACE models (e.g., Cholesky model fit to two phenotypes) and multilevel behavior genetic models. Is there a formal comparison somewhere?

My sense is that multivariate ACE models conflate within and between pair variance. A and C could capture within and between effects while E captures the within pair effects. E will map right onto the within pair effect in the multilevel model (MLM). Is all this correct?

Posted on
No user picture. qingwen Joined: 08/15/2019

nuclear twin family design in OpenMx

Hi, all

I am trying to fix a Nuclear Twin Family model in OpenMx.

Genetic and environmental parameter estimates in my model are obtained from four pieces of information, the covariance between MZ twins, the covariance between DZ twins and the covariance between parents and the covariance between parents and children.

Posted on
No user picture. Aurina Joined: 01/05/2016

Assumptions bivariate ACE twin model

Hello,

I am planning to do bivariate ACE modelling using IQ and brain imaging derived measures and I would be very grateful if you could help me answer a few questions regarding the general assumptions for this type of modelling.

Posted on
No user picture. JuanJMV Joined: 07/20/2016

Saturated vs ACE model- different correlations

Hi all,

I am trying to fit a univariate model.

I fitted the saturated model first and I got correlations of 0.64 and 0.34 for MZ and DZ respectively.

However, when I fitted the ACE model, I got these results:
A=0.04
C=0.52
E=0.44

So, I decided to check the correlations from the ACE model and I got 0.56 and 0.54 for MZ and DZ respectively.
I know that correlations may change from the ACE to the saturated model. However, these differences are huge.