diff --git a/r.scripts/sc-TissueMapper_RUN.PdPbPc.Pb.R b/r.scripts/sc-TissueMapper_RUN.PdPbPc.Pb.R index abb3b1ea95253ae7c4b85371731c475dfd54b84b..e42948f9950bd51205a4c2d9656b8848c9fce924 100644 --- a/r.scripts/sc-TissueMapper_RUN.PdPbPc.Pb.R +++ b/r.scripts/sc-TissueMapper_RUN.PdPbPc.Pb.R @@ -107,7 +107,7 @@ option_list=list( make_option("--st",action="store",default=TRUE,type='logical',help="Remove stressed cells?"), make_option("--stg",action="store",default="dws",type='character',help="Geneset to use for stress ID"), make_option("--cut.stress",action="store",default=0.9,type='numeric',help="Cutoff for stress score"), - make_option("--res.poststress",action="store",default=0.5,type='numeric',help="Resolution to cluster, post-stress"), + make_option("--res.poststress",action="store",default=1,type='numeric',help="Resolution to cluster, post-stress"), make_option("--cut.ne",action="store",default=0.999,type='numeric',help="Cutoff for NE score") ) opt=parse_args(OptionParser(option_list=option_list)) @@ -116,7 +116,7 @@ if (opt$lm==0){opt$lm=-Inf} sc10x <- scLoad("Pb") -sc10x <- scSubset(sc10x,i="Tumor",g="Tumor") +sc10x <- scSubset(sc10x,i="Glandular",g="Glandular") if (opt$cc==TRUE){ results <- scCellCycle(sc10x)