From 70f103736ce9a86055dba424c2390270f4ae83e4 Mon Sep 17 00:00:00 2001 From: wukong Date: Sun, 7 Dec 2025 10:07:44 -0800 Subject: commented out weighted mean evaluation in lin_reg; --- lin_reg.awk | 14 +++++++------- 1 file changed, 7 insertions(+), 7 deletions(-) (limited to 'lin_reg.awk') 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])) } } } -- cgit v1.2.3