cdflib 0.4.1

Pure-Rust port of CDFLIB: cumulative distribution functions (CDF) associated with common probability distributions
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
73
74
75
76
77

#![cfg(not(miri))]

//! Reference-table tests for the Beta family of distributions: Beta,
//! Student's *t*, Fisher–Snedecor (F). Each table was produced by the
//! Fortran `cumbet`/`cumt`/`cumf` routines from
//! `tests/regenerate/refs/cdflib.f90`.

mod common;

use cdflib::{Beta, ContinuousCdf, FisherSnedecor, StudentsT};
use common::{assert_close_eps, read_csv, DISTRIBUTION_ABS_TOL, DISTRIBUTION_REL_TOL};

#[test]
fn beta_cdf_matches_cumbet_reference() {
    for row in read_csv("tests/data/beta_cdf.csv") {
        let [a, b, x, expected_cdf, expected_sf] = row[..] else {
            panic!("width");
        };
        let d = Beta::new(a, b);
        assert_close_eps(
            d.cdf(x),
            expected_cdf,
            DISTRIBUTION_REL_TOL,
            DISTRIBUTION_ABS_TOL,
        );
        assert_close_eps(
            d.ccdf(x),
            expected_sf,
            DISTRIBUTION_REL_TOL,
            DISTRIBUTION_ABS_TOL,
        );
    }
}

#[test]
fn students_t_cdf_matches_cumt_reference() {
    for row in read_csv("tests/data/students_t_cdf.csv") {
        let [df, t, expected_cdf, expected_sf] = row[..] else {
            panic!("width");
        };
        let d = StudentsT::new(df);
        assert_close_eps(
            d.cdf(t),
            expected_cdf,
            DISTRIBUTION_REL_TOL,
            DISTRIBUTION_ABS_TOL,
        );
        assert_close_eps(
            d.ccdf(t),
            expected_sf,
            DISTRIBUTION_REL_TOL,
            DISTRIBUTION_ABS_TOL,
        );
    }
}

#[test]
fn f_cdf_matches_cumf_reference() {
    for row in read_csv("tests/data/f_cdf.csv") {
        let [dfn, dfd, fx, expected_cdf, expected_sf] = row[..] else {
            panic!("width");
        };
        let d = FisherSnedecor::new(dfn, dfd);
        assert_close_eps(
            d.cdf(fx),
            expected_cdf,
            DISTRIBUTION_REL_TOL,
            DISTRIBUTION_ABS_TOL,
        );
        assert_close_eps(
            d.ccdf(fx),
            expected_sf,
            DISTRIBUTION_REL_TOL,
            DISTRIBUTION_ABS_TOL,
        );
    }
}