ci-adjust: CI Adjust

In Rariant: Identification and Assessment of Single Nucleotide Variants through Shifts in Non-Consensus Base Call Frequencies

Description

Multiple testing adjustment of confidence levels, as proposed by Benjamini and Yekutieli.

Usage

  ciAdjustLevel(eta0, conf_level)

Arguments

eta0	Estimated fraction of tests that are consistent with the null hypothesis.
conf_level	Unadjusted confidence level

Value

The adjusted confidence level.

References

Benjamini, Yoav, and Daniel Yekutieli. False Discovery Rate–adjusted Multiple Confidence Intervals for Selected Parameters. Journal of the American Statistical Association 100, no. 469 (2005): 71–81.

Examples

conf_level = 0.95
eta0 = seq(0, 1, by = 0.02)

conf_level_adj = ciAdjustLevel(eta0, conf_level)

plot(eta0, conf_level_adj, pch = 20, ylim = c(conf_level, 1))

Rariant documentation built on Nov. 8, 2020, 6:56 p.m.