census-proteomics 0.3.3

Rust library for working with proteomics data quantified by the Census algorithm
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

use std::u32;

/// Calaculate the sum of a slice
#[inline]
pub fn sum(slice: &[f64]) -> f64 {
    slice.iter().fold(0.0f64, |acc, x| acc + x)
}

/// Calculate the mean value of a slice
#[inline]
pub fn mean(slice: &[u32]) -> f64 {
    slice.iter().sum::<u32>() as f64 / slice.len() as f64
}

/// Calculate the mean value of a slice
#[inline]
pub fn meanf(slice: &[f64]) -> f64 {
    slice.iter().sum::<f64>() as f64 / slice.len() as f64
}

/// Return the maximum value of a slice
#[inline]
pub fn max(slice: &[u32]) -> u32 {
    slice
        .iter()
        .fold(u32::MIN, |acc, &x| if x > acc { x } else { acc })
}

/// Calculate the standard deviation (population) of a slice
#[inline]
pub fn stddev(slice: &[u32]) -> f64 {
    let mean = mean(slice);
    (slice
        .iter()
        .fold(0.0f64, |acc, x| acc + (*x as f64 - mean).powi(2))
        / slice.len() as f64)
        .sqrt()
}

/// Calculate the standard deviation (population) of a slice
#[inline]
pub fn stddevf(slice: &[f64]) -> f64 {
    let mean = meanf(slice);
    (slice
        .iter()
        .fold(0.0f64, |acc, x| acc + (*x - mean).powi(2))
        / slice.len() as f64)
        .sqrt()
}

/// Calculate the standard error (population) of a slice
#[inline]
pub fn stderr(slice: &[u32]) -> f64 {
    stddev(slice) / (slice.len() as f64).sqrt()
}

pub fn cv(slice: &[u32]) -> f64 {
    stddev(slice) / mean(slice)
}