hypors 0.3.0

Hypothesis Testing with Rust
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

#[cfg(test)]
mod tests_anova {
    use hypors::anova::{anova, f_sample_size};

    const EPSILON: f64 = 0.001; // Tolerance for floating-point comparisons

    #[test]
    fn test_anova() {
        let data1 = vec![2.0, 3.0, 3.0, 5.0, 6.0];
        let data2 = vec![3.0, 4.0, 4.0, 6.0, 8.0];
        let data3 = vec![5.0, 6.0, 7.0, 8.0, 9.0];

        let result = anova(&[data1, data2, data3], 0.05).unwrap();

        let expected_f_statistic = 4.261;
        let expected_p_value = 0.039;
        let expected_null_hypothesis = "H0: µ1 = µ2 = µ3";
        let expected_alt_hypothesis = "Ha: At least one group mean is different";

        assert!((result.test_statistic - expected_f_statistic).abs() < EPSILON);
        assert!((result.p_value - expected_p_value).abs() < EPSILON);
        assert_eq!(result.reject_null, true);

        assert_eq!(result.null_hypothesis, expected_null_hypothesis);
        assert_eq!(result.alt_hypothesis, expected_alt_hypothesis);
    }

    #[test]
    fn test_anova_no_rejection() {
        let data1 = vec![2.0, 3.0, 4.0, 5.0, 6.0];
        let data2 = vec![3.0, 4.0, 5.0, 6.0, 7.0];
        let data3 = vec![4.0, 5.0, 6.0, 7.0, 8.0];

        let result = anova(&[data1, data2, data3], 0.05).unwrap();

        let expected_f_statistic = 2.0;
        let expected_p_value = 0.177;
        let expected_null_hypothesis = "H0: µ1 = µ2 = µ3";
        let expected_alt_hypothesis = "Ha: At least one group mean is different";

        assert!((result.test_statistic - expected_f_statistic).abs() < EPSILON);
        assert!((result.p_value - expected_p_value).abs() < EPSILON);
        assert_eq!(result.reject_null, false);

        assert_eq!(result.null_hypothesis, expected_null_hypothesis);
        assert_eq!(result.alt_hypothesis, expected_alt_hypothesis);
    }

    #[test]
    fn test_f_sample_size() {
        let effect_size = 0.25; // Cohen's f
        let alpha = 0.05; // 5% significance level
        let power = 0.80; // 80% power
        let num_groups = 3; // Number of groups

        let n = f_sample_size(effect_size, alpha, power, num_groups);
        let expected_sample_size = 21717.0;

        assert!((n - expected_sample_size).abs() <= 1.0);
    }
}