perGeneQValue: Summarize per-exon p-values into per-gene q-values.

In DEXSeq: Inference of differential exon usage in RNA-Seq

Description

The use case for this function is the following analysis: given per-exon p-values for null hypothesis H0, we can determine the number of genes in which at least for one exon H0 is rejected. What is the associated false disovery rate?

Usage

perGeneQValue(object, p = "pvalue", method = perGeneQValueExact)

Arguments

object	An DEXSeqResults object.
p	A character string indicating the name of the columns in DEXSeqResults from which to take the per-exon p-values.
method	Use the default value. This is for debugging only.

Details

Details

Value

A named numeric vector, values are per-gene q-values, names are gene.

Examples

  data(pasillaDEXSeqDataSet, package="pasilla")
  dxd <- estimateSizeFactors( dxd )
  dxd <- estimateDispersions( dxd )
  dxd <- testForDEU( dxd )
  dxr <- DEXSeqResults( dxd )

  perGeneQValue(dxr)

DEXSeq documentation built on Nov. 8, 2020, 5:11 p.m.