diff options
| author | wukong <wukong@longaeva> | 2019-01-10 23:41:14 -0800 |
|---|---|---|
| committer | wukong <wukong@longaeva> | 2019-01-10 23:41:14 -0800 |
| commit | 05661d00ca79b009d96dc677a7157377ae272793 (patch) | |
| tree | a64ebf2c7103f9aef84e4dce9d1964edda889051 | |
| parent | 6333d2e42921b99bf0af537649d47def0b3f4c32 (diff) | |
added totals to mean and mean_avg
Diffstat (limited to '')
| -rw-r--r-- | lin_reg.awk | 2 | ||||
| -rw-r--r-- | mean.awk | 4 | ||||
| -rw-r--r-- | mean_avg.awk | 6 |
3 files changed, 6 insertions, 6 deletions
diff --git a/lin_reg.awk b/lin_reg.awk index 7e71804..aeb5dfc 100644 --- a/lin_reg.awk +++ b/lin_reg.awk @@ -75,7 +75,7 @@ NF > 0 { if (a_err_den[x,y]) a_err[x,y] = sqrt(sum2[x]/count[x,y])*b_err[x,y] - # weighted mean, from HP-20S manual + # weighted mean, from HP-20S manual, pg 60 xw[x,y] = sum_xy[x,y]/sum[y] yw[x,y] = b[x,y]*xw[x,y] + a[x,y] xw_dist[x,y] = (xw[x,y] - mean[x]) @@ -29,7 +29,7 @@ NF > 0 { }
END {
- print "mean", "std_err", "std_dev", "min", "max", "count"
- print mean, sqrt(var/count), sqrt(var), min, max, count
+ print "mean", "std_err", "std_dev", "min", "max", "total", "count"
+ print mean, sqrt(var/count), sqrt(var), min, max, (mean*count), count
}
diff --git a/mean_avg.awk b/mean_avg.awk index 4bf7051..dc39118 100644 --- a/mean_avg.awk +++ b/mean_avg.awk @@ -36,12 +36,12 @@ NF > 0 { } END { - print "col", "mean", "std_err", "std_dev", "min", "max", "count" + print "col", "mean", "std_err", "std_dev", "min", "max", "total", "count" for (n=1; n<=max_nf; n++) { if (header[n]) - print header[n], mean[n], sqrt(var[n]/count[n]), sqrt(var[n]), min[n], max[n], count[n] + print header[n], mean[n], sqrt(var[n]/count[n]), sqrt(var[n]), min[n], max[n], mean[n]*count[n], count[n] else - print "col" n, mean[n], sqrt(var[n]/count[n]), sqrt(var[n]), min[n], max[n], count[n] + print "col" n, mean[n], sqrt(var[n]/count[n]), sqrt(var[n]), min[n], max[n], mean[n]*count[n], count[n] } } |
