R: 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

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

[Package GenABEL version 1.8-0 Index]