Question: spike-in Single cell RNA-seq
0
gravatar for mikysyc2016
2 days ago by
mikysyc201670
mikysyc201670 wrote:

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!

sequencing rna-seq R • 64 views
ADD COMMENTlink modified 2 days ago • written 2 days ago by mikysyc201670
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1717 users visited in the last hour