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

Description Usage Arguments Details Value Examples

View source: R/perGeneQValue.R

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

1
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

1
2
3
4
5
6
7
  data(pasillaDEXSeqDataSet, package="pasilla")
  dxd <- estimateSizeFactors( dxd )
  dxd <- estimateDispersions( dxd )
  dxd <- testForDEU( dxd )
  dxr <- DEXSeqResults( dxd )

  perGeneQValue(dxr)

areyesq89/DEXSeq documentation built on June 26, 2020, 2:31 a.m.