From 6d4794a3863468a5094f57befb84bc3a68ba3875 Mon Sep 17 00:00:00 2001
From: "Gervaise H. Henry" <gervaise.henry@utsouthwestern.edu>
Date: Wed, 15 Aug 2018 16:21:07 -0500
Subject: [PATCH] Fix downsampling logic for QuSAGE
---
r.scripts/sc-TissueMapper.R | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/r.scripts/sc-TissueMapper.R b/r.scripts/sc-TissueMapper.R
index ecc5c2b..d287162 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]
--
GitLab