R/get_bh_threshold.R

In LPRelevance: Relevance-Integrated Statistical Inference Engine

Defines functions get_bh_threshold

Documented in get_bh_threshold

get_bh_threshold <-
function(pvals, alpha, mtests = length(pvals)){
  m <- length(pvals)
  pvals <- sort(pvals)
  prejected <- which(pvals <= (1:m)/mtests*alpha)
  ifelse(length(prejected)==0, 0, pvals[prejected[which.max(prejected)]])
}

Try the LPRelevance package in your browser