I'm interested in finding the interaction impact of 3 genes together on a quantitative phenotype with adjustment of 3 covariates . I want to study the interaction effect but not the main effect. I'm fitting this model:
model<-lm(quantitative trait~ covariate 1+ covariate 2+ covariate 3+gene1:gene2:gene3,data)
So my questions are: Is it correct to look at the interaction effect without looking at the main effect? Is the model correct? Your help is very appreciated.