powered by
This function calculates Q-value based on Storey.
calculateQvalue(df, lambda = seq(0.05, 0.5, 0.01))
data.frame
a data frame containng p-values
a vector of the lambda values utilized to obtain pi0.lambda
Atsushi Fukushima
Storey JD, The Annals of Statistics 31:2013-2035 (2003)
data(example.df) p.value <- example.df$p.value df <- data.frame(p.value = p.value) res <- calculateQvalue(df)
Run the code above in your browser using DataLab