Awk Notes: Difference between revisions
From Federal Burro of Information
Jump to navigationJump to search
No edit summary |
No edit summary |
||
| Line 3: | Line 3: | ||
standard deviation with awk: | standard deviation with awk: | ||
< | <pre> | ||
awk '{for(i=1;i<=NF;i++) {sum[i] += $i; sumsq[i] += ($i)^2}} | awk '{for(i=1;i<=NF;i++) {sum[i] += $i; sumsq[i] += ($i)^2}} | ||
END {for (i=1;i<=NF;i++) { | END {for (i=1;i<=NF;i++) { | ||
printf "%f %f \n", sum[i]/NR, sqrt((sumsq[i]-sum[i]^2/NR)/NR)} | printf "%f %f \n", sum[i]/NR, sqrt((sumsq[i]-sum[i]^2/NR)/NR)} | ||
}' file.dat >> aver-std.dat | }' file.dat >> aver-std.dat | ||
</ | </pre> | ||
where file.dat is one column of numbers. | |||
[[Category:Computers]] | [[Category:Computers]] | ||
Revision as of 15:40, 2 November 2022
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
where file.dat is one column of numbers.
