Help with the output from mnnCorrect
0
0
Entering edit mode
@sharvari-gujja-6614
Last seen 5 weeks ago
United States

Hi,

I am trying to use mnnCorrect to identify shared cell population between two batches from different conditions - WT and KO.

classic.out <- batchelor::mnnCorrect(sce1, sce2, subset.row=chosen.hvgs , correct.all = TRUE)


Using HVGs as input to mnnCorrect function, I get:

 > classic.out
class: SingleCellExperiment
dim: 9608 2542
assays(1): corrected
rownames(9608): ENSMUSG00000000056.7 ENSMUSG00000000078.7 ... NA NA
Error in rownames<-(*tmp*, value = c("ENSMUSG00000000056.7", "ENSMUSG00000000078.7",  :
missing values not allowed in rownames


There are no NAs in the individual single cell objects. Can you please help with the error.

Thanks

mnnCorect MNN • 525 views
0
Entering edit mode

Do all of the values in chosen.hvgs appear in both rownames(sce1) and rownames(sce2)?

Which is to say, the output of this operation should be empty:

setdiff(chosen.hvgs, intersect(rownames(sce1), rownames(sce2)))


is it?

0
Entering edit mode

In the future, it would help to tag this question with the relevant package name, i.e., batchelor. Otherwise the maintainer (i.e., me) will not even see it.