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
|
package de.lmu.ifi.dbs.elki.math.statistics.distribution;
import static org.junit.Assert.assertEquals;
import static org.junit.Assert.assertTrue;
/**
* Abstract base class for distribution unit testing.
*
* @author Erich Schubert
*/
public class AbstractDistributionTest {
public void checkPDF(Distribution d, double[] x, double[] expected, double err) {
int maxerrlev = Integer.MIN_VALUE;
for(int i = 0; i < x.length; i++) {
double val = d.pdf(x[i]);
if(val == expected[i]) {
continue;
}
double diff = Math.abs(val - expected[i]);
final int errlev = (int) Math.ceil(Math.log10(diff / expected[i]));
maxerrlev = Math.max(errlev, maxerrlev);
if(diff < err || diff / expected[i] < err) {
continue;
}
assertEquals("Error magnitude: 1e" + errlev + " at " + x[i], expected[i], val, err);
}
int given = (int) Math.floor(Math.log10(err * 1.1));
// if (given > maxerrlev) {
// System.err.println("PDF Error for "+d+" magnitude is not tight: expected "+maxerrlev+" got "+given);
// }
assertTrue("Error magnitude is not tight: expected " + maxerrlev + " got " + given, given <= maxerrlev);
}
public void checkCDF(Distribution d, double[] x, double[] expected, double err) {
int maxerrlev = Integer.MIN_VALUE;
for(int i = 0; i < x.length; i++) {
double val = d.cdf(x[i]);
if(val == expected[i]) {
continue;
}
double diff = Math.abs(val - expected[i]);
final int errlev = (int) Math.ceil(Math.log10(diff / expected[i]));
maxerrlev = Math.max(errlev, maxerrlev);
if(diff < err || diff / expected[i] < err) {
continue;
}
assertEquals("Error magnitude: 1e" + errlev + " at " + x[i], expected[i], val, err);
}
int given = (int) Math.floor(Math.log10(err * 1.1));
// if (given > maxerrlev) {
// System.err.println("CDF Error for "+d+" magnitude is not tight: expected "+maxerrlev+" got "+given);
// }
assertTrue("Error magnitude is not tight: expected " + maxerrlev + " got " + given, given <= maxerrlev);
}
public void checkQuantile(Distribution d, double[] x, double[] expected, double err) {
int maxerrlev = Integer.MIN_VALUE;
for(int i = 0; i < x.length; i++) {
double val = d.quantile(x[i]);
if(val == expected[i]) {
continue;
}
double diff = Math.abs(val - expected[i]);
final int errlev = (int) Math.ceil(Math.log10(diff / expected[i]));
maxerrlev = Math.max(errlev, maxerrlev);
if(diff < err || diff / expected[i] < err) {
continue;
}
assertEquals("Error magnitude: 1e" + errlev + " at " + x[i], expected[i], val, err);
}
int given = (int) Math.floor(Math.log10(err * 1.1));
// if (given > maxerrlev) {
// System.err.println("Probit Error for "+d+" magnitude is not tight: expected "+maxerrlev+" got "+given);
// }
assertTrue("Error magnitude is not tight: expected " + maxerrlev + " got " + given, given <= maxerrlev);
}
}
|