View source: R/posthocBySimes.R
posthocBySimes0Rcpp | R Documentation |
Lower bound on the number of correct rejections using Simes' reference family
posthocBySimes0Rcpp(p, select, alpha)
posthocBySimes(p, select, alpha, Rcpp = FALSE, verbose = FALSE)
p |
A numeric vector of |
select |
A vector of indices in |
alpha |
A numeric value, the significance level of the test procedure. |
Rcpp |
A boolean value: use Rcpp version of the implementation? Defaults to FALSE. |
verbose |
If |
If (R_k)_k provides jFWER control at level \alpha
then with
probability greater than 1-\alpha
, |H_0 cap R| \leq \min_k {|R
\cap (R_k)^c|+k-1}
A bit better: |H_0 cap R| \leq (\min_{k<= |R|} {|R
\cap R_k^c|+k-1}) \wedge |R|
A integer value, Simes's lower bound on the number of correct rejections within the selected hypotheses
posthocBySimes0Rcpp()
: Rcpp version
posthocBySimes()
: R version
Gilles Blanchard, Pierre Neuvial and Etienne Roquain
Blanchard, G., Neuvial, P., & Roquain, E. (2020). Post hoc confidence bounds on false positives using reference families. Annals of Statistics, 48(3), 1281-1303.
Goeman, J. J., & Solari, A. (2011). Multiple testing for exploratory research. Statistical Science, 26(4), 584-597.
m <- 1e3
m1 <- 200
p <- 1-pnorm(c(rnorm(m1, mean=4), rnorm(m-m1, mean=0)))
R <- union(1:10, sample(m, 10))
alpha <- 0.10
if (require("cherry")) {
hom <- hommelFast(p)
pickSimes(hom, R, silent=TRUE, alpha = alpha)
}
posthocBySimes(p, R, alpha=alpha)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.