diff --git a/r.scripts/sc-TissueMapper.R b/r.scripts/sc-TissueMapper.R index ecc5c2ba400f179ba4535013f0d9baf67f15e481..d28716294af4f36352d03160c0af39b2c272b747 100644 --- a/r.scripts/sc-TissueMapper.R +++ b/r.scripts/sc-TissueMapper.R @@ -473,7 +473,7 @@ scQuSAGE <- function(sc10x,gs,res.use=0.1,ds=25000,nm="Lin",folder="lin"){ number.clusters <- length(unique(sc10x@ident)) labels <- paste0("Cluster_",as.vector(factor(sc10x@ident))) - if ((ncol(sc10x@data)<ds & ds!=0)){ + if ((ncol(sc10x@data)>ds & ds!=0)){ rnd <- sample(1:ncol(sc10x@data),ds) data <- sc10x@data[,rnd] labels <- labels[rnd] @@ -618,7 +618,7 @@ scQuSAGEsm <- function(sc10x,gs,ds=25000,nm="Lin",folder="lin"){ clusters <- unique(sc10x@ident) labels <- as.vector(factor(sc10x@ident)) - if ((ncol(sc10x@data)<ds & ds!=0)){ + if ((ncol(sc10x@data)>ds & ds!=0)){ rnd <- sample(1:ncol(sc10x@data),ds) data <- sc10x@data[,rnd] labels <- labels[rnd] @@ -722,7 +722,7 @@ scQuSAGElg <- function(sc10x,gs,ds=25000,nm="Lin",folder="lin"){ clusters <- unique(sc10x@ident) labels <- as.vector(factor(sc10x@ident)) - if ((ncol(sc10x@data)<ds & ds!=0)){ + if ((ncol(sc10x@data)>ds & ds!=0)){ rnd <- sample(1:ncol(sc10x@data),ds) data <- sc10x@data[,rnd] labels <- labels[rnd]