shell 求某一列的中值，标准差，最大值，最小值

pdzcowboy

我学shell快一个月了，主要想用来处理一些文本，确实挺方便的。但我用shell来计算中值，标准差，最大值，最小值这些常见的统计量，发现挺不方便的。
像我如果要计算第二列的这些统计量，有什么简单的方法吗。
cat data|awk 'BEGIN {max = 0} {if ($2>max) max=$2 fi} END {print "Max=", max}'[/code]
这个令从网上抄下来的，这个“fi”是什么意思？其他几个统计量怎么求？有内置的函数之类吗，有可能问题比较白痴，希望大家能谅解，我也一直在恶补相关知识，谢谢了。

162616.5347586470       87367.4899084031
85181.9193508478        115419.9635060868
128583.9942762910       102857.0688335032
171259.8270212160       130301.9875994774
148783.9678689920       128483.8809052336
93168.6444210467        129336.0618793606
110496.2277431580       123883.6186766724
174233.2178257970       150691.5762701073
136478.9755860420       157927.8243340223
161632.3055155560       121767.1137681215
133953.8764284690       176559.0927623106
168696.1419453560       185411.2571383122
129933.4663808970       148160.8502728328
283201.8596277380       205012.4679922750
147516.3518521880       197898.2265649752

Herowinter

有现成的函数吗？没有的话，标准差这种
只能通过平均数自己算了，用Shell做这种
复杂的运算，合适吗？

jason680

回复 1# pdzcowboy

$ awk '{v[NR]=$2;sum+=$2}{if($2>max)max=$2;if(min>$2)min=$2}NR==1{max=min=$2}END{avg=sum/NR;for(n=1;n<=NR;n++)sd+=(v[n]-avg)**2;sd=sqrt(sd/NR);printf("Max=%f\nMin=%f\nSum=%f\nAvg=%f\nSD=%f\nItem=%d\n",max,min,sum,avg,sd,NR)}' data
Max=205012.467992
Min=87367.489908
Sum=2161078.480412
Avg=144071.898694
SD=33583.438884
Item=15


   
Note: fi for shell script with if function, not for awk

这个冬天不冷

回复 1# pdzcowboy


   这么复杂的统计 还是找个脚本 语言处理，python  ，都很好，

另外你说fi是 shell中 有这种语法 if fi    成对出现 ，至于awk的 if没见过这么写的，貌似也没别的啥意思

wslm31

awk '{sum=sum+$2}END{printf sum/NR "\n"}' date
显示平均值
cat date|awk '{printf $2"\n"}'|sort|sed -n '1p'
显示最大值，排序一下，取第一行值
cat date|awk '{printf $2"\n"}'|sort|sed -n '$p'
显示最小值，先排序一下，取最尾值

li0924

至于标准差忘记是什么了。这里求最大值和最小值

1. awk 'NR==1{max=min=$2;next}{max=max>$2?max:$2;min=min<$2?min:$2}END{print max,min}'

复制代码

你其中的fi；那个是shell的语法；跟awk的语法弄混淆了。

yestreenstars

标准差这些这么专业的就应该找专业的工具处理，awk没有这类函数~

这个冬天不冷

至于呢个fi经过大家讨论，是一个 无效的 变量，对于程序无影响 ，可以换成 其它字符串

rdcwayx

@jason680@pdzcowboy

刚刚用 ideone 这个在线工具测了一下，这个是结果

http://ideone.com/o320AX

Herowinter

回复 6# li0924

如果我没记错的话，标准差应该这样的。
假设有3个数据a1, a2, a3, 平均值avg
标准差=sqrt( ((a1-avg)^2+(a2-avg)^2+(a3-avg)^2)/3)
用来表示样本离散程度的一个参数。

不过，我也不大确定了，本科做物理实验计算
误差是用得比较多；但那个时候我们算出来数据
误差过大就回头改数据，想想真是欢乐。