Description Usage Arguments Value Author(s) References Examples
Calculates adjusted p-values of Hommel's method efficiently.
1 | hommel (p, simes = TRUE)
|
p |
A vector of p-values. |
simes |
If |
Returns a hommel-class
object.
Jelle Goeman, Rosa Meijer, Thijmen Krebs.
Hommel, G. (1988). A stagewise rejective multiple test procedure based on a modified Bonferroni test. Biometrika, 75(2), 383-386.
Meijer, R. J., Krebs, T. J., & Goeman, J. J. (2019). Hommel's procedure in linear time. Biometrical Journal, 61(1), 73-82.
Goeman, J., Meijer, R., Krebs, T., & Solari, A. (2016). Simultaneous control of all false discovery proportions in large-scale multiple hypothesis testing. arXiv preprint arXiv:1611.06739.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 | #Generate a vector of pvalues
set.seed(102)
m <- 10
pvalues <- c(runif(0.5*m,0,0.02), runif(0.5*m,0,1))
# First step: create a hommel object.
# Use Simes'test (simes = TRUE) or Hommel's robust test (simes = FALSE)
hom <- hommel(pvalues, simes = TRUE)
hom
summary(hom)
# retrieve familywise error adjusted p-values
p.adjust(hom)
# Find lower confidence bounds
# for the number of false null hypotheses (discoveries)
# In any subset of interest
discoveries(hom, 1:2)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.