Search
Question: 2 factors interaction with Deseq2
0
gravatar for negroni
14 days ago by
negroni0
negroni0 wrote:

Dear Michael,

 

I have a full factorial  experimental design with 2 factors that have 2 levels each, a total of 4 combinations.

I am interested in the interaction between those two factor and that is the reason why I used the following command:

dds <- DESeqDataSetFromMatrix(countData= reads,colData = Design, design = ~ Factor1 + Factor2 + Factor1: Factor2).

Afterward I got a list of contigs that are significantly interacting between the two factors.

My questions are:

1) is there a way to test within the same model the independent effect of each factor or shall I rerun another model without the interaction?

2) with this complete model I can get lists of differentially expressed contigs from the comparison between two level of one factor within each level of the other factor separately and reciprocally (total of 4 comparisons). This looks like a post hoc multiple comparison test (which is allowed only if the interaction is significant). So shall I substract from those lists of contigs the ones that are absent from the list of interacting contigs or not?

 

Than you for your help.

 

Matteo

ADD COMMENTlink modified 9 days ago by Michael Love18k • written 14 days ago by negroni0
0
gravatar for Michael Love
9 days ago by
Michael Love18k
United States
Michael Love18k wrote:

(1) Depending on what you mean by independent, then no. The interaction model means that the fitted effects for factors 1 and 2 take into account their interaction (you obtain the effects for the reference level of the other factor). See the interactions diagram in the DESeq2 vignette. 

(2) Different people approach this different ways, but my preference would be to present results from three tables separately, with multiple correction within each table. The three tables of interest in such a setup in my opinion are, e.g. the condition effect for group 1, for group 2, and the interaction (whether the condition effect is different). Since I know ahead of time that I will present all three tables, I'm not doing any changes to the analysis based on what results I obtain, and I correct for multiple testing within each table, so the reported FDRs can be interpreted appropriately.

ADD COMMENTlink written 9 days ago by Michael Love18k
Please log in to add an answer.

Help
Access

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