Computes Benjamini-Hochberg (95) q-value

qvaluebh95 {GenABEL}R Documentation

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

Author(s)

Yurii Aulchenko

Examples

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

[Package GenABEL version 1.6-7 Index]