hommelFast {cherry} | R Documentation |
Calculates adjusted p-values of Hommel's method efficiently.
hommelFast (pvalues, simes = TRUE)
pvalues |
A vector of p-values. |
simes |
If |
Returns a hommel
object.
Rosa Meijer: r.j.meijer@lumc.nl
#Generate a vector of pvalues
set.seed(1)
n <- 1000
pvalues <- c(runif(0.50*n,0,0.001), runif(0.50*n,0,1))
#Create an hommel object in which the adjusted p-values are stored, based on a Simes'test
#(or Hommel's test, by choosing simes = FALSE):
hom <- hommelFast(pvalues, simes = TRUE)
#Retrieve the first 10 adjusted p-values by using the \code{p.adjust} method
# from the hommel package. Note that they are not sorted
hommel::p.adjust(hom)[1:10]