mutoss (version 0.1-13)

rom: Rom's (1990) step-up-procedure.

Description

Rom's step-up-procedure is applied to pValues. The procedure controls the FWER in the strong sense if the pValues are stochastically independent.

Usage

rom(pValues, alpha, silent=FALSE)

Value

A list containing:

adjPValues

A numeric vector containing the adjusted pValues

rejected

A logical vector indicating which hypotheses are rejected.

criticalValues

A numeric vector containing critical values used in the step-up-down test.

errorControl

A Mutoss S4 class of type errorControl, containing the type of error controlled by the function and the level alpha.

Author

Marsel Scheer

Arguments

pValues

pValues to be used. They are assumed to be stochastically independent.

alpha

the level at which the FWER shall be controlled.

silent

if true any output on the console will be suppressed.

Details

This function calculates the critical values by the formula given in Finner, H. and Roters, M. (2002) based on the joint distribution of order statistics. After that a step-up test is performed to reject hypotheses associated with pValues.

Since the formula for the critical values is recursive, the calculation of adjusted pValues is far from obvious and is not implemented here.

References

Rom, D. M. (1990). A sequentially rejective test procedure based on a modified Bonferroni inequality. Biometrika 77, 663-665.

Finner, H. and Roters, M. (2002). Multiple hypotheses testing and expected type I errors. Ann. Statist. 30, 220-238.

Examples

r <- c(runif(50), runif(50, 0, 0.01))
result <- rom(r, 0.05)
result <- rom(r, 0.05, silent = TRUE)