CummeRbund is driving me crazy. Apparently, Cuffdiff labels gene_ids according to the GTF file provided. My genes are all labelled with the Ensemble ids. I am not willing to rerun cuffdiff because my bam files are so large that it took more than 24hrs to run on 48 processors.
I wrote a script that annotates my cuffdiff output with the gene common names and descriptions for the >3500 significantly differentially expressed genes, so I have that problem solved.
but now my boss wants some expression level bar graphs from CummeRbund for approximately 20 genes (a simple request). I have spent a whole day trying to change the labels on the graphs from "ENSMUST.....etc" to a more understandable gene name, and have succumb to taking screenshots of the graphs and changing the labels on paint. This is ridiculous, does anyone know a better way to change labels?
I wrote a script that annotates my cuffdiff output with the gene common names and descriptions for the >3500 significantly differentially expressed genes, so I have that problem solved.
but now my boss wants some expression level bar graphs from CummeRbund for approximately 20 genes (a simple request). I have spent a whole day trying to change the labels on the graphs from "ENSMUST.....etc" to a more understandable gene name, and have succumb to taking screenshots of the graphs and changing the labels on paint. This is ridiculous, does anyone know a better way to change labels?
Comment