Question: diffbind dba.plotvolcano bLabels option use genames instead of numbers
0
4 months ago by
Germany

Hi all, So am using the diffbind R package and I was doing some plots after analyzing my samples. using the dba.plotVolcano and the option

bLabels = TRUE

I can have the location of those 50 genes that are on the top of the list.

That is great but is not that informative

Is it possible instead of the numbers to plot the gene names?

P.S. is it possible to import the table to another package and plot the volcano plotting function, as described in the following link? https://www.gettinggeneticsdone.com/2016/01/repel-overlapping-text-labels-in-ggplot2.html what do I need to transform or change?

diffbind • 138 views
modified 12 weeks ago by Rory Stark2.9k • written 4 months ago by theodore.georgomanolis0
0
12 weeks ago by
Rory Stark2.9k
CRUK, Cambridge, UK
Rory Stark2.9k wrote:

DiffBind currently does not support adding labels such as gene names to the differentially bound sites.

As you suggest, the best bet is to export the table, and annotate and plot it using other tools. You can get the full table either using dba.report() with th=1, or use only the differentially bound sites by assigning the return value from the call to dba.plotVolcano(). The relevant fields are the FDR and Fold columns in the table.

Thank you Roy for the answer, do you plan to implement an option as such?

This is on my list of possible features to add, but there are a number of things ahead of it! The main thing it requires is an interface to add an annotation to the peakset, as most ChIP peak intervals do not correspond directly to annotated genes.