qvaluebh95: Computes Benjamini-Hochberg (95) q-value

In GenABEL: genome-wide SNP association analysis

Description

Computes Benjamini-Hochberg (95) q-value

Usage

qvaluebh95(p, fdrate=0.1)

Arguments

p	vector containing p-values
fdrate	desired FRD

Value

A list containing

pass	Is true if this P-value passed specified FDR
qvalue	qvalue

Author(s)

Yurii Aulchenko

Examples

require(GenABEL.data)
data(srdta)
a<-qtscore(qt2,data=srdta)
qv <- qvaluebh95(a[,"P1df"])
plot(a[,"Position"],-log10(qv$qvalue))

Example output

Loading required package: MASS
Loading required package: GenABEL.data

GenABEL documentation built on May 30, 2017, 3:36 a.m.