Re: [Bioc-devel] SummarizedExperiments not equal after serialisation

2019-05-16 Thread Kasper Daniel Hansen
Interesting detective work. This is nasty. Best, Kasper On Thu, May 16, 2019 at 2:19 AM Pages, Herve wrote: > Let's try to go to the bottom of this. But let's leave > SummarizedExperiment objects out of the picture for now and focus on what > happens with a very simple reference object. > > Whe

Re: [Bioc-devel] SummarizedExperiments not equal after serialisation

2019-05-15 Thread Pages, Herve
Let's try to go to the bottom of this. But let's leave SummarizedExperiment objects out of the picture for now and focus on what happens with a very simple reference object. When you create 2 instances of a reference class with the same content: A <- setRefClass("A", fields=c(stuff="ANY"))

Re: [Bioc-devel] SummarizedExperiments not equal after serialisation

2019-05-11 Thread Aaron Lun
I would say it's much worse than mismatching class definitions. https://github.com/Bioconductor/SummarizedExperiment/issues/16 -A On 5/11/19 5:07 AM, Martin Morgan wrote: I think it has to do with the use of reference classes in the assay slot, which have different environments se = Summa

Re: [Bioc-devel] SummarizedExperiments not equal after serialisation

2019-05-11 Thread Martin Morgan
I think it has to do with the use of reference classes in the assay slot, which have different environments se = SummarizedExperiment() saveRDS(se, fl <- tempfile()) se1 = readRDS(fl) and then > all.equal(se@assays, se1@assays) [1] "Class definitions are not identical" > all.equal(se@assa

[Bioc-devel] SummarizedExperiments not equal after serialisation

2019-05-11 Thread Laurent Gatto
I would appreciate some background about the following: > suppressPackageStartupMessages(library("SummarizedExperiment")) > set.seed(1L) > m <- matrix(rnorm(16), ncol = 4, dimnames = list(letters[1:4], LETTERS[1:4])) > rowdata <- DataFrame(X = 1:4, row.names = letters[1:4]) > se1 <- SummarizedExpe