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

Fix downsampling logic for QuSAGE

parent 0edebb0a
No related merge requests found
...@@ -473,7 +473,7 @@ scQuSAGE <- function(sc10x,gs,res.use=0.1,ds=25000,nm="Lin",folder="lin"){ ...@@ -473,7 +473,7 @@ scQuSAGE <- function(sc10x,gs,res.use=0.1,ds=25000,nm="Lin",folder="lin"){
number.clusters <- length(unique(sc10x@ident)) number.clusters <- length(unique(sc10x@ident))
labels <- paste0("Cluster_",as.vector(factor(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) rnd <- sample(1:ncol(sc10x@data),ds)
data <- sc10x@data[,rnd] data <- sc10x@data[,rnd]
labels <- labels[rnd] labels <- labels[rnd]
...@@ -618,7 +618,7 @@ scQuSAGEsm <- function(sc10x,gs,ds=25000,nm="Lin",folder="lin"){ ...@@ -618,7 +618,7 @@ scQuSAGEsm <- function(sc10x,gs,ds=25000,nm="Lin",folder="lin"){
clusters <- unique(sc10x@ident) clusters <- unique(sc10x@ident)
labels <- as.vector(factor(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) rnd <- sample(1:ncol(sc10x@data),ds)
data <- sc10x@data[,rnd] data <- sc10x@data[,rnd]
labels <- labels[rnd] labels <- labels[rnd]
...@@ -722,7 +722,7 @@ scQuSAGElg <- function(sc10x,gs,ds=25000,nm="Lin",folder="lin"){ ...@@ -722,7 +722,7 @@ scQuSAGElg <- function(sc10x,gs,ds=25000,nm="Lin",folder="lin"){
clusters <- unique(sc10x@ident) clusters <- unique(sc10x@ident)
labels <- as.vector(factor(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) rnd <- sample(1:ncol(sc10x@data),ds)
data <- sc10x@data[,rnd] data <- sc10x@data[,rnd]
labels <- labels[rnd] labels <- labels[rnd]
......
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