如何使用 awk 打印每个 ID 的测量值的标准差？

Question

我需要打印每个唯一 ID (1 美元) 的测量标准差 (2 美元)。

数据如下所示：

101 560  
101 460
101 530  
101 480  
104 600
104 510
104 500
107 450
107 490
107 550
107 500

Answer 1

R 或 datamash 可能是更好的选择！

根据标准差的定义，我们可以：

$ cat my-sd
#!/usr/bin/awk -f

    { s[$1]["sum"] += $2 ; 
      n = s[$1]["oco"] ++; 
      v[$1][n]=$2  }

END { for(x in s){ 
         m=s[x]["sum"]/s[x]["oco"]; 
         s1=0;
         for(y in v[x]){
            s1 += (v[x][y]-m)*(v[x][y]-m);}
         print x, sqrt(s1/s[x]["oco"])}
    }

$ my-sd example
101 39.6074
104 44.9691
107 35.6195