Circumventing statistical non-significance on a triple interaction using edgeR
1
0
Entering edit mode
David R ▴ 90
@david-rengel-6321
Last seen 10 months ago
European Union

Hi,

I have already discussed in the past this experimental design in which I deal with three factors (inoc: Bact or Mock; tissue: IN or OUT and hpi: 0,1,2,4 and 6h). Indeed, I am interested in the triple interaction, written on edgeR as:

 ((Bact.IN.Tx - MOCK.IN.Tx) – (Bact.IN.T0 - MOCK.IN.T0)) –

((Bact.OUT.Tx- MOCK.OUT.Tx) - (Bact.OUT.T0- MOCK.OUT.T0))

Now, I am aware that when all fitted counts on one side of the contrast are 0, then the contrast itself will be non-significant because of the indetermination of the zero on zero likelihood. I can circumvent this in the case of a double interaction. Let’s say I am looking at the double interaction between inoc and time. A given gene has 0 counts at T0, whatever Bact or Mock, but not at Tx. Hence, the double interaction is non-significant. However, I may say that, biologically, it makes sense to consider this gene in interaction if there is not inoc effect at T0 (obvious in this case) AND there is a significant inoc effect at Tx. Does this make sense to you?

Well, if that approach with the double sounded reasonable, when dealing with the triple interaction this may get trickier. The question would be: which contrasts should be considered to be taken separately at T0 and Tx when the statistical triple interaction gives non-significance because of the all-zeros at T0? The double interactions at Tx maybe?

Thanks a lot for your help.

David

rnaseq edgeR makecontrasts interactions • 1.1k views
ADD COMMENT
0
Entering edit mode
Aaron Lun ★ 28k
@alun
Last seen 8 hours ago
The city by the bay

However, I may say that, biologically, it makes sense to consider this gene in interaction if there is not inoc effect at T0 (obvious in this case) AND there is a significant inoc effect at Tx. Does this make sense to you?

No, it doesn't. The whole point of this particular phenomenon is that with two zeroes, you have no information to determine whether a gene is DE or not. You can't even say that there is no inoculation effect at T0. For example, what if T0 is sequenced at lower depth than Tx? Or if the gene is expressed at much lower levels in T0 compared to Tx? For all we know, the gene has the same true log-fold change upon inoculation in T0 and Tx, i.e., no interaction - we just don't have deep enough coverage to know.

In this situation, the lack of significance from edgeR is the correct behaviour, as there's not enough evidence to say otherwise. In contrast, your approach would happily consider a gene as exhibiting an interaction effect, even when the interaction is truly zero. If you want to identify genes that are significantly DE upon inoculation at Tx and not at T0, you can do so, e.g., by intersecting the sets of DE and non-DE genes - but don't call it an interaction effect.

ADD COMMENT
0
Entering edit mode

Thanks Aaron for your remarks, with which I completely agree. I understand the "depth" problem etc, i.e. the eventual lack of information at T0. However, it is also true that, in fact, many genes are not expressed at T0 and they are only expressed upon bacterial inoculation, whose effects are seen at Tx. My approach was intended to put those genes back in. And I agree that the term "interaction" is tricky here. However, I used it because I was thinking of using significant double interactions at Tx between inoc and tissue in order to fish out those genes that are not significant for the triple interaction.

Thanks again.,

D.

 

ADD REPLY

Login before adding your answer.

Traffic: 1025 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6