[BioC] IRanges problem: findOverlaps

Tim Triche, Jr. tim.triche at gmail.com
Tue Jul 10 18:46:11 CEST 2012


how about

resize(genes,1, fix='start')
findOverlaps(peaks, genes)

however,

+1 for Julie Zhu's ChIPpeakAnno package since, as was just pointed
out, it is meant for such things!




On Tue, Jul 10, 2012 at 9:43 AM, Jonathan Cairns
<Jonathan.Cairns at cancer.org.uk> wrote:
> Hi,
>
> >From ?findOverlaps:
>
> "If ‘type’ is ‘start’ or ‘end’, the intervals are required to have matching starts or ends, respectively."
>
> This doesn't sound like it's what you're after, as type = "start" considers only the starts of peaks, as well as the starts of whatever annotation you're using.
>
> Have you looked at the package ChIPpeakAnno? It is designed for annotating ChIP-seq peaks and might be useful.
>
> Jonathan
>
> ________________________________________
> From: bioconductor-bounces at r-project.org [bioconductor-bounces at r-project.org] On Behalf Of Nicolas DESCOSTES [descostes at ciml.univ-mrs.fr]
> Sent: 10 July 2012 16:40
> To: bioconductor at r-project.org
> Subject: [BioC] IRanges problem: findOverlaps
>
> Dear members,
>
> I have a set of chip-Seq peaks which I want to find the overlap with a bench of annotations downloaded from ucsc.
>
> When I am doing the overlap between my peaks and annotations list, I am finding 5635 positive matches. When adding "type = "start"", I have no results returned. However, when I am visualizing the 5635 intervals, it overlaps many start site of my genes annotations.
>
> I tried to update IRanges but I am still getting no results.
>
> Any idea?
>
> Thanks.
>
>
>         [[alternative HTML version deleted]]
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
>
> NOTICE AND DISCLAIMER
> This e-mail (including any attachments) is intended for ...{{dropped:16}}
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor



-- 
A model is a lie that helps you see the truth.

Howard Skipper



More information about the Bioconductor mailing list