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
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
use crate::evaluator::*;
#[derive(Clone, Debug)]
pub struct PercentDifferenceMagnitudePercentile {
quantile: f32,
name: String,
description: String,
}
lazy_info!(
PERCENT_DIFFERENCE_MAGNITUDE_PERCENTILE_INFO,
size: 1,
min_ts_length: 1,
t_required: false,
m_required: true,
w_required: false,
sorting_required: false,
);
impl PercentDifferenceMagnitudePercentile {
pub fn new(quantile: f32) -> Self {
assert!(
(quantile > 0.0) && (quantile < 0.5),
"quantiles should be between zero and half"
);
Self {
quantile,
name: format!(
"percent_difference_magnitude_percentile_{:.0}",
100.0 * quantile
),
description: format!(
"ratio of inter-percentile {:.3e}% - {:.3e}% range of magnitude to its mdeian",
100.0 * quantile,
100.0 * (1.0 - quantile),
),
}
}
pub fn set_name(&mut self, name: String) {
self.name = name;
}
#[inline]
pub fn default_quantile() -> f32 {
0.05
}
}
impl Default for PercentDifferenceMagnitudePercentile {
fn default() -> Self {
Self::new(Self::default_quantile())
}
}
impl<T> FeatureEvaluator<T> for PercentDifferenceMagnitudePercentile
where
T: Float,
{
fn eval(&self, ts: &mut TimeSeries<T>) -> Result<Vec<T>, EvaluatorError> {
self.check_ts_length(ts)?;
let nominator =
ts.m.get_sorted().ppf(1.0 - self.quantile) - ts.m.get_sorted().ppf(self.quantile);
let denominator = ts.m.get_median();
if nominator.is_zero() & denominator.is_zero() {
Err(EvaluatorError::ZeroDivision("median magnitude is zero"))
} else {
Ok(vec![nominator / denominator])
}
}
fn get_info(&self) -> &EvaluatorInfo {
&PERCENT_DIFFERENCE_MAGNITUDE_PERCENTILE_INFO
}
fn get_names(&self) -> Vec<&str> {
vec![self.name.as_str()]
}
fn get_descriptions(&self) -> Vec<&str> {
vec![self.description.as_str()]
}
}
#[cfg(test)]
#[allow(clippy::unreadable_literal)]
#[allow(clippy::excessive_precision)]
mod tests {
use super::*;
use crate::tests::*;
eval_info_test!(
percent_difference_magnitude_percentile_info,
PercentDifferenceMagnitudePercentile::default()
);
feature_test!(
percent_difference_magnitude_percentile,
[
Box::new(PercentDifferenceMagnitudePercentile::default()),
Box::new(PercentDifferenceMagnitudePercentile::new(0.05)),
Box::new(PercentDifferenceMagnitudePercentile::new(0.1)),
],
[4.85, 4.85, 4.6],
[
80.0_f32, 13.0, 20.0, 20.0, 75.0, 25.0, 100.0, 1.0, 2.0, 3.0, 7.0, 30.0, 5.0, 9.0,
10.0, 70.0, 80.0, 92.0, 97.0, 17.0
],
);
}