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

Remove low end cutoff for cell with low gene content

parent ec340b65
Branches
Tags
2 merge requests!4Merge ds,!3Merge completed D17 downsample in to develop
...@@ -94,7 +94,7 @@ if (!dir.exists("./analysis/pseudotime")){ ...@@ -94,7 +94,7 @@ if (!dir.exists("./analysis/pseudotime")){
option_list=list( option_list=list(
make_option("--p",action="store",default="DPrF",type='character',help="Project Name"), make_option("--p",action="store",default="DPrF",type='character',help="Project Name"),
make_option("--g",action="store",default="ALL",type='character',help="Group To analyze"), make_option("--g",action="store",default="ALL",type='character',help="Group To analyze"),
make_option("--lg",action="store",default=500,type='integer',help="Threshold for cells with minimum genes"), make_option("--lg",action="store",default=0,type='integer',help="Threshold for cells with minimum genes"),
make_option("--hg",action="store",default=3000,type='integer',help="Threshold for cells with maximum genes"), make_option("--hg",action="store",default=3000,type='integer',help="Threshold for cells with maximum genes"),
make_option("--lm",action="store",default=0,type='numeric',help="Threshold for cells with minimum %mito genes"), make_option("--lm",action="store",default=0,type='numeric',help="Threshold for cells with minimum %mito genes"),
make_option("--hm",action="store",default=0.1,type='numeric',help="Threshold for cells with maximum %mito genes"), make_option("--hm",action="store",default=0.1,type='numeric',help="Threshold for cells with maximum %mito genes"),
...@@ -114,6 +114,7 @@ option_list=list( ...@@ -114,6 +114,7 @@ option_list=list(
) )
opt=parse_args(OptionParser(option_list=option_list)) opt=parse_args(OptionParser(option_list=option_list))
rm(option_list) rm(option_list)
if (opt$lg==0){opt$lg=-Inf}
if (opt$lm==0){opt$lm=-Inf} if (opt$lm==0){opt$lm=-Inf}
sc10x.data <- Read10X(data.dir="./analysis/DATA/10x/filtered_gene_bc_matrices/GRCh38/") sc10x.data <- Read10X(data.dir="./analysis/DATA/10x/filtered_gene_bc_matrices/GRCh38/")
......
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