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

View of /FigKernelScripts/figfam_subsystems.pl

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.2 - (download) (as text) (annotate)
Thu Oct 8 18:54:01 2009 UTC (10 years, 1 month ago) by arodri7
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_2010_0118, 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.1: +21 -3 lines
commit changes

#!/usr/bin/env /home/arodri7/FIGdisk/bin/run_perl

use Data::Dumper;
use Carp;
use FIG_Config;
use FIG;
use FF;

my $usage = "figfam_subsystems $famsDirectory\n";
my $famsD = $ARGV[0];

use FigFams;
my $fig = new FIG;
my $figfams = new FigFams($fig,$famsD);

open (FH, "$famsD/family.functions");
while (my $line = <FH>){
    next if ($line =~ /not subsystem-based/);
    chomp ($line);
    my ($fam_id, $function) = split(/\t/, $line);
    my $ff_obj = new FF($fam_id,$famsD);
    my $subs={};
    foreach my $id (@{$ff_obj->pegs_of})
    {
	my @subs1 = $fig->subsystems_for_peg($id);
	my $count = scalar (@{$subs1[0]});
	my $i=0;
	foreach my $sub (@{$subs1[0]})
	{
	    last if ($count-1 == $i);
	    $subs->{$sub}++;
	    $i++;
	}
    }
    my @subsystems = keys %$subs;
    if (scalar @subsystems > 0)
    {
	#print join ("\t", ($fam_id, $function, $fig->function_to_subsystems($function))) . "\n" if (! $fig->function_to_subsystems($function));
	print join ("\t", ($fam_id, $function, @subsystems)) . "\n";
    }
}
close FH;


MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3