小程序救助，希望大神能帮忙指点，小弟跪谢~~~~

albertkingking

求大神指导~~~跪谢！！！！

现有文件格式如下

>seq1        12
CAGCTAGCTACGTCG
>seq2        22
CACGATGCTACGTAGCAGCT
>seq3        2
CAGCTACGTAGCGAAG
>seq4        4
CAGCTAGCGAGCTACG
>seq5        6
CACGATCGAGCGTAGCTAGC

现在想得到不同长度的序列 他们ID行的后面的数字的和；
好比如下结果：

length=15 number=12
length=20 number=28
length=16 number=6

长度是15的序列，number就是原文件里的12
长度是20的序列，number就是原文件里的22+6=28
长度是16的序列，number就是原文件里的2+4=6


求大神指点，万分感谢~~~~~

huang6894

回复 1# albertkingking

1. my %h;
   
2. $/=">";
   
3. 
   
4. <DATA>;
   
5. while(<DATA>){
   
6.         chomp;
   
7.         my ($name,@seq) = (split /\n/);
   
8.         my $key = join("",@seq);
   
9.         my $value = (split /\s+/,$name)[-1];
   
10.         $h{length($key)} += $value;
    
11. }
    
12. $/="\n";
    
13. foreach my $k (sort keys %h){
    
14.         print "length=".$k."\tnumber=".$h{$k}.$/;
    
15. }
    
16. __DATA__
    
17. >seq1        12
    
18. CAGCTAGCTACGTCG
    
19. >seq2        22
    
20. CACGATGCTACGTAGCAGCT
    
21. >seq3        2
    
22. CAGCTACGTAGCGAAG
    
23. >seq4        4
    
24. CAGCTAGCGAGCTACG
    
25. >seq5        6
    
26. CACGATCGAGCGTAGCTAGC

复制代码

albertkingking

好厉害 我咋就想不出来呢 谢谢~~~~~回复 2# huang6894


   

yestreenstars

1. #!/usr/bin/perl
   
2. use strict;
   
3. use warnings;
   
4. 
   
5. my($num, %hash);
   
6. while(<DATA>){
   
7.         chomp;
   
8.         if($. % 2){
   
9.                 $num = (split)[1];
   
10.         }else{
    
11.                 $hash{length $_} += $num;
    
12.         }
    
13. }
    
14. printf "length=%d number=%d\n", $_, $hash{$_} for keys %hash;
    
15. 
    
16. __DATA__
    
17. >seq1        12
    
18. CAGCTAGCTACGTCG
    
19. >seq2        22
    
20. CACGATGCTACGTAGCAGCT
    
21. >seq3        2
    
22. CAGCTACGTAGCGAAG
    
23. >seq4        4
    
24. CAGCTAGCGAGCTACG
    
25. >seq5        6
    
26. CACGATCGAGCGTAGCTAGC

复制代码

清泉一边

我也写了一个
#!/usr/bin/perl -w
$/=">";<DATA>;$/="\n";
my %hash;
while (<DATA>){
my ($a,$number) = split /\s+/;
$/=">";
my $seq = <DATA>;
  $/="\n";
$seq =~ s/>$//;
$seq =~ s/\n//g;
my $length = length($seq);
push @{$hash{$length}},$number;
}
foreach my $k (sort keys %hash){
my $sum = 0;
    $sum += $_ for @{$hash{$k}};
print "length=$k\tnumber=$sum\n";
}

__DATA__
>seq1        12
CAGCTAGCTACGTCG
>seq2        22
CACGATGCTACGTAGCAGCT
>seq3        2
CAGCTACGTAGCGAAG
>seq4        4
CAGCTAGCGAGCTACG
>seq5        6
CACGATCGAGCGTAGCTAGC回复 1# albertkingking


   

albertkingking

醉了  太牛叉了  谢谢~~~~回复 4# yestreenstars


   

albertkingking

都好牛~~~~ 学习了~~~~回复 5# 清泉一边