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

View of /FigKernelScripts/generate_ilits.pl

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.3 - (download) (as text) (annotate)
Fri Sep 17 20:41:39 2010 UTC (9 years, 5 months ago) by olson
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_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, myrast_33, rast_rel_2011_0928, mgrast_dev_04052011, mgrast_dev_02222011, mgrast_dev_10262011, HEAD
Changes since 1.2: +4 -2 lines
fix for FFs changes

########################################################################

use FIG;
my $fig = new FIG;
my $rdbH  = $fig->db_handle;

use FFs;
use FF;

my $ff_data = $fig->get_figfams_data;

my $ffs = new FFs($ff_data);

#v usage: generate_ilits  > ilits";
#
#  ilits are [PEG-ilit,PEG-dlit,PubMed]


my $tuples = $rdbH->SQL("select DISTINCT md5_hash,pubmed
                                FROM dlits
                                WHERE status = 'D'");
foreach my $tuple (@$tuples)
{
    my($hash,$pubmed) = @$tuple;
    my @pegsD = $fig->pegs_with_md5($hash);
    if (@pegsD > 0)
    {
	my $pegD = $pegsD[0];
	my $fam;
	if ($fam = $ffs->family_containing_peg($pegD))
	{
#	    my $famO = new FF($fam,$ff_data);
#	    foreach my $pegI ($famO->list_members)

	    foreach my $pegI (@{$ffs->family_pegs($fam)})
	    {
		if ($pegI ne $pegD)
		{
		    print join("\t",($pegI,$pegD,$pubmed)),"\n";
		}
	    }
	}
    }
}



MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3