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

View of /FigKernelScripts/check_hope_info.pl

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.8 - (download) (as text) (annotate)
Mon Jan 25 21:15:20 2010 UTC (9 years, 9 months ago) by dejongh
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.7: +11 -1 lines
updates to check for invalid maps and reaction/ec mismatches

use strict;
use FIG;
use model;

my $fig = new FIG;
model::set_fig($fig);
my @subsystems = @ARGV;

if (scalar @subsystems == 0)
{
    my $ss_to_superset = model::load_superset_file;
    @subsystems = keys %$ss_to_superset;
}

foreach my $ss (sort @subsystems)
{
    print "Checking hope info for $ss\n";
    my ($invalid_roles, $invalid_reactions, $invalid_classification, $invalid_kegg_maps, $reaction_ec_mismatch) = model::check_hope_info($ss);
    print "\tInvalid Classification: $invalid_classification\n" if ($invalid_classification);

    foreach my $role (keys %$invalid_roles)
    {
	print "\tRole no longer exists: $role [@{$invalid_roles->{$role}}]\n";
    }

    foreach my $reaction (keys %$invalid_reactions)
    {
	print "\tInvalid reaction: $reaction [@{$invalid_reactions->{$reaction}}]\n";
    }

    foreach my $id (keys %$invalid_kegg_maps)
    {
	print "\tInvalid kegg map id: $id\n";
    }

    foreach my $reaction (keys %$reaction_ec_mismatch)
    {
	print "\tMismatched EC for reaction: $reaction [@{$reaction_ec_mismatch->{$reaction}}]\n";
    }
}


MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3