diff options
Diffstat (limited to 'lin_reg1.awk')
| -rw-r--r-- | lin_reg1.awk | 12 |
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]) { |