I am quite new to R, but using the DESeq2 package for analyzing total RNA-seq. data.
I would like to do a heatmap of only a subset of the differentially expressed genes, i.e. only protein-coding genes.
Thus, how can I inpute information on gene_biotype into the S4 object (results of DESeq2) or directly into the dds?
The only solution I have found involve converting res to a dataframe and merging with a dataframe containing the info on gene_biotype, but as this does not subset the s4 object data I am pretty stuck.
I would very much appreciate any suggestions as intensive google-ing and reading up on various posts over many days have not yet brought me closer to a solution.