[Bioc-devel] MultiAssayExperiment silently converting RangedSummarizedExperiments into SummarizedExperiments.
Charles Plessy
ch@r|e@@p|e@@y @end|ng |rom o|@t@jp
Wed Aug 28 07:28:55 CEST 2019
Hello,
I am using MultiAssayExperiment as a base class in the CAGEr package.
In Bioconductor 3.10, CAGEr is broken because updates of
RangeSummarizedExperiment experiments trigger their silent conversion
to SummarizedExperiment objects (basically discarding the ranges). This
only happens when there is another experiment slot that contains a plain
SummarizedExperiment.
Here is a toy example to reproduce the issue:
```{r}
library("MultiAssayExperiment")
(DF <- DataFrame(sample1=c(1,2), sample2=c(1,3)))
(rRanges <- GPos("chr1", 1:2, "+"))
(cData <- DataFrame(type=c("control", "treatment"), row.names =
c("sample1", "sample2")))
(rSE <- SummarizedExperiment(SimpleList(DF), rowRanges = rRanges,
colData = cData))
(SE <- SummarizedExperiment(SimpleList(DF), colData = cData))
(m <- MultiAssayExperiment(ExperimentList(exp1=rSE)))
(m[[1]] <- m[[1]])
m # No change
(m <- c(m, exp2=SE))
(m[[1]] <- m[[1]])
m # RangedSummarizedExperiment replaced with SummarizedExperiment
sessionInfo()
```
Am I misusing the MultiAssayExperiment objects, or is this a bug ?
Have a nice day,
Charles
--
Charles Plessy - - ~ ~ ~ ~ ~ ~~~~ ~ ~ ~ ~ ~ - - charles.plessy using oist.jp
Okinawa Institute of Science and Technology Graduate University
Staff scientist in the Luscombe Unit - ~ - https://groups.oist.jp/grsu
Toots from work - ~ ~~ ~ - https://mastodon.technology/@charles_plessy
More information about the Bioc-devel
mailing list