bounding.seq {HCTR}

Bounding Sequence

Description

Calculates bounding sequence of higher crticism for proportion estimator using p-values

Usage

bounding.seq(p.value, alpha)

Arguments

p.value

A matrix of p-values from permutation: row is from each permutation; column is from each variable.

alpha

Probability of Type I error for bounding sequence, the default value is 1 / sqrt(log(p)), where p is number of p-values in each permutation.

Value

A bounding value of higher criticism with (1 - alpha) confidence.

References

Jeng XJ, Zhang T, Tzeng J (2019). “Efficient Signal Inclusion With Genomic Applications.” Journal of the American Statistical Association, 1–23.

Examples

set.seed(10)
X <- matrix(runif(n = 10000, min = 0, max = 1), nrow = 100)
result <- bounding.seq(p.value = X)