[Bio] / FigKernelScripts / get_coupling_values.pl Repository:
ViewVC logotype

Diff of /FigKernelScripts/get_coupling_values.pl

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.1, Thu Aug 12 16:15:05 2010 UTC revision 1.3, Tue Sep 14 18:39:54 2010 UTC
# Line 1  Line 1 
1  use strict;  use strict;
2  use SeedEnv;  use SeedEnv;
3  use Data::Dumper;  use Data::Dumper;
4    use FIG_Config;
5    
6  my $usage = "get_coupling_values families.2c  [URL]";  my $usage = "get_coupling_values families.2c  [URL]";
7  my($fam2c);  my($fam2c);
# Line 9  Line 10 
10  )  )
11      || die $usage;      || die $usage;
12    
13    my $genome_sets = "$FIG_Config::global/genome.sets";
14    
15  my $url = shift @ARGV;  my $url = shift @ARGV;
16  my $sapO = new SAPserver(  url => $url );  my $sapO = new SAPserver(  url => $url );
17    
# Line 19  Line 22 
22    
23  # my $otuMaps    = $sapO->representative_genomes;  # my $otuMaps    = $sapO->representative_genomes;
24  # my $to_otu     = $otuMaps->[0];  # my $to_otu     = $otuMaps->[0];
25  my $to_otu = { map { ($_ =~ /^(\d+)\t(\d+\.\d+)/) ? ($2 => $1) : () } `cat genome.sets` };  
26    open(GS, "<", $genome_sets) or die "Cannot open $genome_sets: $!";
27    my $to_otu = { map { ($_ =~ /^(\d+)\t(\d+\.\d+)/) ? ($2 => $1) : () } <GS> };
28    close(GS);
29    
30  my %otus = map { $to_otu->{$_} => 1 } keys(%$to_otu);  my %otus = map { $to_otu->{$_} => 1 } keys(%$to_otu);
31    
32  $_ = keys(%otus);  $_ = keys(%otus);
# Line 151  Line 158 
158      return abs($max - $min);      return abs($max - $min);
159  }  }
160    
 sub process_pair {  
     my($sapO,$peg1,$peg2,$to_ff) = @_;  
   
 }  

Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3