stat.deseq: Statistical testing with DESeq

Description Usage Arguments Value Author(s) Examples

View source: R/metaseqr.stat.R

Description

This function is a wrapper over DESeq statistical testing. It accepts a matrix of normalized gene counts or an S4 object specific to each normalization algorithm supported by metaseqR.

Usage

1
2
    stat.deseq(object, sample.list, contrast.list = NULL,
        stat.args = NULL)

Arguments

object

a matrix or an object specific to each normalization algorithm supported by metaseqR, containing normalized counts. Apart from matrix (also for NOISeq), the object can be a SeqExpressionSet (EDASeq), CountDataSet (DESeq) or DGEList (edgeR).

sample.list

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

contrast.list

a named structured list of contrasts as returned by make.contrast.list or just the vector of contrasts as defined in the main help page of metaseqr.

stat.args

a list of DESeq statistical algorithm parameters. See the result of get.defaults("statistics", "deseq") for an example and how you can modify it. It is not required when the input object is already a CountDataSet from DESeq normalization as the dispersions are already estimated.

Value

A named list of p-values, whose names are the names of the contrasts.

Author(s)

Panagiotis Moulos

Examples

1
2
3
4
5
6
require(DESeq)
data.matrix <- counts(makeExampleCountDataSet())
sample.list <- list(A=c("A1","A2"),B=c("B1","B2","B3"))
contrast <- "A_vs_B"
norm.data.matrix <- normalize.deseq(data.matrix,sample.list)
p <- stat.deseq(norm.data.matrix,sample.list,contrast)

metaseqR documentation built on Nov. 8, 2020, 5:57 p.m.