Awk Notes

From Federal Burro of Information
Revision as of 15:39, 2 November 2022 by David (talk | contribs)
Jump to navigationJump to search

I love you awk

standard deviation with awk:

   awk '{for(i=1;i<=NF;i++) {sum[i] += $i; sumsq[i] += ($i)^2}} 
         END {for (i=1;i<=NF;i++) {
         printf "%f %f \n", sum[i]/NR, sqrt((sumsq[i]-sum[i]^2/NR)/NR)}
        }' file.dat >> aver-std.dat


```

   awk '{for(i=1;i<=NF;i++) {sum[i] += $i; sumsq[i] += ($i)^2}} 
         END {for (i=1;i<=NF;i++) {
         printf "%f %f \n", sum[i]/NR, sqrt((sumsq[i]-sum[i]^2/NR)/NR)}
        }' file.dat >> aver-std.dat

```