Re: DiffBind annotation

Rory Stark <[email protected]> Fri, 12 Sep 2014 14:20:19 +0000
Newsgroups gmane.science.biology.informatics.conductor
Message-ID <D038BDE6.49CCE%[email protected]>
Hi Igor-

We avoided adding annotation to DiffBind because of the other packages
that specialize in this. They don't take a DiffBind object directly, but
you can always retrieve any peak set as a GRanges object, which is
commonly used in Bioconductor packages. The main packages I use for
annotation and gene set enrichment are ChIPpeakAnno and chipenrich.

Peaksets are returned as GRanges objects by default by dba.peakset() (with
bRetrieve=TRUE), dba.report() (for differentially bound peaks), and
dba.overlap() (for peaks that do or don't overlap between peak sets; you
can also use dbaPlotVenn() with bReturnPeaksets=TRUE).

Cheers-
Rory

On 11/09/2014 21:22, "Dolgalev, Igor" <[email protected]> wrote:

>Hello.
>
>I am getting started with DiffBind. I am really enjoying all the options
>it provides. Unfortunately, it does not seem to do any annotation. For
>example, it would be great to get the closest gene or distance to TSS/TTS
>for the binding sites. I know there are tools available for this purpose.
>Is there a specific one that you would recommend and that would integrate
>with DiffBind (if it could take DBA object directly, for example)?
>
>Any assistance would be much appreciated.
>
>Thank you very much.
>
>Igor
>

_______________________________________________
Bioconductor mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/bioconductor
Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor