diff options
Diffstat (limited to '')
| -rw-r--r-- | cov.awk | 1 | ||||
| -rw-r--r-- | lin_reg.awk | 6 |
2 files changed, 3 insertions, 4 deletions
@@ -14,7 +14,6 @@ BEGIN { NR == 1 { for (y=1; y<=NF; y++) ($y ~ number) ? header[y] = "col" y : header[y] = $y - printf(header[y]) } NF > 0 { 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]) } } } |