Learn R Programming

sommer (version 1.1)

fdr: False Discovery Rate calculation

Description

This function calculates the FDR rate of 0.05 and is calculated using the qvalue package (Storey and Tibshirani 2003). The p-value corresponding to a q-value of 0.05 is determined by interpolation. The input is just a vector of p-values.

Usage

fdr(p, fdr.level = 0.05, type = "-log10")

Arguments

p
a vector of -log(p-values) or p-values
fdr.level
the level of false discovery rate desired
type
the type of p-values we are entering, either "-log10" or "raw" p-values

Value

  • [object Object]

References

Storey and Tibshirani. 2003. Statistical significance for genome-wide studies. PNAS 100:9440-9445

Examples

Run this code
pp <- rnorm(100,15,2) # -log10(p-values) obtained by GWAS
fdr(pp)

Run the code above in your browser using DataLab