From aaa36569ded69399f73848440e18329f7457354b Mon Sep 17 00:00:00 2001
From: Andy Eschbacher <andy.eschbacher@gmail.com>
Date: Mon, 26 Sep 2016 16:26:34 -0400
Subject: [PATCH] first add

---
 src/pg/sql/18_outliers.sql | 52 ++++++++++++++++++++++++++++++++++++++
 1 file changed, 52 insertions(+)
 create mode 100644 src/pg/sql/18_outliers.sql

diff --git a/src/pg/sql/18_outliers.sql b/src/pg/sql/18_outliers.sql
new file mode 100644
index 0000000..48588ce
--- /dev/null
+++ b/src/pg/sql/18_outliers.sql
@@ -0,0 +1,52 @@
+CREATE OR REPLACE FUNCTION CDB_StaticOutlier(attr numeric, threshold numeric)
+RETURNS numeric
+AS $$
+BEGIN
+
+  RETURN attr > threshold;
+
+END;
+$$ LANGUAGE plpgsql;
+
+
+CREATE OR REPLACE FUNCTION CDB_PercentOutlier(attr numeric[], outlier_fraction numeric, ids int[])
+RETURNS TABLE(outlier boolean, rowid int)
+AS $$
+DECLARE
+  avg_val numeric;
+  out_vals boolean[];
+BEGIN
+
+  SELECT avg(i) INTO avg_val FROM unnest(attr) As x(i);
+
+  SELECT array_agg( i > avg_val * outlier_fraction) INTO out_vals
+  FROM unnest(attr) As x(i);
+
+  RETURN QUERY
+  SELECT unnest(out_vals) As outlier,
+         unnest(ids) As rowid;
+
+END;
+$$ LANGUAGE plpgsql;
+
+CREATE OR REPLACE FUNCTION CDB_StdDevOutlier(attrs numeric[], num_deviations numeric, ids int[])
+RETURNS TABLE(outlier boolean, rowid int)
+AS $$
+DECLARE
+  stddev_val numeric;
+  avg_val numeric;
+  out_vals boolean[];
+BEGIN
+
+  SELECT stddev(i), avg(i) INTO stddev_val, avg_val
+  FROM unnest(attrs) As x(i);
+
+  SELECT array_agg(abs(i - avg_val) / stddev_val > num_deviations) INTO out_vals
+  FROM unnest(attrs) As x(i);
+
+
+  RETURN QUERY
+  SELECT unnest(out_vals) As outlier,
+         unnest(ids) As rowid;
+END;
+$$ LANGUAGE plpgsql;