循环不执行

明瞳狐狸

回复 8# flw
不好意思啊，我正在学，不过有些地方我自己实在想不明白，又找不到人问，才上来发帖的，可能我的问题比较幼稚，让你见笑了


   

Okelani

hg19.txt

1. 1        104
   
2. 2        225
   

复制代码

obj.txt

1. 1        29
   
2. 1        78
   
3. 1        13
   
4. 1        67
   
5. 1        49
   
6. 1        10
   
7. 1        32
   
8. 1        71
   
9. 1        52
   
10. 1        78
    
11. 1        14
    
12. 1        36
    
13. 1        60
    
14. 1        95
    
15. 1        42
    
16. 1        67
    
17. 1        20
    
18. 1        0
    
19. 1        97
    
20. 1        10
    
21. 2        220
    
22. 2        166
    
23. 2        54
    
24. 2        187
    
25. 2        92
    
26. 2        115
    
27. 2        132
    
28. 2        172
    
29. 2        16
    
30. 2        133
    
31. 2        68
    
32. 2        125
    
33. 2        159
    
34. 2        56
    
35. 2        192
    
36. 2        14
    
37. 2        75
    
38. 2        89
    
39. 2        87
    
40. 2        206
    

复制代码

newout.txt

1. 1        1        3
   
2. 1        2        3
   
3. 1        3        1
   
4. 1        4        2
   
5. 1        5        3
   
6. 1        6        1
   
7. 1        7        3
   
8. 1        8        2
   
9. 1        9        0
   
10. 1        10        2
    
11. 2        1        2
    
12. 2        2        0
    
13. 2        3        2
    
14. 2        4        4
    
15. 2        5        1
    
16. 2        6        4
    
17. 2        7        0
    
18. 2        8        3
    
19. 2        9        2
    
20. 2        10        2
    

复制代码

code

1. #!/usr/bin/perl
   
2. 
   
3. open my $A, 'hg19.txt';
   
4. open my $B, 'obj.txt';
   
5. open my $C, '>', 'newout.txt';
   
6. my $N = 10;
   
7. my %pos;
   
8. 
   
9. while (<$B>) {
   
10.     my ( $chr, $pos ) = split;
    
11.     push @{ $pos{$chr} }, $pos;
    
12. }
    
13. 
    
14. while (<$A>) {
    
15.     my ( $chr, $len ) = split;
    
16.     my $step = $len / $N;
    
17.     my @steps = map { $step * $_ } 1 .. $N;
    
18.     $steps[-1] = $len;
    
19.     my @result = (0) x $N;
    
20.     my $begin  = 0;
    
21. 
    
22.     for my $P ( sort { $a <=> $b } @{ $pos{$chr} } ) {
    
23.         for my $i ( $begin .. $#steps ) {
    
24.             $P <= $steps[$i] ? do { $result[$i]++; last } : $begin++;
    
25.         }
    
26.     }
    
27. 
    
28.     my $zone = 1;
    
29.     print $C "$chr\t", $zone++, "\t$_\n" for @result;
    
30. }
    

复制代码

明瞳狐狸

回复 12# Okelani
非常感谢你~
里边有几个用法我是第一次见到，需要好好琢磨一下。