| qvaluebh95 {GenABEL} | R Documentation |
Computes Benjamini-Hochberg (95) q-value
qvaluebh95(p, fdrate=0.1)
p |
vector containing p-values |
fdrate |
desired FRD |
A list containing
pass |
Is true if this P-value passed specified FDR |
qvalue |
qvalue |
Yurii Aulchenko
require(GenABEL.data) data(srdta) a<-qtscore(qt2,data=srdta) qv <- qvaluebh95(a[,"P1df"]) plot(a[,"Position"],-log10(qv$qvalue))