Splitting single cell data
0
0
Entering edit mode
8 days ago
otieno43 ▴ 40

I merged single cell RNA seq data but now I want to separate them. how do I do this? Below is what I did..

sce.rep1 <- read10xCounts(fname, col.names=TRUE)
sce.rep2 <- read10xCounts(fname1, col.names=TRUE)

I the combined the two

All_scell <- list(sce.rep1, sce.rep2)

I then order the data

All_scell <- mget(ls(pattern = "^sce.rep[0-9]"))

Removed rRNA

rRNA_Scell <- import.gff3("/Users/TrypDB.gff") %>% 
as_tibble() %>%
filter(type=="rRNA") %>%
select(ID, description) %>%
mutate(ID = gsub("\\.1", "", ID))
for (i in 1:length(All_scell)) {
All_scell[[i]] <- All_scell[[i]][!rownames(All_scell[[i]])%in%rRNA_Scell$ID,]}

Removed empty cells

set.seed(123456)
All_scell <- lapply(All_scell, function(x){
    e.out <- emptyDrops(counts(x), lower = 1000)
    x <- x[ ,which(e.out$FDR<=0.001)]})

Now I want to separate the filtered two datasets (sce.rep1 , sce.rep2) as a seurat object. some help please, thanks.

ScRNA • 300 views
ADD COMMENT

Login before adding your answer.

Traffic: 1330 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6