Behavioral Genetics Models

I am running a multivariate model (3 dependent vars) with sex (female MZ, male MZ, female DZ, male DZ, F&M DZ, M&F DZ) the model runs fine with the SLSQP optimizer.

No SEs are suspect and the return code is zero.

However when I try the command:
FitCholMod <- mxRun(CholMod, intervals =TRUE)

It runs as before but I receive this error message.

Error in runHelper(model, frontendStart, intervals, silent, suppressWarnings, :
NLOPT fatal error -1

So the confidence intervals are not produced.

Hi,

I'm wondering that in which condition we should choose ACE model, or ADE? Is there any criteria to be referred to? For example, in one trait of my own data, the ACE model indicate the CE model is the best fitted model, while the ADE model indicate the DE model is the best fitted model, is this trait heritable? And in the ADE model, how to calculate the heritability? Only the part of A, or the part of A and D, or only the D?

Hi there,

Now I have a problem about CI estimation in a common pathway model. I can not estimate the CI of the path coeffcients of latent factors. I added following codes to the script:

matIl <- mxMatrix( type="Iden", nrow=nf, ncol=nf, name="Il")
invSDl <- mxAlgebra( expression=solve(sqrt(Il*V)), name="iSDl")
sta <- mxAlgebra(iSDl %*% al, name="sta") # standarized path coefficients of A
stc <- mxAlgebra(iSDl %*% cl, name="stc") # standarized path coefficients of C
ste <- mxAlgebra(iSDl %*% el, name="ste") # standarized path coefficients of E

Hello, everyone.

I'm new to behavioral genetics. I've learn some knowledge about this area by myself and finish the univariate genetic analysis about twin data. But it's a little difficult for me to find updated bivariate model script. The newest one I found is written on 2012, some lines cannot be read neither by R or me.

I wish someone can throw a web which contain updated bivariate genetic model script.

Hi all,

I am trying to run a bivariate ACE model. My cross-twin-cross trait correlations are as following:

MZM -.36 (-.44 , -.27) DZM -.23(-.25 ,-.09) DOS -.14(-.21, -.07)
MZF -.24 (-.33, -.14) DZF -.16(-.28 , -.04)

My standardized estimates, however, do not seem to match these correlations. Namely, C is estimated much lower than expected based on the twin correlations (esp. for females).

Males A 0.60(0.38, 0.81) C 0.25(0.07, 0.43) E 0.16(0.14, 0.19)
Females A 0.82(0.48, 1.15) C 0.01 (-.30, 0.30) E 0.18(0.10, 0.27)

Hello

Does anyone have an ACE script for enabling the variance for twins 1 and 2 to be estimated separately (and separately for MZ and DZ twins) to deal with heterogeneity in their variance?

(or even just the adjustments that need to be made to the standard script?)

Thanks so very much

Karen

Hi there,

I am working with twin data. My current analysis involves data that has a lot of variation, and I am using small numbers n=90 MZ twins and n=48 DZ twins. I am working off a univariate ACE script and I am getting the CE model as the best fit for variables where the MZ correlation is significantly higher than the DZ, I would expect some heritable component.

Hi,

I would appreciate some advice on what model I need to use in order to test some specific hypotheses. I have suggested some ways to do this, but not sure if it is doable or correct

I have three traits, with phenotypic correlations of around .40.

Hypotheses:
1- large genetic correlation between trait 1 and trait 2 and 3
2- Much of the genetic variance component of trait 2 and trait 3 can be explained by the interaction between environmental variable B and common genetic factors of trait 1&2&3.

I thought I could maybe try the following analyses:

Hello! This is a follow-up post to some guidance I got on this board awhile ago that was incredibly helpful. I am trying to run a bivariate saturated model that controls for age and sex. The two variables were assessed at different time points, so one variable needs to be adjusted for age1 and sex, the other variable for age2 and sex.