[Bioc-devel] Bug in setdiff() for signature 'IRangesList, IRangesList'?

[Michael Lawrence]

I fixed this in 2.9.6 and 2.8.1.

It brings to light an interesting pecularity of the recycling rule.

Thanks,
Michael

On Mon, Nov 7, 2016 at 9:46 AM, Leonard Goldstein
<goldstein.leonard at gene.com> wrote:
> Dear Hervé et al.
>
> I noticed that setdiff(x, y) where 'x' is an IRangesList containing an
> empty IRanges results in an error (see below). I would have expected the
> function to return the IRangesList with an empty IRanges?
>
> Many thanks for your help.
>
> Leonard
>
> --
>> x <- IRangesList(IRanges())
>> y <- IRangesList(IRanges())
>>
>> setdiff(x, y)
> Error in ans[] <- x : replacement has length zero
>>
>> sessionInfo()
> R version 3.3.0 Patched (2016-05-09 r70594)
> Platform: x86_64-pc-linux-gnu (64-bit)
> Running under: Red Hat Enterprise Linux Server release 6.6 (Santiago)
>
> locale:
>  [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C
>  [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8
>  [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8
>  [7] LC_PAPER=en_US.UTF-8       LC_NAME=C
>  [9] LC_ADDRESS=C               LC_TELEPHONE=C
> [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
>
> attached base packages:
> [1] stats4    parallel  stats     graphics  grDevices utils     datasets
> [8] methods   base
>
> other attached packages:
> [1] IRanges_2.8.0       S4Vectors_0.12.0    BiocGenerics_0.20.0
>
> loaded via a namespace (and not attached):
> [1] zlibbioc_1.20.0 XVector_0.14.0
>
>         [[alternative HTML version deleted]]
>
> _______________________________________________
> Bioc-devel at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/bioc-devel