Behavioral Genetics Models
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?
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.
- Read more about Bivariate Cholesky: error in do.call
- 14 comments
- Log in or register to post comments
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.
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?
- Read more about Residuals in ACE model with covariates
- 6 comments
- Log in or register to post comments
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
- Read more about how to model Twins and non twins siblings
- 5 comments
- Log in or register to post comments
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.
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.
- Read more about estimate CI
- 4 comments
- Log in or register to post comments
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) :
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).
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
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