I have 4 samples (Adar1-null stimulated, Adar1-null unstimulated, Control Stimulated, and Control Unstimulated). For each of these 4 samples, there are 3 replicates.
I have performed the differential gene expression analysis by doing pairwise comparisons and I am now trying to produce a heatmap for the differentially expressed genes for Adar1-null stimulated vs. Control stimulated.
I have stored the result of the comparison in a variable, but I don't know how to use this variable to generate a heatmap. The only heatmap I was able to generate was the one across all samples, but I am only interested in generating the heatmap across these two samples. How do I do this using DESeq2 and R?