生物信息，长度分布统计，求大神指导一下！！！

旋转小马

>a1;30266 K: 25 length: 7764
>a2;21197 K: 25 length: 1303
>a3;19639 K: 25 length: 6523
>a4;21559 K: 25 length: 1051
>a5;9715 K: 25 length: 815
>a6;60114 K: 25 length: 49
>a7;9601 K: 25 length: 1387
>a8;11933 K: 25 length: 1662
>a9;28016 K: 25 length: 96
>a10;5612 K: 25 length: 1879
>a11;12713 K: 25 length: 916
>a12;13947 K: 25 length: 442
>a13;5887 K: 25 length: 2657
>a14;2346 K: 25 length: 9605
>a15;9119 K: 25 length: 405
>a16;2507 K: 25 length: 4692
>a17;18192 K: 25 length: 100
>a18;8527 K: 25 length: 1909
>a19;2477 K: 25 length: 6703
>a20;12117 K: 25 length: 949
>a21;27978 K: 25 length: 49
>a22;10051 K: 25 length: 1189
>a23;6358 K: 25 length: 470
>a24;5987 K: 25 length: 3481
>a25;15254 K: 25 length: 94
>a26;2738 K: 25 length: 1960
>a27;6649 K: 25 length: 603
>a28;11581 K: 25 length: 543
>a29;4613 K: 25 length: 1895
>a30;6662 K: 25 length: 811
>a31;22138 K: 25 length: 35
>a32;8808 K: 25 length: 372 现在我需要做的是统计最后一列的分布个数，以10为一单位，
比如说，最后输出一文件，是按照0~10之间出现的个数，11~20之间出现的个数、、、、、、、、、以此类推，最后输出一个文档！！！

旋转小马

万分感谢，又学到了好多东西！！！回复 17# 飞越地平线


   

飞越地平线

回复 16# 旋转小马

while书上有介绍。我只知道它是一个循环结构。在这里是一行一行读入。
你多看看小骆驼书，还有楼上几位写的代码。

1. #!/usr/bin/perl -w
   
2. use strict;
   
3. my %hash;
   
4. my $x;
   
5. my $y;
   
6. #my $number;
   
7. while(<DATA>){
   
8.         chomp;
   
9.         my @line=split/\s+/,$_;
   
10.         #$number=int($line[-1]/10);
    
11.         if(exists $hash{int($line[-1]/10)}){
    
12.            $hash{int($line[-1]/10)}++;
    
13.         }else{
    
14.            $hash{int($line[-1]/10)}=1;
    
15.         }
    
16. 
    
17. }
    
18. map{
    
19.         $x=$_*10;
    
20.         $y=$x+9;
    
21.         print "$x~$y\t".$hash{$_}."\n";
    
22. } sort {$a<=>$b} keys %hash;
    
23. 
    
24. __DATA__
    
25. asdf 123
    
26. dfasdf 456
    
27. adsfads 457
    
28. asd 12
    
29. aefd 98
    
30. asdf 1
    
31. asdf 2
    
32. asdf 9
    
33. asdf 10

复制代码

旋转小马

原来是这样，while我还不大理解什么意思呢，继续学习，再次感谢，我想让区间从小到大排序如何实现呢？回复 15# 飞越地平线


   

飞越地平线

回复 14# 旋转小马

那个要统计完个数存在%hash之后，然后打印出结果来的啊

1. #!/usr/bin/perl -w
   
2. use strict;
   
3. my %hash;
   
4. my $a;
   
5. my $b;
   
6. my $c;
   
7. #my $number;
   
8. while(<DATA>){
   
9.         chomp;
   
10.         my @line=split/\s+/,$_;
    
11.         #$number=int($line[-1]/10);
    
12.         if(exists $hash{int($line[-1]/10)}){
    
13.            $hash{int($line[-1]/10)}++;
    
14.         }else{
    
15.            $hash{int($line[-1]/10)}=1;
    
16.         }
    
17. 
    
18. }
    
19. map{
    
20.         $a=$_;
    
21.         $b=$a*10;
    
22.         $c=$b+9;
    
23.         print "$b~$c\t".$hash{$_}."\n";
    
24. } sort keys %hash;
    
25. 
    
26. __DATA__
    
27. asdf 123
    
28. dfasdf 456
    
29. adsfads 457
    
30. asd 12
    
31. aefd 98
    
32. asdf 1
    
33. asdf 2
    
34. asdf 9
    
35. asdf 10

复制代码

旋转小马

  1 #!/usr/bin/perl -w
  2 use strict;
  3 open IN,"len.xls"or die"Can not open IN.\n";
  4 open OUT,">result.xls"or die"Can not open endfile.\n";
  5 my %hash;
  6 my $a;
  7 my $b;
  8 my $c;
  9 #my $number;
10 while(<IN>){
11         chomp;
12         my @line=split/\s+/,$_;
13         #$number=int($line[-1]/10);
14         if(exists $hash{int($line[-1]/10)}){
15            $hash{int($line[-1]/10)}++;
16         }else{
17            $hash{int($line[-1]/10)}=1;
18         }
19         $a=$hash{int($line[-1]/10)};
20         $b=$a*10;
21         $c=$b+9;
22         print "$b~$c\t$a\n";
23 }
24 close IN;
25 close OUT;这是我写的完整代码，最后无法实现呢，能不能帮我看看，是哪里出错了呢？回复 9# 飞越地平线


   

旋转小马

嗯嗯，好的，明天到公司试一试，这样的更简单，谢啦~~~回复 12# rdcwayx


   

rdcwayx

by awk

1. awk '{a[int($NF/10)]++}END{for (i in a) print i*10 "-" i*10+9,a[i]}' infile |sort -n
   
2. 30-39 1
   
3. 40-49 2
   
4. 90-99 2
   
5. 100-109 1
   
6. 370-379 1
   
7. 400-409 1
   
8. 440-449 1
   
9. 470-479 1
   
10. 540-549 1
    
11. 600-609 1
    
12. 810-819 2
    
13. 910-919 1
    
14. 940-949 1
    
15. 1050-1059 1
    
16. 1180-1189 1
    
17. 1300-1309 1
    
18. 1380-1389 1
    
19. 1660-1669 1
    
20. 1870-1879 1
    
21. 1890-1899 1
    
22. 1900-1909 1
    
23. 1960-1969 1
    
24. 2650-2659 1
    
25. 3480-3489 1
    
26. 4690-4699 1
    
27. 6520-6529 1
    
28. 6700-6709 1
    
29. 7760-7769 1
    
30. 9600-9609 1
    

复制代码

rubyish

BiRu:

1. #!/usr/bin/perl
   
2. use strict;
   
3. use warnings;
   
4. 
   
5. my %range;
   
6. 
   
7. while (<DATA>) {
   
8.     my @line = split;
   
9.     my $number = int( ( $line[-1] - 1 ) / 10 );
   
10. 
    
11.     # buhao:    my $number = int( $line[-1] / 10 );
    
12.     # biru:     11 => int 11/10 = 1   OK
    
13.     # danshi:   10 => int 10/10 = 1   NG
    
14.     # yinwei:   10 => range 0, 11 => range 1
    
15.     # suoyi:    10 => int( ( 10 - 1 ) / 10 ) = 0  OK
    
16.     # yiji:     11 => int( ( 11 - 1 ) / 10 ) = 1  OK
    
17.     # yiji:     0  => int( ( 0 - 1 ) / 10 )  = 0  OK
    
18. 
    
19.     $range{$number}++;
    
20. }
    
21. 
    
22. # xianzai %range: { 1 => 2, 4 => 1, 0 => 4 }
    
23. 
    
24. my $max = 0;
    
25. for ( keys %range ) {
    
26.     $max = $_ if $max < $_;
    
27. }
    
28. 
    
29. # print "max: $max\n";
    
30. # max: 4
    
31. 
    
32. for ( 0 .. $max ) {
    
33.     my $begin = $_ ? ( $_ * 10 ) + 1 : 0;
    
34. 
    
35.     # buhao:    my $begin = $_ * 10 + 1;
    
36.     # biru:     1 => 1 * 10 + 1 = 11   OK: 11-
    
37.     # danshi:   0 => 0 * 10 + 1 = 1    NG: 1-, OK: 0-
    
38. 
    
39.     my $end   = ( $_ + 1 ) * 10;
    
40. 
    
41.     # buhao:    my $end = $begin + 9;
    
42.     # biru:     11 + 9 = 20    OK: 11-20
    
43.     # danshi:   0 + 9 = 9      NG: 0-9, OK: 0-10
    
44.     # huo:      my $end = ( $_ + 1 ) . '0';
    
45.     # huo:      my @end = ( $_ + 1, 0 );
    
46.     # huo:      my $end = $begin + $_ ? 9 : 10;
    
47. 
    
48.     my $count = $range{$_} // 0;
    
49.    
    
50.     # buhao:    my $count = $range{$_};
    
51.     # biru:     21-30
    
52.     # Use of uninitialized value $count in print....
    
53.     # 11-20 2
    
54.     # 21-30 HERE EMPTY
    
55.    
    
56.     print $begin, '-', $end, "\t", $count, "\n";
    
57. }
    
58. 
    
59. __DATA__
    
60. a 1
    
61. b 2
    
62. c 0
    
63. d 10
    
64. e 11
    
65. f 12
    
66. g 44

复制代码

1. 0-10    4
   
2. 11-20   2
   
3. 21-30   0
   
4. 31-40   0
   
5. 41-50   1

复制代码

旋转小马

嗯哪，明天在公司试试，住的地方电脑不好弄，不会再请教！！！回复 9# 飞越地平线