cummeRbund getGenes gene order

[Enrique Zudaire]

Hi all,

I'm using cummeRbund expressionBarplot to plot genes that I retrieve using getGenes

myGeneIds<-c('Atg12', 'Atg16l1','Atg4a')
myGenes<-getGenes(cuff, myGeneIds)
expressionBarplot(myHematopoieticGenes, showErrorbars=TRUE)

The bar plot shows the genes but not in the order I specified in the string vector but in ascending order of gene Ids (ENSMUSG000000XXXXX). Is there a way to keep the same order I had in the string vector?

Thanks in advance for the help
K

[millerma1]

Did you ever figure it out? I'm having the same problem with no idea how to fix it.

[Enrique Zudaire]

I kind of forgot about this, but yes, I believe I figured something using ggplot2.

I ended up writing a small function like this:

genesExpressionBarPlot<-expressionBarplot(

      genes, showErrorbars=myShowErrorBars, replicates=myReplicates)  +
      
      ...

      scale_x_discrete(limits=genesTable[,1], 
                        labels=genesTable[,3])
      ...
   }

I believe the trick is to have the gene_id from the cuff database (my genesTable[,1]) be "limits". The "labels" are optional but useful for further customization.

Let me know if it works, I haven't used this in a while.

[millerma1]

I just tried it but I'm obviously missing the object 'genesTable' what did you use for this?

Also, I thought that the reordering of the genes happened in the getGenes step. Did you bypass or change this in some way?

[Enrique Zudaire]

geneTable is a table I build with info I need from the cuff database and other sources. In the example above geneTable [,1] contains gene_id and geneTable[,3] contains gene labels.
In the example above genes is the result of getGenes() and goes before building the plot.

K