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. Did anyone know how I can do this step in R now? Thanks in advance! Yachen
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")
• 86 views
•
link
0 answers
No answers yet.
Log in to answer this question.