From 2482727a6902e44e6a68236f878f5f9bf7947bd2 Mon Sep 17 00:00:00 2001
From: wukong <wukong@longaeva>
Date: Sun, 17 Jun 2018 16:38:38 -0700
Subject: added covariance matrix calculation

---
 lin_reg1.awk | 11 +++++------
 1 file changed, 5 insertions(+), 6 deletions(-)
 mode change 100755 => 100644 lin_reg1.awk

(limited to 'lin_reg1.awk')

diff --git a/lin_reg1.awk b/lin_reg1.awk
old mode 100755
new mode 100644
index cf3e69a..8b03e2f
--- a/lin_reg1.awk
+++ b/lin_reg1.awk
@@ -4,7 +4,7 @@
 # simple linear regression between columns
 
 BEGIN {
-    OFS = "%.18g"
+    OFMT = "%.18g"
     sign = "[+-]?"
     decimal = "[0-9]+[.]?[0-9]*"
     fraction = "[.][0-9]*"
@@ -13,13 +13,11 @@ BEGIN {
 }
 
 NR == 1 {
-    header_nf = NF
-    for (n=1; n<=NF; n++) {
-        ($n !~ number) ? header[n] = $n : header[n] = "col" n
-    }
+    for (n=1; n<=NF; n++)
+        ($n ~ number) ? header[n] = "col" n : header[n] = $n
 }
 
-NF != 0 {
+NF > 0 {
     if (NF > nf_max)
         nf_max = NF
 
@@ -44,6 +42,7 @@ NF != 0 {
                 sum_xy[x,y] += $x*$y
                 sum_delta_xy[x,y] += delta[x]*delta[y]
 
+                # correlation
                 r_den[x,y] = sqrt(sum_delta2[x]*sum_delta2[y])
                 (r_den[x,y]) ? r[x,y] = sum_delta_xy[x,y]/r_den[x,y] : r[x,y] = 0
 
-- 
cgit v1.2.3