DEXSeq (version 1.18.4)

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

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.

Value

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

Details

Details

Examples

Run this code

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

  perGeneQValue(dxr)


Run the code above in your browser using DataLab