pub fn z_score(value: f64, mean: f64, stddev: f64) -> f64
Calculate the z-score: how many standard deviations a value is from the mean