求个perl程序，生物信息方面的！

conn2011

1. #!/usr/bin/perl -w
   
2. use strict;
   
3. use Data::Dumper;
   
4. 
   
5. my %in;
   
6. open(IN, "file2path") || die "Cannot open file int.txt\n";
   
7. while(<IN>){
   
8.             chomp;
   
9.             my @line = split(' ',$_);
   
10.             $in{$line[1]}{'text'} = $line[3];
    
11.             push @{$in{$line[1]}{'value'}},$line[5];
    
12.         }
    
13. close IN;
    
14. 
    
15. my %info;
    
16. open(INFO,"file1path") || die "Cannot open file info.txt\n";
    
17. while(<INFO>){
    
18.             chomp;
    
19.             my @line = split(' ',$_);
    
20.             $info{$line[0]} = $line[2];
    
21.         }
    
22. close INFO;
    
23. 
    
24. foreach my $key (keys %in){
    
25.             my $sum = 0;
    
26.             foreach my $value (@{$in{$key}{'value'}}){
    
27.                         $sum += $value;
    
28.                     }
    
29.             $in{$key}{'avg'} = $sum/scalar(@{$in{$key}{'value'}});
    
30.             $in{$key}{'time'} = scalar(@{$in{$key}{'value'}});
    
31.             my $tot = $info{$key} ? $info{$key} : 0;
    
32.             my $margin = $tot - $in{$key}{'avg'};
    
33.             printf "+ %s\t%s\t%10s\t%10s\t%s\n",($key,$in{$key}{'text'},$in{$key}{'avg'},$margin,$in{$key}{'time'});
    
34.         }

复制代码