Behavioral Genetics Models

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

Multiple testing

One of my Phd students recently received reviewer comments for a submitted paper in which she presented 9 univariate ACE models on separate outcome measures. One of the reviewers asked: “Please indicate how correction for multiple comparisons was handled for the genetic modeling. Given that several brain regions were being analyzed, what statistical threshold was used?”
First of all, is any such correction really needed for 9 outcomes?
Posted on
No user picture. JuanJMV Joined: 07/20/2016

doubt with binary variables

Hi everyone,

I am trying to do a joint analysis between 2 binary variables and 1 continuous variable I am having some problems to fix the variance for the 2 binary variables equal 1.

I think I should change something here

matUnv <- mxMatrix( type="Unit", nrow=nvo, ncol=1, name="Unv1" )
matOc <- mxMatrix( type="Full", nrow=1, ncol=nv, free=FALSE, values=oc, name="Oc" )
var1 <- mxConstraint( expression=diag2vec(Oc%&%V)==Unv1, name="Var1" )

NVO=2
NV=3

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

Rationale for Bootstrapping CIs

Hi all,

Is there a rationale for bootstrapping 95% CIs for univariate ACE estimates when relying on a relatively small convenience sample of MZ/DZ twins? Any source(s) that would support doing so?

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

Moderation Model

Hi everyone,

I am trying to fit a moderation model. I have modified one script but in this script the moderation variable is the same for twin1 and twin 2 (age). I want to fit a moderation model with 1 variable with different values for twin1 and twin2 (like IQ). I have done all the changes but I am not sure if it is correct. Specially this line:

meanG <- mxAlgebra( expression= cbind((int + DR1),(int + DR2)), name="expMeanG" )

Posted on
No user picture. ECTore Joined: 11/06/2017

Direction of causation model

Hello,

I have a question regarding the Direction of Causation (DOC) model. I have found a very useful script at the following link and I'm trying to modify it according to my dataset: http://ibg.colorado.edu/cdrom2016/verhulst/Causation/DOC%20student.R.

I have 2 observed variables (x and y, both continuous) and I would like to verify if their association is causal (i.e. if x causes y. The inverse direction is not possible as x comes before y).

Posted on
No user picture. Fengxia Joined: 11/02/2017

Bivariate Cholesky - Continuous & Ordinal some errors

Hello,

I am trying run bivariate cholesky model analyses with my data in OpenMx2.5.0 recently. But I've got a couple of questions:

1)When I ran the code "CorFit <-mxRun(CorModel, intervals=F)", I got error message"Error in eval(result, envir = env) : object 'MZ.rMZ' not found".

Posted on
No user picture. JuanJMV Joined: 07/20/2016
Posted on
No user picture. likelycandidate Joined: 09/28/2017

Exclamation points indicating

I am estimating twin data using the standard model (http://openmx.psyc.virginia.edu/docs/OpenMx/2.7.9/GeneticEpi_Path.html).

Posted on
No user picture. koak Joined: 05/10/2012

Deleting attachments from posts

Hello

Does anyone know how I can delete attachments that I put on replies to posts on the forum? There are a couple I would like to delete, but at the moment cannot even view them in the reply posts (so no attachment is up with the post).

Thank you for your assistance

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

Another CI problem w/ ordinal ACE model

Hi all,

It appears my 95% CI problems for a series of ordinal ACE models I'm running is not solved after all. Specifically, while CIs for my AE and CE submodels appear to be fine, the ACE CIs return an error (shown in bold below).