I am using R version 3.4.1, Bioconductor 3.5 and RStudio to use scater package. I have 7 single-cell RNA seq datasets (from a 10x Genomics experiment), each corresponding to a different time point. Each dataset is a SCESet object. I want to merge the data all together. I have found there is the
mergeSCESet function. However, this works perfectly well on some pairs of datasets, but not all. If not, I got this error :
Error in Biobase::combine(exprs(x), exprs(y)) : matrix shared row and column elements differ: Mean relative difference: 1.102234
What is this mean relative difference ? How could I fix this or bypass it ?
Here is an example code of what my program is doing (there is no error when running this code - taken from https://rdrr.io/bioc/scater/man/mergeSCESet.html ) :
library(scater) data("sc_example_counts") data("sc_example_cell_info") pd <- new("AnnotatedDataFrame", data = sc_example_cell_info) example_sceset <- newSCESet(countData = sc_example_counts, phenoData = pd) a <- mergeSCESet(example_sceset[, 1:20], example_sceset[, 21:30])