diff options
| author | wukong <wukong@longaeva> | 2019-01-06 00:58:38 -0800 |
|---|---|---|
| committer | wukong <wukong@longaeva> | 2019-01-06 00:58:38 -0800 |
| commit | 37e92fab3085b52ff76efe00b9c95731d9da8294 (patch) | |
| tree | 71b1c3482d989b716744013fca7b31681732ca63 /lin_reg.awk | |
| parent | ddcb1f73423a9b47d9b457cf3b19f622b0aa567f (diff) | |
added weighted mean estimate, avg x per y, to lin_reg output
Diffstat (limited to 'lin_reg.awk')
| -rw-r--r-- | lin_reg.awk | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/lin_reg.awk b/lin_reg.awk index aa94b0a..b81ff0f 100644 --- a/lin_reg.awk +++ b/lin_reg.awk @@ -48,7 +48,6 @@ NF > 0 { # covariance #(count[x,y] > 1) ? cov[x,y] = sum_delta_xy[x,y]/(count[x,y] - 1) : cov[x,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] = 1 @@ -86,10 +85,11 @@ END { for (y=1; y<=nf_max; y++) { for (x=1; x<=nf_max; x++) { if (x != y && r[x,y]) { - printf(OFMT OFS "(%s)" OFS " = (" OFMT " +/- " OFMT ")(%s)" OFS " + (" OFMT " +/- " OFMT ")" ORS, + printf(OFMT OFS "(%s)" OFS " = (" OFMT " +/- " OFMT ")(%s)" OFS " + (" OFMT " +/- " OFMT ") \t [" OFMT " (%s/%s)]" ORS, (r[x,y]*r[x,y]), header[y], b[x,y], b_err[x,y], - header[x], a[x,y], a_err[x,y]) + header[x], a[x,y], a_err[x,y], + sum_xy[x,y]/sum[y], header[x], header[y]) } } } |