does this contradict the doc in some way?  the second range of your result
has width zero, which i think is correct.

On Sun, Nov 11, 2012 at 12:16 PM, Kasper Daniel Hansen <
kasperdanielhansen@gmail.com> wrote:

> Example
>
> library(GenomicRanges)
> gr1 = GRanges(seqnames = c("chr1", "chr1"), ranges = IRanges(start =
> c(1,5), width = 7))
> gr2 = GRanges(seqnames = "chr1", ranges = IRanges(start = 1, end = 4))
> pintersect(gr1, rep(gr2, 2))
>
> GRanges with 2 ranges and 0 metadata columns:
>       seqnames    ranges strand
>          <Rle> <IRanges>  <Rle>
>   [1]     chr1    [1, 4]      *
>   [2]     chr1    [5, 4]      *
>   ---
>   seqlengths:
>    chr1
>      NA
>
> The second element in the return GRanges should be empty.  I would
> therefore expect that the return length of the pintersect would have
> length 1, despite the fact that the first object has length 2.
>
> Kasper
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor@r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives:
> http://news.gmane.org/gmane.science.biology.informatics.conductor
>

	[[alternative HTML version deleted]]

