geopagg 0.3.1

An implementation of the geopagg algorithm for p-value aggregation of sgRNAs across gene groupings
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
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72

use bon::Builder;

#[derive(Debug, Clone, Builder)]
pub struct GeneResult {
    pub gene: String,
    pub wgm: f64,
    pub logfc: f64,
    #[builder(default = false)]
    pub amalgam: bool,
    #[builder(default = 1.0)]
    pub empirical_fdr: f64,
    #[builder(default = 1.0)]
    pub adjusted_empirical_fdr: f64,
    #[builder(skip = -logfc * wgm.ln())]
    pub product: f64,
}

pub struct GeoPAGGResults {
    pub genes: Vec<String>,
    pub wgms: Vec<f64>,
    pub logfcs: Vec<f64>,
    pub amalgams: Vec<bool>,
    pub empirical_fdr: Vec<f64>,
    pub adjusted_empirical_fdr: Vec<f64>,
}
impl GeoPAGGResults {
    #[must_use]
    pub fn from_vec(mut gene_results: Vec<GeneResult>) -> Self {
        let mut genes = vec![];
        let mut wgms = vec![];
        let mut logfcs = vec![];
        let mut amalgams = vec![];
        let mut empirical_fdr = vec![];
        let mut adjusted_empirical_fdr = vec![];

        // Sort the results by gene name
        gene_results.sort_unstable_by(|a, b| a.gene.cmp(&b.gene));

        for gene_result in gene_results {
            genes.push(gene_result.gene);
            wgms.push(gene_result.wgm);
            logfcs.push(gene_result.logfc);
            amalgams.push(gene_result.amalgam);
            empirical_fdr.push(gene_result.empirical_fdr);
            adjusted_empirical_fdr.push(gene_result.adjusted_empirical_fdr);
        }

        Self {
            genes,
            wgms,
            logfcs,
            amalgams,
            empirical_fdr,
            adjusted_empirical_fdr,
        }
    }

    pub fn pprint(&self) {
        println!("Gene\tWGM\tLogFC\tAmalgam\tEmpirical FDR\tAdjusted Empirical FDR");
        for i in 0..self.genes.len() {
            println!(
                "{}\t{}\t{}\t{}\t{}\t{}",
                self.genes[i],
                self.wgms[i],
                self.logfcs[i],
                self.amalgams[i],
                self.empirical_fdr[i],
                self.adjusted_empirical_fdr[i]
            );
        }
    }
}