I have RNA-seq data that I would like to visualize via a GOterm heat map between time points. I have successfully been able to run DESeq2 on my transcriptome, but the output file with logpvalue/pvalue is my entire transcriptome, not compartmentalized by time point (0, 6, 12, 24, 48hr). I was able to make a heat map from this via time point but again, I have A LOT of genes so wanted to group by GOterm.
I have examined GO enrichment with programs like GOMWU, but again this examines the entire transcriptome, not each time point to compare between samples. Is there a way to use DESeq2 results and merge with GO terms to create a heat map?
I'm hoping my wording for this questions isn't too disorganized. Any suggestions on packages or programs that might help with my goal? I'm fairly new to all of this so apologies in advance.