求助：刚入门，想请教一个挺简单的问题

CJ妹子

各位大侠 我刚接触PERL 也没学过C 各种不懂呀
我想计算一条序列上各种氨基酸的数量 然后这个写的 很多错误呀
希望大家帮我改改
谢啦 本人在线等

1. #!/usr/bin/perl
   
2. use strict;
   
3. use warnings;
   
4. 
   
5. my $seq = 'ACLKRASFAKSVIGTPEFMAP';
   
6. my %AA;
   
7. foreach  $seq {
   
8.         for( $i=0;$i<length($seq);$i++) {
   
9.                 my $aminoAcid= substr($seq,$i,1);
   
10.                 my $Key=$aminoAcid."$i";
    
11.                 if(exists $AA{$Key}){
    
12.                         $AA{$Key} +=1;
    
13.                 }
    
14.                 else{
    
15.                 $AA{$Key}= 1;
    
16.                 }
    
17.         }
    
18. }
    
19. while (( $key, $value) = each %AA) {
    
20.         print "$key=$value\n";
    
21. }

复制代码

CJ妹子

本人在线等呀

zhlong8

$seq =~ /(.)(?{$AA{$1}++})/g;

%AA 中就包含所有计数了

CJ妹子

回复 3# zhlong8


      能写一下不

CJ妹子

额 上面那个表情居然这么谄媚
其实  我都跪谢了T T

zhlong8

use Data::Dumper;

()= $seq =~ /(.)(?{$AA{$1}++})/g;

print Dumper \%AA

CJ妹子

大师 那个没学过
T T
如果是用一个哈希
key是单词 value是出现次数
每读进一行 判断哈希中是否包含这个key 包含的话对应的value+1 否则将这个单词作为key加到哈希中 value初始化为1
这样做呢？

sellie

楼主能先解释一下怎样 计算一条序列上各种氨基酸的数

sellie

楼主不用解释了. 我看懂你的代码了. 我觉的你的代码写的蛮好的,我一眼看不出来什么问题.
你用哈希判断键值是否存在的思路很正确的.估计有什么语法错误.
这个foreach 大循环应该不需要的
foreach  $seq {

}

CJ妹子

回复 9# sellie


    FOREACH这个循环 我想的是 因为我想统计$seq上所有字母出现的频率 所以就用foreach了