This is a test version of Biostars. For the public version, visit https://www.biostars.org.
spike-in Single cell RNA-seq

Hi all, I need to define control features (genes) - ERCC spike-ins and mitochondrial genes when I do single cell RNA-seq analysis. while isSpike in SingleCellExperiment package defunct now.

isSpike(umi, "ERCC") <- grepl("^ERCC-", rownames(umi))
isSpike(umi, "MT") <- rownames(umi) %in% 
c("ENSG00000198899", "ENSG00000198727", "ENSG00000198888",
"ENSG00000198886", "ENSG00000212907", "ENSG00000198786",
"ENSG00000198695", "ENSG00000198712", "ENSG00000198804",
"ENSG00000198763", "ENSG00000228253", "ENSG00000198938",
"ENSG00000198840")

Did anyone know how I can do this step in R now? Thanks in advance!

r sequencing rna-seq

0 answers

No answers yet.

Log in to answer this question.