GenABEL (version 1.8-0)

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

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

Examples

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