Login

News

GxE model for binary variable

Posted on
No user picture. diana Joined: 04/17/2019
Forums
Behavioral Genetics Models
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.

Many many thanks!

Replied on Tue, 05/14/2019 - 19:59
Picture of user. AdminRobK Joined: 01/24/2014

What, specifically, are the phenotype and moderator? If you'd rather not say, then are either of the following statements plausible concerning the moderator?:
1. Its covariance with the phenotype is equally due to _A_, _C_, and _E_.
2. All of the variance in the moderator is _entirely_ due to either _A_, _C_, or _E_.

Another question: would you feel more comfortable treating the ordinal moderator as (1) an approximately interval-scale variable, or (2) a categorical variable with unordered categories?

I will say that I know of no OpenMx example script for a binary phenotype and a multicategory ordinal moderator. Indeed, I do not believe I have ever read a paper fitting such a model. In fact, OpenMx may simply not be capable of a completely satisfactory solution to that data-analysis problem....

...Are you familiar with Markov-Chain Monte Carlo Bayesian methods?

Log in or register to post comments
Replied on Fri, 05/17/2019 - 10:40
Picture of user. tbates Joined: 07/31/2009

Have a go with `umx::umxGeE`

Or perhaps `umx::umxGxE_window`

Log in or register to post comments
Replied on Tue, 08/13/2019 - 19:55
No user picture. xiyuesenlinyu Joined: 05/24/2019

hello
I just notice that you are tying to fit a GxE model with a binary phenotype variable, what is i am going to do.

I have read the previous posts and related literature, but still can not change the code by myself.

I would like to ask if your problem has been solved? Could you share your chaged code with me.

Thank your!

Log in or register to post comments
Replied on Thu, 08/15/2019 - 02:48
No user picture. xiyuesenlinyu Joined: 05/24/2019

Thanks so much for your help!
My moderators include both continuous and categorical variables, and the phenotype variable is a binary variable.
I will try the code first. If there are any questions, I'll consult you again.
Thank you again!
Log in or register to post comments
  Login

News