diff --git a/R/bio/cummerutils.R b/R/bio/cummerutils.R
index ed55dd22647e5ac6b5ee1f13e0289a8370742b88..6522de5dbe55bd6887524d797ae99a66532551a8 100644
--- a/R/bio/cummerutils.R
+++ b/R/bio/cummerutils.R
@@ -1,14 +1,18 @@
 
-require.auto(cummeRbund)
+require(cummeRbund)
 
 
-
-getExpressedGenes <- function(cuff, minFPKM=1, logMode=F){
+getExpressedGenes <- function(cuff, ...){
     fpkmMat<-repFpkmMatrix(genes(cuff))
 
+    rownames(filterByExpression(fpkmMat, ...))
+}
+
+
+filterByExpression <- function(fpkmMat, minFPKM=1, logMode=F){
     if(logMode) fpkmMat<-log10(fpkmMat+1) ## add a pseudocount
 
     geneMax <- apply(fpkmMat, 1, max)
 
-    rownames(fpkmMat[geneMax>minFPKM,])
+    fpkmMat[geneMax>minFPKM,]
 }