Skip to content
Snippets Groups Projects
Commit b5c34682 authored by Gervaise Henry's avatar Gervaise Henry :cowboy:
Browse files

Fix sample name error in sc_Demultiplex.R

parent ccbb5d14
No related merge requests found
......@@ -6,7 +6,7 @@ library(readr)
#retrive command line options
option_list=list(
make_option("--p",action="store",default="Pr",type='character',help="Project Name"),
make_option("--p",action="store",default="ALL",type='character',help="Project Name"),
make_option("--d",action="store",default=4,type='integer',help="Demultiplex Group Number"),
make_option("--mc",action="store",default=3,type='integer',help="Minimum Cells"),
make_option("--mg",action="store",default=200,type='integer',help="Minimum Genes")
......@@ -24,7 +24,7 @@ cell.codes <- as.data.frame(sc10x@raw.data@Dimnames[[2]])
colnames(cell.codes) <- "barcodes"
rownames(cell.codes) <- cell.codes$barcodes
cell.codes$lib.codes <- as.factor(gsub(pattern=".+-",replacement="",cell.codes$barcodes))
cell.codes$samples <- as.vector(sc10x.aggr$library_id[cell.codes$lib.codes])
cell.codes$samples <- sc10x.aggr$library_id[match(cell.codes$lib.codes,as.numeric(rownames(sc10x.aggr)))]
sc10x <- CreateSeuratObject(raw.data=sc10x.data,meta.data=cell.codes["samples"],min.cells=opt$mc,min.genes=opt$mg,project=Project.Name)
rm(sc10x.data)
rm(cell.codes)
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment