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