Behavioral Genetics Models
Benny
Joined: 05/15/2020
valentinav
Joined: 06/15/2020
ACE output: what is the p
I run the ACE model and got a summary table which reports the p among others.
The p is also not clear to me. Apparently it is higher for the winning model, but what exactly is this p?
Thanks,
Valentina
- Read more about ACE output: what is the p
- 1 comment
- Log in or register to post comments
Leo
Joined: 01/09/2020
ACE RAM model with ordinal manifest
I've been trying everything to get an ordinal manifest working in RAM models but I can't get it done. Also, didn't find anything online about this, as most examples have matrix specification. Would be very appreciated to get some help. Besides from the solution below, I have also tried every possible way with mxthreshold and umxthresholdmatrix.
- Read more about ACE RAM model with ordinal manifest
- 12 comments
- Log in or register to post comments
valentinav
Joined: 06/15/2020
really similar AIC
I run the ACE model and submodels, and the AIC are really similar.
ACE 2579
AE 2577.5
CE 2577.3
E 2584
Overall, I always get a difference of 1 between the AIC.
The best winning model is always coherent with the results of the correlation analysis, so I tend to trust it, however, the difference is really little.
Is there a way to establish what is a significant difference between AIC?
Thank you all!
Valentina
- Read more about really similar AIC
- 7 comments
- Log in or register to post comments
lemire
Joined: 12/13/2019
ACE: negative variance component from H.Maes’ published script and slides
could you help provide some insights as to why running the following script, which uses data(twinData), gives negative estimates for Vc in an ACE model?
https://ibg.colorado.edu/cdrom2020/maes/twinModeling/oneACEvc.R
These are the most relevant lines related to the model:
inooradd
Joined: 03/11/2020
Step 1: Testing basic assumptions (equal means & variances for twin 1/twin 2 and MZ/DZ pairs)
I am following the steps laid out in [this](https://vipbg.vcu.edu/media/course/HGEN619_2015/TwinModels.pdf ) openmx presentation and could use help with step 1.
1. Use data to test basic assumptions (equal means & variances for twin 1/twin 2 and MZ/DZ pairs)
valentinav
Joined: 06/15/2020
ACE model lower CI negative number
I am new to ACE model estimation so I am trying to follow some tutorials I found online (script attached).
The results seems fine, and are also coherent with another script I used. However, I got negative lower CI: ACEvc.VarC[1,4] -0.2443061 0.23255753 0.7557620.
Does this make sense to you? As far as I know when CI include 0 we cannot conclude that the value is different from 0...
I would really appreciate whether someone could help me with this.
Thank you in advance!
Best,
Valentina
- Read more about ACE model lower CI negative number
- 7 comments
- Log in or register to post comments
inooradd
Joined: 03/11/2020
Comparing ACE and ADE models
I have a couple of questions regarding comparing ACE and ADE models.
1. Is there a quick method of comparing the fit of both of these models? As of right now I am running the scripts for each model, then comparing the two models using mxcompare:
<%
mxCompare( fitADENOTW , fitACENOTW ) %>
I know typically when we compare nested, we run something to the effect of...
- Read more about Comparing ACE and ADE models
- 7 comments
- Log in or register to post comments
Benny
Joined: 05/15/2020
GxE with latent moderator
- Read more about GxE with latent moderator
- 2 comments
- Log in or register to post comments
Leo
Joined: 01/09/2020
Bivariate Cholesky model with one continuous and ordinal variable
Particularly, I am trying to fit a model with one continuous and one ordinal variable. I have found this script but it is modeling only two ordinal variables:
http://ibg.colorado.edu/cdrom2016/maes/UnivariateAnalysis/twoa/twoACEoa.R
Is there a shortcut to only apply the threshold to one variable?
Thanks in advance.