统计数据：请教

huazhorg

因为经常要处理类似于这样的文件，要统计这个文件中第二列与第三列都为0时，第一列对应的区间长度。
例如一个名为data的文件，cat data:
  0.96  0.171  0.000           
  0.97  0.111  0.000           
  0.98  0.127  0.000           
  0.99  0.271  0.000           
  1.00  0.000  0.000           
  1.01  0.000  0.000            
  1.34  0.000  0.000            
  1.35  0.000  0.000            
  1.36  0.000  0.000            
  1.37  0.000  0.000            
  1.38  0.000  0.312            
  1.39  0.268  0.154            
  1.40  0.550  0.228            
  1.41  0.781  0.000               
  1.42  0.130  0.000           
  1.96  0.000  0.000            
  1.97  0.000  0.000            
  2.00  0.000  0.000           
  2.01  0.000  0.000            
  2.02  0.053  0.000            
  2.03  0.188  0.110            
  2.04  0.849  0.760
第二列与第三列都为0时，对应的第一列区间为：  1.00～1.37，1.96～2.01
所以对应的第一列的区间长度为1.37-0.99+2.01-1.95=0.38+0.06=0.44
请问，如何用命令操作这个文件，得到这个结果呢？谢谢！

qinggzha

给你提供思路：

如果第二三列都为0是，打印出 行数a 和此时第一列的数字b
然后a++，判断，如果a++ 行数的第二三列还未0，继续a++。 知道a++行数时，二三列中一个不为0时候，打印出行数 c ，和此时的第一列数字d ，
然后 d- b-0.01

seesea2517

1. [seesea2517@UC ~] #awk -f stat.awk data
   
2. 0.44
   
3. 
   
4. [seesea2517@UC ~] #cat stat.awk
   
5. BEGIN {
   
6.     reset()
   
7.     result = 0
   
8. }
   
9. 
   
10. {
    
11.     if ($2==0 && $3==0)
    
12.     {
    
13.         min = ($1 < min ? $1 : min)
    
14.         max = ($1 > max ? $1 : max)
    
15. 
    
16.         next
    
17.     }
    
18. 
    
19.     if (max < 0)
    
20.     {
    
21.         next
    
22.     }
    
23. 
    
24.     result += max - min + 0.01
    
25.     reset()
    
26. }
    
27. 
    
28. END {
    
29.     print result
    
30. }
    
31. 
    
32. function reset()
    
33. {
    
34.     min = 9999
    
35.     max = -1
    
36. }
    

复制代码

huazhorg

seesea2517 发表于 2012-11-14 11:12

首先，谢谢，这个脚本没有问题，就是稍显长了点，不知能否更精简点？

huazhorg

awk 'BEGIN{c=0} {if($2==0 && $3==0){c=$1;b=[ ]; result=c-b+0.01 }} END{print result}' data

其中的b=[ ] 我想表达为第一个$2==0 && $3==0时的$1,怎么表达？

huazhorg

qinggzha 发表于 2012-11-14 10:25
给你提供思路：

如果第二三列都为0是，打印出 行数a 和此时第一列的数字b

思路比较清晰，不过我想知道，一定要打印行数吗？有什么作用呢

qinggzha

没必要行数，我只是用来做递增的时候根据行数来统计，所以打印了 回复 6# huazhorg


   

huazhorg

qinggzha 发表于 2012-11-14 16:14
没必要行数，我只是用来做递增的时候根据行数来统计，所以打印了 回复 6# huazhorg

能否把你写的贴出来看看？

seesea2517

huazhorg 发表于 2012-11-14 15:18
首先，谢谢，这个脚本没有问题，就是稍显长了点，不知能否更精简点？

本人无能，不能精简。

huazhorg

seesea2517 发表于 2012-11-14 16:34
本人无能，不能精简。

感谢！那就这样了。