菜鸟求助用脚本求标准偏差，谢谢！

tengfei0311

我现在有很多个文件，要求第五列值的标准偏差，用shell脚本该怎么实现，谢谢！

假如其中一个文件名为10.txt，里面部分内容如下：
                                           （5列）
98.7    27.7    102.7   25.1    3.096   0.04    1
98.7    27.7    101.9   25.7    3.038   0.04    1
100.2   26.5    102.7   25.1    3.067   0.04    1
100.7   23.5    98.5    24.7    3.077   0.04    1
98.7    27.7    102.7   25.1    3.129   0.04    1
98.7    24.9    104.2   28.1    3.149   0.04    1
98.5    25.2    104.2   28.1    2.931   0.04    1
98.9    25.8    104.2   28.1    3.023   0.04    1
102.8   23.7    103.6   25.4    3.106   0.04    1
103.6   25.4    101.0   22.8    2.89    0.04    1
103.6   25.4    99.7    27.8    3.055   0.04    1
101.0   22.8    104.2   28.1    3.19    0.04    1
101.0   22.8    104.2   28.1    3.011   0.04    1
102.8   24.1    104.2   28.1    3.152   0.04    1
104.2   28.1    102.0   23.6    3.007   0.04    1
98.1    24.1    103.7   27.3    3.021   0.04    1
98.6    24.4    103.7   27.3    2.99    0.04    1
98.6    24.4    102.9   26.9    2.962   0.04    1
98.6    24.4    103.7   27.3    2.998   0.04    1
100.2   26.9    104.2   28.1    3.049   0.04    1
99.9    26.1    104.2   28.1    3.193   0.04    1
99.9    22.6    104.3   24.9    3.101   0.04    1

jason680

回复 1# tengfei0311

$ awk -vt=5 '{s+=$t;a[NR]=$t}END{N=NR;X=s/N;for(n=1;n<=N;++n){S=S+(a[n]-X)**2}print sqrt(S/N)}' 10.txt
0.0778859


$ awk -vt=5 '{s+=$t;a[NR]=$t}END{N=NR;X=s/N;for(n=1;n<=N;++n){S=S+(a[n]-X)**2}print sqrt(S/(N-1))}' 10.txt
0.0797187

   

jason680

回复 3# 时时20160307

http://baike.baidu.com/view/663819.htm
   
标准偏差(Std Dev,Standard Deviation) -统计学名词。...

tengfei0311

回复 2# jason680

谢谢！

是对的，我用execl算了一遍。