Merge pull request jtleek#4 from andrewejaffe/master

fixed bug with vfilter
DarwinAwardWinner · Feb 25, 2015 · fdcf505 · fdcf505
2 parents 21cee1c + 39deb94
commit fdcf505
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/R/sva.R b/R/sva.R
@@ -35,7 +35,7 @@ sva <- function(dat, mod, mod0 = NULL,n.sv=NULL,controls=NULL,method=c("irw","tw
       stop(paste("sva error: the number of genes used in the analysis must be between 100 and",dim(dat)[1],"\n"))
     }
     tmpv = rowVars(dat)
-    ind = which(rank(-tmpv) < vfilter)
+    ind = which(rank(-tmpv) <= vfilter)
     dat = dat[ind,]
   }