filter.genes: Filter gene expression based on gene counts

Description Usage Arguments Value Author(s) Examples

View source: R/metaseqr.filter.R

Description

This function performs the gene expression filtering based on gene read counts and a set of gene filter rules. For more details see the main help pages of metaseqr.

Usage

1
2
    filter.genes(gene.counts, gene.data, gene.filters,
    sample.list)

Arguments

gene.counts

a matrix of gene counts, preferably after the normalization procedure.

gene.data

an annotation data frame usually obtained with get.annotation containing the unique gene accession identifiers.

gene.filters

a named list with gene filters and their parameters. See the main help page of metaseqr for details.

sample.list

the list containing condition names and the samples under each condition.

Value

a named list with three members. The first member (result is a named list whose names are the gene filter names and its members are the filtered rownames of gene.data. The second member (cutoff is a named list whose names are the gene filter names and its members are the cutoff values corresponding to each filter. The third member is a matrix of binary flags (0 for non-filtered, 1 for filtered) for each gene. The rownames of the flag matrix correspond to gene ids.

Author(s)

Panagiotis Moulos

Examples

1
2
3
4
5
6
7
8
9
data("mm9.gene.data",package="metaseqR")
gene.counts <- mm9.gene.counts
sample.list <- sample.list.mm9
gene.counts <- normalize.edger(as.matrix(gene.counts[,9:12]),
    sample.list)
gene.data <- get.annotation("mm9","gene")
gene.filters <- get.defaults("gene.filter","mm9")
filter.results <- filter.genes(gene.counts,gene.data,
    gene.filters,sample.list)

pmoulos/metaseqr documentation built on Dec. 29, 2020, 5:56 a.m.