qvalue.cal(p, p0, version = 1)
version=2
, the original version of
the q-value, i.e. min{pFDR}, will be computed. if
version=1
, min{FDR} will be used in the computation
of the q-value.p
containing the q-values
corresponding to the p-values in p
.
version = 1
in qvalue.cal
corresponds to setting
robust = FALSE
in the function qvalue
of John Storey's
R package qvalue, while version = 2
corresponds to
robust = TRUE
.
Storey, J.D., and Tibshirani, R. (2003). Statistical Significance for Genome-wide Studies. PNAS, 100, 9440-9445.
pi0.est
,SAM-class
,sam