diff options
Diffstat (limited to 'lin_reg.awk')
| -rw-r--r-- | lin_reg.awk | 14 |
1 files changed, 7 insertions, 7 deletions
diff --git a/lin_reg.awk b/lin_reg.awk index 0213c7f..5a9df56 100644 --- a/lin_reg.awk +++ b/lin_reg.awk @@ -79,10 +79,10 @@ NF { a_err[x,y] = sqrt(sum2[x]/count[x,y])*b_err[x,y] # 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]) - yw_dist[x,y] = b[x,y]*(xw[x,y] - mean[x]) + # 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]) + # yw_dist[x,y] = b[x,y]*(xw[x,y] - mean[x]) } } } @@ -99,11 +99,11 @@ END { (r[x,y]*r[x,y]), header[y], b[x,y], b_err[x,y], header[x], a[x,y], a_err[x,y]) - printf("[" OFMT "," OFMT "][" OFMT "," OFMT "]" OFS "[" OFMT "," OFMT "]" OFS, + printf("[" OFMT "," OFMT "][" OFMT "," OFMT "]" OFS "[" OFMT "," OFMT "]" ORS, 0, a[x,y], (-1.0*a[x,y]/b[x,y]), 0, mean[x], b[x,y]*(mean[x]) + a[x,y]) - printf("[" OFMT "," OFMT "]" OFS, xw[x,y], yw[x,y]) - printf("[" OFMT "]" ORS, sqrt(xw_dist[x,y]*xw_dist[x,y] + yw_dist[x,y]*yw_dist[x,y])) + # printf("[" OFMT "," OFMT "]" OFS, xw[x,y], yw[x,y]) + # printf("[" OFMT "]" ORS, sqrt(xw_dist[x,y]*xw_dist[x,y] + yw_dist[x,y]*yw_dist[x,y])) } } } |