Behavioral Genetics Models

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?

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.
- Read more about Bivariate Cholesky: error in do.call
- 14 comments
- Log in or register to post comments

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.

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?
- Read more about Residuals in ACE model with covariates
- 6 comments
- Log in or register to post comments

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
- Read more about how to model Twins and non twins siblings
- 5 comments
- Log in or register to post comments

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.

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.
- Read more about estimate CI
- 4 comments
- Log in or register to post comments

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) :

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).

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
2 HomoAce AE 5 1024.610 194 636.6104 1.710576 1 1.909107e-01
- Read more about Question about univariate ACE model selection
- 1 comment
- Log in or register to post comments
Pagination
- Previous page
- Page 21
- Next page