Behavioral Genetics Models

Posted on
No user picture. JuliaJoplin Joined: 10/12/2014

Bivariate Cholesky with Covariates - Best Practice?

I'm doing a bivariate Cholesky decomposition between two variables that were measured at two different time points. I know that it's normally best practice to adjust variables for age and sex effects, and I see many different threads on this board that have discussed how to add covariates to an ACE model. Two questions:

1. For my purposes, should I regress out sex and age at time point 1 for the first variable and then regress out sex and age at time point 2 for the second variable?

Posted on
No user picture. Daxide Joined: 02/15/2016

Bivariate Cholesky: error in do.call

Hi! I am a novice, trying to move from univariate to multivariate analysis. So I downloaded the script cholesky.R from Sarah Medland's webpage (https://www.genepi.qimr.edu.au/staff/sarahMe/Workshop13.html), along with the data file. I tried to adapt the script to my data but I ran into an error. So I ran the script as I downloaded it, on the supplied simulated data, without modifying it (see eattachment) and I still encounter the same error, at this point.

Posted on
No user picture. k.j.kan Joined: 01/29/2016

Multivariate heritability in an independent pathway model

Dear Guinea-pig lovers,

Heritability is defined on the basis of phenotypic variance. One may regard a so-called common pathway model as a model that contains one or more latent phenotypes that mediate certain genetic and environmental effects on observed variables that function as indicators of those latent phenotypes. As a result the common pathway model enables us to calculate the heritability of both the mediating latent phenotype as well as the observed phenotypes.

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

Residuals in ACE model with covariates

Hello all,

I am really new in twin modeling as well as in R, therefore I would be very grateful if someone could advise me on a few issues regarding the script.

I am trying to run an ACE model with covariates (sex and age) on some twin height data. It would be very helpful to know whether the script (attached: test_height.R) does not have any major flaws and does what is expected of it. And one more thing, I would like to plot the residuals after performing the regression on the data. Could you please help me to find the most straight forward way to do that?

Posted on
No user picture. yassinebha Joined: 01/02/2016

how to model Twins and non twins siblings

Anyone knows how to use OpenMx for twins study including non twins siblings, like family model?
Thanks

Posted on
No user picture. jbizzle09 Joined: 10/02/2015

Bivariate model for ordinal and continuous data

I would like to run a bivariate model for an ordinal variable and a continuous variable. The ordinal variable was measured in childhood and the continuous variable was measured in adulthood - therefore it makes sense that I should use a Cholesky decomposition testing whether the ACE influences on the ordinal (childhood) variable also influence the continuous (adulthood) variable. However in a liability threshold model for continuous and ordinal data, the ordinal variable should go last.

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

estimate CI

Hi everyone,

I'm recently learning the univariate twin analysis. According to the documents, I got the a2, c2, and e2. How could I estimate the confidence intervals for each value? Thanks.

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

Equating covariances as a test for sex-limitation effects?

Hi all,

I have a question about sex-limitation model testing. I came across an article which mentioned the following (http://www.ncbi.nlm.nih.gov/pubmed/15685427) :

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

Bivariate analysis: why best-fitting model's ACE estimates don't correspond well with the univariate estimates?

Hi all,

I have a question about bivariate cholesky ACE model. For model fitting, I compared the reduced models (selectively dropping a21, c21 and e21) to the full bivariate ACE model and the best-fitting model was selected based on lowest AIC and nonsignificant likelihood ratio chi-square test.

However, I noticed that the best-fitting model’s ACE estimates do not seem to match with what I obtained from fitting each of the two variables in univariate ACE models. For example, in univariate model fitting I got:
First trait: C = 60% (CI: 0.46-0.71) and E = 40% (CI: 0.29-0.54).

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

Question about univariate ACE model selection

Hi all,

I have a quick question to ask about model selection.

Let’s say I am interested to find out the ACE for years of music training, and for the full model I got the following estimates: A = 51% (CI: 0.18-0.89), C = 33.6% (CI: 0-0.65), E = 15.4% (CI: 0.10-0.24).
The fit statistics of the submodels are as follows:


base comparison ep minus2LL df AIC diffLL diffdf p
1 HomoAce 6 1022.900 193 636.8998 NA NA NA
2 HomoAce AE 5 1024.610 194 636.6104 1.710576 1 1.909107e-01