R/internal-nonzeroGenes.R

In roryk/bcbioRnaseq: RNA-Seq Utilities

.nonzeroGenes <- function(object) {
    raw <- counts(object, normalized = FALSE)
    nonzero <- rowSums(raw) > 0L
    genes <- rownames(raw[nonzero, , drop = FALSE])
    message(paste(length(genes), "non-zero genes detected"))
    sort(genes)
}