r2rs-stats 0.1.1

Statistics programming for Rust based on R's stats package
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

// "Whatever you do, work at it with all your heart, as working for the Lord,
// not for human masters, since you know that you will receive an inheritance
// from the Lord as a reward. It is the Lord Christ you are serving."
// (Col 3:23-24)

use std::str::FromStr;

use strafe_testing::{
    r::{RString, RTester},
    r_assert_relative_equal,
};
use strafe_type::FloatConstraint;

use crate::funcs::{
    ppoints::ppoints,
    var::{sd, variance},
};

pub fn sd_test_inner(x: &[f64]) {
    let rust_sd = sd(x);
    let r_sd = RTester::new()
        .set_script(&RString::from_string(format!(
            "x={};",
            RString::from_f64_slice(x),
        )))
        .set_display(&RString::from_str("sd(x)").unwrap())
        .run()
        .expect("R running error")
        .as_f64()
        .expect("R running error");
    r_assert_relative_equal!(rust_sd, r_sd);
}

pub fn variance_test_inner(x: &[f64]) {
    let rust_variance = variance(x);
    let r_variance = RTester::new()
        .set_script(&RString::from_string(format!(
            "x={};",
            RString::from_f64_slice(x),
        )))
        .set_display(&RString::from_str("var(x)").unwrap())
        .run()
        .expect("R running error")
        .as_f64()
        .expect("R running error");
    r_assert_relative_equal!(rust_variance, r_variance);
}

pub fn ppoints_test_inner(n: usize) {
    let rust_ppoints = ppoints(n, None);
    let r_ppoints = RTester::new()
        .set_script(&RString::from_string(format!(
            "n={};",
            RString::from_f64(n as f64),
        )))
        .set_display(&RString::from_str("ppoints(n)").unwrap())
        .run()
        .expect("R running error")
        .as_f64_vec()
        .expect("R running error");
    for (rust_ret, r_ret) in rust_ppoints.into_iter().zip(r_ppoints.into_iter()) {
        r_assert_relative_equal!(rust_ret.unwrap(), r_ret);
    }
}