Behavioral Genetics Models

Posted on
No user picture. martonandko Joined: 02/19/2015

Obtain chi-square fit index of univariate ACE model

I'm trying to get chi-square fit index of my univariate ACE model. I tried to use mxRefModels(ACEFit) and then use mxCompare between the mxRefModels output and my ACE model (ACEFit), but it returned NA for the p value

base comparison ep minus2LL df AIC diffLL diffdf p
Saturated twinACE 0 NA 0 NA NA NA NA
Saturated twinACE twinACE 6 879.7042 178 523.7042 NA 178 NA
Independence twinACE 0 NA 0 NA NA NA NA
Independence twinACE twinACE 6 879.7042 178 523.7042 NA 178 NA

Posted on
No user picture. Yuan Joined: 02/13/2015

Running CholAE

When I was running the CholAE model, it warns that:
In model 'CholAE' Optimizer returned a non-zero status code 1. The final iterate satisfies the optimality conditions to the accuracy requested, but the sequence of iterates has not yet converged. Optimizer was terminated because no further improvement could be made in the merit function (Mx status GREEN).

What does this mean?

Thank you

Posted on
No user picture. slukow Joined: 02/25/2015

Calculation of BIC discrepancy?

Hello - I'm having some trouble with discrepancies in the calculation of BIC in comparing output from Mx with OpenMx on both univariate and bivariate ACE models. I am currently using OpenMx 2.0.1 and R version 3.1.2. Below is the univariate script I used. It seems from the output that the problem may be coming from the fact that mxData objects are pulling number of observations from the total number of rows in mzData and dzData, which includes some rows in which both variables are missing.

Posted on
No user picture. Yuan Joined: 02/13/2015

bivariate genetic model

Hi,
I am new to Twin data analysis and would like to get help on the bivariate genetic analysis. I want to calculate the genetic correlations between two traits in twins. However, I only found the scripts about the univariate twin analysis in documentations. How can I get the bivariate ACE model in OpenMX?

Thank you

Posted on
No user picture. GRBS Joined: 02/13/2015

Constrain correlation matrices

I have fit two saturated models and have created correlation matrices from the expected covariance matrices using the cov2cor function. I would like to constrain the correlation matrices to equality between these models without constraining the covariance matrices.

Posted on
No user picture. K Ram Joined: 08/21/2012

Correlated factors ok to use even when phentypic varaiance not significant?

We have 7 variables, of which one of the variable is significantly heritable - modeled using univariate analyses. Is it theoretically congruent to test for shared genetics between all 7 variables using the correlated factors model despite the fact that 6 of the variables are not significantly heritable?

Posted on
No user picture. jarrode28 Joined: 03/01/2010

Correlated Genetic Factors in Independent Pathways Model

Hi,

I'm attempting to run a BG script of an independent pathways model on 9 dichotomous (symptom) manifest variables. That is, the model includes biometric factors specific to each manifest variable and common to all manifest variables. I've been able to get this model to run without any problems.

Posted on
No user picture. luna.clara Joined: 12/23/2012

Theoretical question

Hi all,
I have a theoretical question.
In what ways apparent qualitative and quantitative sex differences are or are not addressed by simply adjusting for sex in a multivariate model?
The thing is that I observed in different univariate models that 3 of my 6 variables (for the multivariate model) had quantitative or qualitative sex differents... then I adjusted for sex and age in the multivariate model, but I'd like to explain in which way those genetic sex differences in some of the variables are addressed or not in the multivariate model just by adjusting for sex. Does it make sense?

Posted on
No user picture. luna.clara Joined: 12/23/2012

How to add a second latent common A factor in an independent pathway model

Dear all,
I have run without problems a multivariate IP model for 6 variables.
The thing is that now I would like to add another latent common A factor that loads only on 3 of the 6 variables to see if this "new" 2-IP factor model (two Ac factors) fits better the data than the 1-IP factor model (just one Ac factor for all the observed variables).
But I think I’m doing something wrong because the model doesn't understand that the second common A factor is loading only in the first 3 observed variables. Can anybody help me? Many many thanks!

Posted on
No user picture. YiTan Joined: 10/08/2014

A question about moderator variable

Hi all,

I'm new to OpenMx and SEM, and I have just read through the various examples on twin modelling. Just wondering, if I have a moderator variable which varies within twin 1 and twin 2 and also across twin pairs, how do I go about specifying it as a definition variable? Or do I have to do bivariate analysis instead (treating the moderator variable as another variable?)

Thank you very much!

Cheers,
Yi