R/incorrect_bh.R

In ndbrown6/MSK-GRAIL-TECHVAL: High-intensity sequencing reveals the sources of plasma circulating cell-free DNA variants

'incorrect_bh' <- function(p_vals)
{
	base_rate = 0.05
	ox = order(p_vals, decreasing=FALSE)
	xthreshold = base_rate*(1:length(ox))/length(ox)
	over_flag = (xthreshold>p_vals[ox])
	return(invisible(over_flag[ox]))
}