summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorwukong <wukong@longaeva>2019-01-10 23:41:14 -0800
committerwukong <wukong@longaeva>2019-01-10 23:41:14 -0800
commit05661d00ca79b009d96dc677a7157377ae272793 (patch)
treea64ebf2c7103f9aef84e4dce9d1964edda889051
parent6333d2e42921b99bf0af537649d47def0b3f4c32 (diff)
added totals to mean and mean_avg
Diffstat (limited to '')
-rw-r--r--lin_reg.awk2
-rw-r--r--mean.awk4
-rw-r--r--mean_avg.awk6
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])
diff --git a/mean.awk b/mean.awk
index b6a794e..3936902 100644
--- a/mean.awk
+++ b/mean.awk
@@ -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]
}
}