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 [6], 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!