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