diff options
| author | wukong <wukong@longaeva> | 2025-12-07 10:07:44 -0800 |
|---|---|---|
| committer | wukong <wukong@longaeva> | 2025-12-07 10:07:44 -0800 |
| commit | 70f103736ce9a86055dba424c2390270f4ae83e4 (patch) | |
| tree | f2ba51664fb62280178374e959d37187afd4d75f | |
| parent | 800f3cc728607b4d097eb05a710a5f44620ced34 (diff) | |
commented out weighted mean evaluation in lin_reg;
| -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])) } } } |