From 52e40995b7adfec8c1a837bed5a84d475f394602 Mon Sep 17 00:00:00 2001
From: JamesP <jlp765@gmail.com>
Date: Thu, 12 Nov 2015 06:33:26 +1000
Subject: [PATCH] added helper procs for standalone calc of stats for array/seq
 Changed derivation to use runningstats, so there is a single point of truth
 for the calculations, and provides more stats functions

---
 lib/pure/stats.nim | 45 +++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 45 insertions(+)

diff --git a/lib/pure/stats.nim b/lib/pure/stats.nim
index 5811e141e6..43d9c35f63 100644
--- a/lib/pure/stats.nim
+++ b/lib/pure/stats.nim
@@ -166,6 +166,51 @@ proc `+`*(a, b: RunningStat): RunningStat =
 proc `+=`*(a: var RunningStat, b: RunningStat) {.inline.} =
   ## add a second RunningStats `b` to `a`
   a = a + b
+# ---------------------- standalone array/seq stats ---------------------
+proc mean*[T](x: openArray[T]): float =
+  var rs: RunningStat
+  rs.push(x)
+  result = rs.mean()
+
+proc variance*[T](x: openArray[T]): float =
+  var rs: RunningStat
+  rs.push(x)
+  result = rs.variance()
+
+proc varianceS*[T](x: openArray[T]): float =
+  var rs: RunningStat
+  rs.push(x)
+  result = rs.varianceS()
+
+proc standardDeviation*[T](x: openArray[T]): float =
+  var rs: RunningStat
+  rs.push(x)
+  result = rs.standardDeviation()
+
+proc standardDeviationS*[T](x: openArray[T]): float =
+  var rs: RunningStat
+  rs.push(x)
+  result = rs.standardDeviationS()
+
+proc skewness*[T](x: openArray[T]): float =
+  var rs: RunningStat
+  rs.push(x)
+  result = rs.skewness()
+
+proc skewnessS*[T](x: openArray[T]): float =
+  var rs: RunningStat
+  rs.push(x)
+  result = rs.skewnessS()
+
+proc kurtosis*[T](x: openArray[T]): float =
+  var rs: RunningStat
+  rs.push(x)
+  result = rs.kurtosis()
+
+proc kurtosisS*[T](x: openArray[T]): float =
+  var rs: RunningStat
+  rs.push(x)
+  result = rs.kurtosisS()
 
 # ---------------------- Running Regression -----------------------------