Function bh_fdr_threshold 

pub fn bh_fdr_threshold(pvalues: &[f64], q: f64) -> f64

Benjamini-Hochberg FDR step-up p-value cutoff for a family of p-values at false-discovery-rate level q. Returns the largest p(k) with p(k) <= (k/m)*q (sorted ascending); every p <= the returned cutoff is a discovery. Returns f64::NEG_INFINITY (reject nothing) for an empty family or when no rank satisfies the criterion.