basic_stats 1.0.0

Lightweight library with basic statistical estimators and hypothesis tests.
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

//! Example of use of [`basic_stats::aok`] traits.
//! Requires features **`aok`** and **`normal`**.

use basic_stats::{
    aok::{AokBasicStats, AokFloat},
    core::{AltHyp, SampleMoments},
    normal::{welch_alt_hyp_ci, welch_p, welch_test},
};

fn main() {
    let x = [14., 15., 15., 15., 16., 18., 22., 23., 24., 25., 25.];
    let y = [
        10., 12., 14., 15., 18., 22., 24., 27., 31., 33., 34., 34., 34.,
    ];

    let moments_x = SampleMoments::from_slice(&x);
    let moments_y = SampleMoments::from_slice(&y);
    let alt_hyp = AltHyp::Gt;

    {
        println!("*** Ok scenario:");

        let alpha = 0.05;

        // Welch function calls below return Ok prior to invocation of aok().

        let p = welch_p(&moments_x, &moments_y, alt_hyp).aok();
        println!("p={p}");
        let ci = welch_alt_hyp_ci(&moments_x, &moments_y, alt_hyp, alpha).aok();
        println!("ci={ci:?}");
        let test_res = welch_test(&moments_x, &moments_y, alt_hyp, alpha).aok();
        println!("test_res={test_res:?}");
    }

    {
        println!("*** Err scenario:");

        let alpha = 1.0;

        // Welch function calls below return Err prior to invocation of aok().

        let p = welch_p(&moments_x, &SampleMoments::default(), alt_hyp).aok();
        println!("p={p}");
        let ci = welch_alt_hyp_ci(&moments_x, &moments_y, alt_hyp, alpha).aok();
        println!("ci={ci:?}");
        let test_res = welch_test(&moments_x, &moments_y, alt_hyp, alpha).aok();
        println!("test_res={test_res:?}");
    }
}