Question: Exclude interaction terms from the design
0
gravatar for polyptwo
29 days ago by
polyptwo0
polyptwo0 wrote:

Good afternoon,

I am aware that this questions is not only relevant to DESeq2 but is more general about Regression Models, but my application right now is in DESeq2.

I have design similar to the examples: 2 genotypes, 2 conditions. We wondered if the genotypes were behaving differently to the condition, so we modelled with an interaction term.

design = ~ condition + genotype + condition:genotype

Turns out that there is no difference, and plotting some transcripts confirms it: they behave the same way after treatment. Knowing that they behave the same way is important biological information for our research that we will use.

At this point, we are therefore interested in knowing the difference between the genotypes. However, is there any benefit in keeping the interaction term in the design? Can we exclude the interaction term from the design while maintaining the condition factor and treating the samples as biological replicates using a simplified design?

Alternatively, is there a way of getting the overall effect between genotypes as currently modelled? genotype_B_vs_A or contrast=c("genotype","B","A")) is the effect between genotypes untreated, and list(c("genotype_B_vs_A","genotypeB.conditionTrt")) is the effect between genotypes treated.

Is then the overall difference between both genotypes across all treatments accessible? If so, should it give the same result as genotype_B_vs_A in the following design without interaction?

design = ~ condition + genotype

Thank you very much for your support, and advice.

deseq2 statistics modelling • 90 views
ADD COMMENTlink modified 28 days ago by Michael Love26k • written 29 days ago by polyptwo0
Answer: Exclude interaction terms from the design
0
gravatar for Michael Love
28 days ago by
Michael Love26k
United States
Michael Love26k wrote:

If you are pretty sure that the effect is consistent across genotype, I'd go with refitting the model with your second design. I don't worry too much about the extra exploration of the data (e.g. garden of forking paths), given the context of high-dimensional analysis.

ADD COMMENTlink written 28 days ago by Michael Love26k
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 16.09
Traffic: 443 users visited in the last hour