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