The editor has been updated to markdown! Please see more info at: Tutorial: Updated Support Site Editor

Question: EdgeR model matrix: interaction term
gravatar for fcamus
11 months ago by
fcamus0 wrote:

Hi all,

I'm in the process of analysing my RNASeq data and had a couple of question about setting up the model matrix to get the comparisons I want. 

First of all my experimental setup is looking at genes that respond to diet in males and females. 
> Protein-Female
> Carbohydrate-Female
> Protein-Male

the model matrix I have setup is:
> des1

   Fem.P Male.P Fem.C Male.C
1      1      0     0      0
2      1      0     0      0
3      1      0     0      0
4      0      1     0      0
5      0      1     0      0
6      0      0     1      0
7      0      0     1      0
8      0      0     1      0
9      0      0     0      1
10     0      0     0      1
11     0      0     0      1

 to look at differences between the diets within each sex I have done the contrast:

FemC-FemP and MaleC-MaleP in separate analyses. 

I am interested in finding out the genes that do different things in males and females in response to diet (sex x diet interaction)
What contrast do I do to get to this question? 
I was thinking something like:
makeContrast((Fem.C-Fem.P)-(Male.C-Male.P), level=des1)

is this correct?






edger • 238 views
ADD COMMENTlink modified 11 months ago by Aaron Lun22k • written 11 months ago by fcamus0
Answer: EdgeR model matrix: interaction term
gravatar for Aaron Lun
11 months ago by
Aaron Lun22k
Cambridge, United Kingdom
Aaron Lun22k wrote:

Yes, that's correct. Your last contrast will identify the difference in the diet effects between sexes, i.e., the interaction between diet and sex. If no interaction existed, sex and diet should be purely additive, such that the effect of diet in females (Fem.C - Fem.P) should be the same as that in males (Male.C - Male.P).

Note that interpretation of the interaction often requires that you report the individual diet effects in both sexes as well, just to see whether the interaction is due to decreased/opposite/elevated effect in males relative to females.

ADD COMMENTlink modified 11 months ago • written 11 months ago by Aaron Lun22k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 390 users visited in the last hour