Behavioral Genetics Models

Posted on
No user picture. thughe01 Joined: 10/18/2009

ACE model with technical replicates

I'm re-entering the world of SEM after a hiatus and slowly familiarising myself with OpenMx (previous Mx user). I've successfully implemented some simple uni- and bi-variate VC models for some of our latest twin data (microbiome phenotypes) using the umx packages.

Posted on
No user picture. noamm Joined: 10/10/2018

Failure to obtain CI's for a umx model

Hi all,

I have used umx for a bivariate sex-limitation model, using the umxSexLim function, which have been a great help. However, I've failed in trying to compute CI's for my estimates, or any other indication of significance. Attached below are all the options I've tried. I would appreciate any help in how to compute any indication of significance.

Thanks!

Posted on
Picture of user. chancey Joined: 09/09/2017

Modeling Suggestions for Experimental Data

Question for the group: How might one go about modeling experimental data with twin respondents?

I ran a single-factor survey experiment with a small group of twin respondents (MZs and DZs, mainly same sex). There were two conditions. I randomly assigned twins to receive either:

Conditions A A
Conditions B B
Or
Conditions A B (or B A)

Posted on
No user picture. JBosdriesz Joined: 03/07/2017

Age as covariate or moderator

I looked around the forum, and found referrals to a post that might address this issue, but that page wouldn't load anymore, so I thought I'll post my question here.

I'm running a set of univariate ACE models, on a range of outcomes. I have two age groups (3-4 and 7-8), that I would like to analyse separately, by also combined.

Posted on
No user picture. camus Joined: 03/21/2019

Obtain 95%CI for ACE model- script revised from Dorothy Bishop

Hi, I am new with R and OpenMx. Recently I revised Dorothy Bishop's script to run ACE models. However, her original code could not get the 95%CI for of A,C,E estimates, but only 95%CI for Unsquared path estimates. I know mxCI() can obtain 95%CI, but after various attempts, it never works for me. Could any expert have a look of the code and let me know how should I add it?
Thanks!
Posted on
No user picture. diana Joined: 04/17/2019

I am wondering if define variables will decrease E and increase A?

I am wondering if too many define variables will decrease E and increase A?

For instance, I want to estimate the heritability of coronary heart disease (CHD), should I adjust all the risk factors of CHD as much as possible, such as smoking, drinking and BMI? I'm afraid that it will overestimate A and underestimate E since I saw a revelent [paper](https://www.ncbi.nlm.nih.gov/pubmed/15520515/).

Many thanks!

Posted on
No user picture. Qiuzhi Xie Joined: 12/27/2018

About the overlap on the genetic and shared environmental factors

I am conducting the Cholesky model and wondering how to judge if two phenotypes significantly overlap on a genetic or environmental factor. It should be judged by the loading factor and the 95% CI, right? In the following example, only the third phenotype (P3) does not significantly overlap with the others on the A1 (genetic) factor, as the lower boundary of the 95% CI is smaller than .00. Am I right?

Posted on
No user picture. diana Joined: 04/17/2019

GxE model for binary variable

Hello! I am trying to fit a GxE model. The phenotype is a binary variable, and the moderator is measured on a 6-point scale .

I want to use the code posted in [there](http://ibg.colorado.edu/cdrom2016/franic/Moderation/Practical/DO_NOT_OPEN_Moderating_covariances_practical_solved.R), but this code is suit for continuous variable and I don't know how to change it to be suit for binary variable. Could you tell me what lines I should change or add?

Because I'm pressed for time, I am sorry to directly ask for solution instead of learning by myself.

Posted on
No user picture. JuanJMV Joined: 07/20/2016

Multivariate ACE-ADE model

Hello,

I am trying to fit a 5-variate model. But I am not sure about the best way to do it.

I have got these results from the univariate analyses

VAR1
ADE 0.34 0.10 0.57
AE 0.42 * 0.58

VAR2
ADE 0.26 0.16 0.58
AE 0.41 * 0.59

VAR3
ACE 0.29 0.13 0.58
AE 0.45 * 0.55
CE * 0.33 0.67

VAR4
ADE 0.23 0.17 0.59
AE 0.39 * 0.61

Posted on
No user picture. diana Joined: 04/17/2019

Why mxTryHardOrdinal()dosen't report start values?

Hi, everyone!
When i useSatFit <- mxTryHardOrdinal(SatModel, intervals=TRUE,extraTries = 15,finetuneGradient=FALSE,bestInitsOutput=TRUE) ,the result is Solution found! Final fit=22505.957 (started at 22819.562) (16 attempt(s): 16 valid, 0 errors). I don't know why it dosen't report start values from the best fit.