Function storey_qvalues 

pub fn storey_qvalues<T>(p_values: &[T], lambda: T) -> Result<Vec<T>>
where
    T: FloatOps,

Apply Storey’s q-value method for controlling false discovery rate

Storey’s q-value method estimates the proportion of true null hypotheses (π0) and uses this to obtain more powerful FDR control than the BH procedure.

Arguments

• p_values - A slice of p-values to adjust
• lambda - Tuning parameter for π0 estimation (between 0 and 1, typically 0.5)

Returns

• Result<Vec<f64>> - Vector of q-values

Example

let p_values = vec![0.01, 0.03, 0.05];
let qvalues = single_statistics::testing::correction::storey_qvalues(&p_values, 0.5).unwrap();