Browse code

Bugfix when using ssgsea with one single sample in the input expression data matrix. Reported by Christopher Eeles at https://github.com/rcastelo/GSVA/issues/61

Robert Castelo authored on 31/05/2022 07:32:44
Showing2 changed files

... ...
@@ -1,5 +1,5 @@
1 1
 Package: GSVA
2
-Version: 1.45.0
2
+Version: 1.45.1
3 3
 Title: Gene Set Variation Analysis for microarray and RNA-seq data
4 4
 Authors@R: c(person("Justin", "Guinney", role=c("aut", "cre"), email="justin.guinney@sagebase.org"),
5 5
              person("Robert", "Castelo", role="aut", email="robert.castelo@upf.edu"),
... ...
@@ -875,8 +875,7 @@ ssgsea <- function(X, geneSets, alpha=0.25, parallel.sz,
875 875
     print("Normalizing...")
876 876
     ## normalize enrichment scores by using the entire data set, as indicated
877 877
     ## by Barbie et al., 2009, online methods, pg. 2
878
-    # es <- apply(es, 2, function(x, es) x / (range(es)[2] - range(es)[1]), es)
879
-    es <- es[,1:n] / (range(es)[2] - range(es)[1])
878
+    es <- es[, 1:n, drop=FALSE] / (range(es)[2] - range(es)[1])
880 879
   }
881 880
   
882 881
   if (length(geneSets) == 1)