Description Usage Arguments Details Value Author(s) References Examples
Wrapper function to the augmentation methods of the multtest package.
1 | augmentation(adjPValues, newErrorControl, newK, newQ, silent=FALSE)
|
adjPValues |
a vector of p-values that are *already* adjusted for FWER |
newErrorControl |
new error control type to adjust for. One of ("gFWER", "FDX", "FDR") |
newK |
k-parameter if newErrorControl is "gFWER" |
newQ |
q-parameter if newErrorControl is "FDX" |
silent |
if true any output on the console will be suppressed. |
The augmentation method turns a vector of p-values which are already adjusted for FWER control into p-values that are adjusted for gFWER, FDX or FDR. The underlying idea (for gFWER and FDX) is that the set of hypotheses rejected at a given level alpha under FWER can be "augmented" by rejecting some additional hypotheses while still ensuring (strong) control of the desired weaker type I criterion. For FDR, it uses the fact that FDX control for q=alpha=1-sqrt(1-beta) entails FDR control at level beta.
Use of these augmentation methods is recommended only in the situation where FWER-controlled p-values are directly available from the data (using some specific method). When only marginal p-values are available, it is generally prerefable to use other adjustment methods directly aimed at the intended criterion (as opposed to first adjust for FWER, then augment)
Note: in the multtest package, two methods ("restricted" and "conservative") are available for FDR augmentation. Here the 'restricted' method is forced for FDR augmentation since it is in fact always valid and better than "conservative" (M. van der Laan, personal communication) with respect to power.
A list containing:
adjPValues |
A numeric vector containing the new adjusted pValues |
rejected |
A logical vector indicating which hypotheses are rejected. Currently always NULL. |
errorControl |
A Mutoss S4 class of type |
KornRohm, GillesBlanchard
Dudoit, S. and van der Laan, M.J. (2008) Multiple Testing Procedures with Applications to Genomics, Springer. (chapter 6)
1 2 3 4 5 6 | ## Not run: TODO NH (MS) EXAMPLE PROBLEM
p <- c(runif(50), runif(50, 0, 0.01))
fwer_adj <- bonferroni(p) # adjust for FWER using Bonferroni correction
fdx_adj <- augmentation(fwer_adj$adjPValues , "FDX", q=0.1, silent = TRUE) #augment to FDX (q=0.1)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.