webgestalt_lib 0.3.2

Library for computing enrichment for different analytes using ORA or GSEA.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45

struct Carrier {
    p: f64,
    original_order: usize,
}

#[derive(Clone)]
pub enum AdjustmentMethod {
    BH,
    None,
}

pub fn adjust(p_vals: &[f64], method: AdjustmentMethod) -> Vec<f64> {
    match method {
        AdjustmentMethod::BH => benjamini_hochberg(p_vals),
        AdjustmentMethod::None => p_vals.to_vec(),
    }
}

fn benjamini_hochberg(p_vals: &[f64]) -> Vec<f64> {
    let mut carriers: Vec<Carrier> = p_vals
        .iter()
        .enumerate()
        .map(|(i, p)| Carrier {
            p: *p,
            original_order: i,
        })
        .collect();
    carriers.sort_by(|a, b| a.p.partial_cmp(&b.p).unwrap());
    let m = carriers.len();
    let mut fdr_vals = vec![0.0; m];
    let mut prev_fdr = 1.0;
    for (i, carrier) in carriers.iter().enumerate().rev() {
        let mut fdr = carrier.p * m as f64 / (i + 1) as f64;
        if fdr > 1.0 {
            fdr = 1.0;
        }
        if fdr > prev_fdr {
            fdr = prev_fdr;
        } else {
            prev_fdr = fdr;
        }
        fdr_vals[carrier.original_order] = fdr;
    }
    fdr_vals
}