Commit c8a38d4a authored by Gervaise Henry's avatar Gervaise Henry 🤠
Browse files

Split cutoff figs by sample names

parent e2689018
......@@ -703,7 +703,7 @@ scScore <- function(sc10x.l,score,geneset,cut.pt=0.9,anchor=FALSE){
sc10x.temp <- subset(sc10x,cells=cells.remove,invert=TRUE)
thresh <- list()
thresh[["all"]] <- scThresh(list(all=sc10x.temp),feature=paste0(score,"1"),sub=score)
thresh[["all"]] <- scThresh(list(all=sc10x.temp),feature=paste0(score,"1"),sub=paste0(score,".",i))
cut.x <- thresh$all$all$threshold[thresh$all$all$method=="RenyiEntropy"]
} else if (cut.pt == "triangle"){
h <- hist(data.frame(sc10x[[paste0(score,"1")]])[,paste0(score,"1")],breaks=1000,plot=FALSE)
......@@ -712,7 +712,7 @@ scScore <- function(sc10x.l,score,geneset,cut.pt=0.9,anchor=FALSE){
sc10x.temp <- subset(sc10x,cells=cells.remove,invert=TRUE)
thresh <- list()
thresh[["all"]] <- scThresh(list(all=sc10x.temp),feature=paste0(score,"1"),sub=score)
thresh[["all"]] <- scThresh(list(all=sc10x.temp),feature=paste0(score,"1"),sub=paste0(score,".",i))
cut.x <- thresh$all$all$threshold[thresh$all$all$method=="Triangle"]
} else if (cut.pt == "minerror"){
h <- hist(data.frame(sc10x[[paste0(score,"1")]])[,paste0(score,"1")],breaks=1000,plot=FALSE)
......@@ -721,7 +721,7 @@ scScore <- function(sc10x.l,score,geneset,cut.pt=0.9,anchor=FALSE){
sc10x.temp <- subset(sc10x,cells=cells.remove,invert=TRUE)
thresh <- list()
thresh[["all"]] <- scThresh(list(all=sc10x.temp),feature=paste0(score,"1"),sub=score)
thresh[["all"]] <- scThresh(list(all=sc10x.temp),feature=paste0(score,"1"),sub=paste0(score,".",i))
cut.x <- thresh$all$all$threshold[thresh$all$all$method=="MinErrorI"]
} else {
cut.x <- quantile(cdf,probs=cut.pt)
......
Markdown is supported
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