R/fdr.r.R

Defines functions fdr.r

Documented in fdr.r

fdr.r = function(r, hatpai, fdr){
  res = sapply(hatpai, function(x){ifelse(x>r,1,0)})
  res = sum(res*(1-hatpai)) / sum(res)

  ifelse(!is.nan(res), return(res-fdr), return(1))
}

Try the BHTSpack package in your browser

Any scripts or data that you put into this service are public.

BHTSpack documentation built on Feb. 8, 2021, 5:07 p.m.