From ddcb1f73423a9b47d9b457cf3b19f622b0aa567f Mon Sep 17 00:00:00 2001
From: wukong <wukong@longaeva>
Date: Sat, 5 Jan 2019 16:33:14 -0800
Subject: added min and max to mean and mean_avg, fixed regex bugs in mean

---
 mean_avg.awk | 10 ++++++----
 1 file changed, 6 insertions(+), 4 deletions(-)

(limited to 'mean_avg.awk')

diff --git a/mean_avg.awk b/mean_avg.awk
index ed6f81a..4bf7051 100644
--- a/mean_avg.awk
+++ b/mean_avg.awk
@@ -24,22 +24,24 @@ NF > 0 {
     for (n=1; n <= NF; n++) {
         if ($n ~ number) {
             count[n] += 1
+            (count[n] == 1 || $n < min[n]) ? min[n] = $n : min[n] = min[n]
+            (count[n] == 1 || $n > max[n]) ? max[n] = $n : max[n] = max[n]
             delta0[n] = $n - mean[n]
             mean[n] += delta0[n]/count[n]
             delta1[n] = $n - mean[n]
             M2[n] += delta0[n]*delta1[n]
-            (count[n] > 1) ? var[n] = M2[n]/(count[n] - 1) : var[n] = ""
+            (count[n] > 1) ? var[n] = M2[n]/(count[n] - 1) : var[n] = "nan"
         }
     }
 }
 
 END {
-    print "col", "mean", "std_err", "std_dev", "count"
+    print "col", "mean", "std_err", "std_dev", "min", "max", "count"
     for (n=1; n<=max_nf; n++) {
         if (header[n])
-            print header[n], mean[n], sqrt(var[n]/count[n]), sqrt(var[n]), count[n]
+            print header[n], mean[n], sqrt(var[n]/count[n]), sqrt(var[n]), min[n], max[n], count[n]
         else
-            print "col" n, mean[n], sqrt(var[n]/count[n]), sqrt(var[n]), count[n]
+            print "col" n, mean[n], sqrt(var[n]/count[n]), sqrt(var[n]), min[n], max[n], count[n]
     }
 }
 
-- 
cgit v1.2.3