I’m using CummeRbound for analyzing Cuffdiff output. I need to focus the analysis of expression dynamics to a restricted gene set (roughly 50-100 genes). My aim would be to construct a heat map where is plotted the expression level of selected genes across the samples.
From the manual they provided a function to "grab a single gene of interest":
myGeneId<-"XLOC_020405" myGene<-getGene(cuff,myGeneId) myGene
How can I extend this analysis from one single gene to a set of genes?
Thank you in advance.