augmentation: Wrapper function to the augmentation methods of the multtest...

Description Usage Arguments Details Value Author(s) References Examples

View source: R/Augmentation.R

Description

Wrapper function to the augmentation methods of the multtest package.

Usage

1
augmentation(adjPValues, newErrorControl, newK, newQ, silent=FALSE)

Arguments

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.

Details

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.

Value

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 errorControl, containing the type of error controlled by the function.

Author(s)

KornRohm, GillesBlanchard

References

Dudoit, S. and van der Laan, M.J. (2008) Multiple Testing Procedures with Applications to Genomics, Springer. (chapter 6)

Examples

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)

mutoss documentation built on May 2, 2019, 5:56 p.m.