I am following this tutorial for DEG analysis of my RNAseq data. However, my data has one factor with three different levels.
How can I create a heatmap including all the DEGs from the pairwise comparisons that DESeq2 has outputted? I have three DESeq2 results files (one for each pairwise comparison between my three groups).
The usual way is to plot each pairwise result individually. You’ll see this in some publications.
For something more interesting, maybe try the other DE tools such as Limma – you can set one factor level as the contrast to see the effect on the other two conditions (factor levels), all together one set of visualization graphs. Using the control condition as the contrast is one place to start.