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

Leonard Goldstein goldstein.leonard at gene.com
Mon Nov 7 20:37:42 CET 2016


Thanks Michael!

Leonard



On Mon, Nov 7, 2016 at 11:35 AM, Michael Lawrence <lawrence.michael at gene.com
> wrote:

> 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
>

	[[alternative HTML version deleted]]



More information about the Bioc-devel mailing list