summaryrefslogtreecommitdiff
path: root/lin_reg1.awk
diff options
context:
space:
mode:
authorwukong <wukong@longaeva>2018-06-13 21:41:39 -0700
committerwukong <wukong@longaeva>2018-06-13 21:41:39 -0700
commit8875b420a873a360a9484778b487394add318a5e (patch)
treed357816354c5d9d4bc3b1db5c43ceda7933b5e3f /lin_reg1.awk
parent6bb6be845b6099550e0a530217633e30522093ad (diff)
added shebang, replaced if-elses with shortcut notation, ported online mean and variance calc to mean.awk
Diffstat (limited to 'lin_reg1.awk')
-rw-r--r--lin_reg1.awk12
1 files changed, 4 insertions, 8 deletions
diff --git a/lin_reg1.awk b/lin_reg1.awk
index 866cab8..de1189f 100644
--- a/lin_reg1.awk
+++ b/lin_reg1.awk
@@ -1,3 +1,5 @@
+#!/usr/bin/awk -f
+
### lin_reg1.awk
# simple linear regression between individual text columns
@@ -12,10 +14,7 @@ BEGIN {
NR == 1 {
header_nf = NF
for (n=1; n<=NF; n++) {
- if ($n !~ number)
- header[n] = $n
- else
- header[n] = sprintf("col_%g", n)
+ ($n !~ number) ? header[n] = $n : header[n] = sprintf("col_%g", n)
}
}
@@ -45,10 +44,7 @@ NF != 0 {
sum_delta_xy[x,y] += delta[x]*delta[y]
r_den[x,y] = sqrt(sum_delta2[x]*sum_delta2[y])
- if (r_den[x,y])
- r[x,y] = sum_delta_xy[x,y]/r_den[x,y]
- else
- r[x,y] = 0
+ (r_den[x,y]) ? r[x,y] = sum_delta_xy[x,y]/r_den[x,y] : r[x,y] = 0
ab_den[x,y] = (count[x,y]*sum2[x] - sum[x]*sum[x])
if (ab_den[x,y]) {