This function calculates Q-value based on Storey.

calculateQvalue(df, lambda = seq(0.05, 0.5, 0.01))

Arguments

df

a data frame containng p-values

lambda

a vector of the lambda values utilized to obtain pi0.lambda

Value

data.frame

References

Storey JD, The Annals of Statistics 31:2013-2035 (2003)

Author

Atsushi Fukushima

Examples

data(example.df)
p.value <- example.df$p.value
df <- data.frame(p.value = p.value)
res <- calculateQvalue(df)