Mercurial > dive4elements > river
comparison artifacts/src/main/java/org/dive4elements/river/artifacts/model/sq/Outlier.java @ 5838:5aa05a7a34b7
Rename modules to more fitting names.
| author | Sascha L. Teichmann <teichmann@intevation.de> |
|---|---|
| date | Thu, 25 Apr 2013 15:23:37 +0200 |
| parents | flys-artifacts/src/main/java/org/dive4elements/river/artifacts/model/sq/Outlier.java@bd047b71ab37 |
| children | 4897a58c8746 |
comparison
equal
deleted
inserted
replaced
| 5837:d9901a08d0a6 | 5838:5aa05a7a34b7 |
|---|---|
| 1 package org.dive4elements.river.artifacts.model.sq; | |
| 2 | |
| 3 import org.dive4elements.river.artifacts.math.GrubbsOutlier; | |
| 4 import org.dive4elements.river.artifacts.math.StdDevOutlier; | |
| 5 | |
| 6 import java.util.ArrayList; | |
| 7 import java.util.List; | |
| 8 | |
| 9 import org.apache.commons.math.MathException; | |
| 10 | |
| 11 import org.apache.log4j.Logger; | |
| 12 | |
| 13 public class Outlier | |
| 14 { | |
| 15 private static Logger log = Logger.getLogger(Outlier.class); | |
| 16 | |
| 17 private static final String GRUBBS = "outlier.method.grubbs"; | |
| 18 | |
| 19 //private static final String STD_DEV = "std-dev"; | |
| 20 | |
| 21 public interface Callback { | |
| 22 | |
| 23 void initialize(List<SQ> sqs) throws MathException; | |
| 24 | |
| 25 double eval(SQ sq); | |
| 26 | |
| 27 void iterationFinished( | |
| 28 double stdDev, | |
| 29 SQ outlier, | |
| 30 List<SQ> remaining); | |
| 31 | |
| 32 } // interface Callback | |
| 33 | |
| 34 public static void detectOutliers( | |
| 35 Callback callback, | |
| 36 List<SQ> sqs, | |
| 37 double stdDevFactor, | |
| 38 String method | |
| 39 ) | |
| 40 throws MathException | |
| 41 { | |
| 42 boolean debug = log.isDebugEnabled(); | |
| 43 | |
| 44 if (method == null) { | |
| 45 method = "std-dev"; | |
| 46 } | |
| 47 | |
| 48 if (debug) { | |
| 49 log.debug("stdDevFactor: " + stdDevFactor); | |
| 50 log.debug("method: " + method); | |
| 51 } | |
| 52 | |
| 53 List<SQ> data = new ArrayList<SQ>(sqs); | |
| 54 | |
| 55 double [] stdDev = new double[1]; | |
| 56 | |
| 57 boolean useGrubbs = method.equals(GRUBBS); | |
| 58 | |
| 59 if (useGrubbs) { | |
| 60 stdDevFactor = Math.max(0d, Math.min(stdDevFactor/100d, 1d)); | |
| 61 } | |
| 62 | |
| 63 List<Double> values = new ArrayList<Double>(data.size()); | |
| 64 | |
| 65 while (data.size() > 2) { | |
| 66 | |
| 67 callback.initialize(data); | |
| 68 | |
| 69 for (SQ sq: data) { | |
| 70 values.add(callback.eval(sq)); | |
| 71 } | |
| 72 | |
| 73 Integer ndx = useGrubbs | |
| 74 ? GrubbsOutlier.findOutlier(values, stdDevFactor, stdDev) | |
| 75 : StdDevOutlier.findOutlier(values, stdDevFactor, stdDev); | |
| 76 | |
| 77 if (ndx == null) { | |
| 78 callback.iterationFinished(stdDev[0], null, data); | |
| 79 break; | |
| 80 } | |
| 81 | |
| 82 SQ outlier = data.remove(ndx.intValue()); | |
| 83 if (debug) { | |
| 84 log.debug("stdDev: " + stdDev[0]); | |
| 85 log.debug("removed " + ndx + | |
| 86 "; S: " + outlier.getS() + " Q: " + outlier.getQ()); | |
| 87 } | |
| 88 callback.iterationFinished(stdDev[0], outlier, data); | |
| 89 values.clear(); | |
| 90 } | |
| 91 } | |
| 92 } | |
| 93 // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf8 : |