OpenSEM Forums

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)