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

View of /FigKernelScripts/server_paper_example4.pl

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.3 - (download) (as text) (annotate)
Wed Mar 31 18:23:02 2010 UTC (9 years, 7 months ago) by parrello
Branch: MAIN
CVS Tags: mgrast_dev_08112011, mgrast_dev_08022011, rast_rel_2014_0912, myrast_rel40, mgrast_dev_05262011, mgrast_dev_04082011, rast_rel_2010_0928, mgrast_version_3_2, mgrast_dev_12152011, mgrast_dev_06072011, rast_rel_2010_0526, rast_rel_2014_0729, mgrast_dev_02212011, rast_rel_2010_1206, mgrast_release_3_0, mgrast_dev_03252011, rast_rel_2011_0119, mgrast_release_3_0_4, mgrast_release_3_0_2, mgrast_release_3_0_3, mgrast_release_3_0_1, mgrast_dev_03312011, mgrast_release_3_1_2, mgrast_release_3_1_1, mgrast_release_3_1_0, mgrast_dev_04132011, mgrast_dev_04012011, rast_rel_2010_0827, myrast_33, rast_rel_2011_0928, mgrast_dev_04052011, mgrast_dev_02222011, mgrast_dev_10262011, HEAD
Changes since 1.2: +2 -3 lines
Fixed to use the correct servers.

#!/usr/bin/perl -w
use strict;
use SAPserver;
use COserver;
use SeedUtils;

my $sapObject = SAPserver->new();

my $genomeHash = $sapObject->all_genomes(-complete => 1);
for my $genome (keys %$genomeHash) {
    my $genomeName = $genomeHash->{$genome};
    my $geneHash = $sapObject->feature_assignments(-genome => $genome,
                                                   -type => 'peg');
    my @hypotheticalGenes = grep { &SeedUtils::hypo($geneHash->{$_}) } sort keys %$geneHash;
    my $couplingHash = $sapObject->conserved_in_neighborhood(-ids => \@hypotheticalGenes,
                                                             -hash => 1);
    for my $gene (@hypotheticalGenes) {
        my $couplingList = $couplingHash->{$gene};
        if (defined $couplingList) {
            my $subHash = $sapObject->ids_to_subsystems(-ids => [ map { $_->[1]} @$couplingList ]);
            my ($bestCoupler, $bestScore, $bestFunction) = (undef, 0, '');
            for my $coupling (@$couplingList) {
                my ($score, $coupler, $function) = @$coupling;
                if ($subHash->{$coupler} && $score > $bestScore) {
                    $bestCoupler = $coupler;
                    $bestScore = $score;
                    $bestFunction = $function;
                }
            }
            if (defined $bestCoupler) {
                print join("\t", $gene, $geneHash->{$gene}, $genome, $genomeName,
                                 $bestCoupler, $bestScore, $bestFunction) . "\n";
            }
        }
    }
}

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3