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

Description Usage Arguments Value Author(s) Examples

Description

Computes Benjamini-Hochberg (95) q-value

Usage

1
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

1
2
3
4
5
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.

Related to qvaluebh95 in GenABEL...