Add no downsampling (NA) as an option for sc_PC.Score.NE.R

r.scripts/sc_PC.Score.NE.R

@@ -105,7 +105,11 @@ plot(plot)
 dev.off()
 
 #subsample all cells (+NE) to better visualize their clustering
-rnd <- sample(1:ncol(sc10x.Group@data),opt$sv)
+if (!is.na(opt$s)){
+  rnd <- sample(1:ncol(sc10x.Group@data),opt$sv)
+} else {
+  rnd <- 1:ncol(sc10x.Group@data)
+}
 sc10x.Group <- SetIdent(object=sc10x.Group,cells.use=sc10x.Group@data@Dimnames[[2]][rnd],ident.use="sample")
 if (clust$size[1]<clust$size[2]){
   clusterCut <- names(clusterCut[clusterCut==1])