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

View of /FigKernelScripts/figfam_size_distribution.pl

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.1 - (download) (as text) (annotate)
Wed Jun 18 15:42:22 2008 UTC (11 years, 5 months ago) by arodri7
Branch: MAIN
CVS Tags: mgrast_dev_08112011, rast_rel_2009_05_18, mgrast_dev_08022011, rast_rel_2014_0912, myrast_rel40, mgrast_dev_05262011, rast_rel_2008_12_18, mgrast_dev_04082011, rast_rel_2008_07_21, rast_rel_2010_0928, rast_2008_0924, mgrast_version_3_2, mgrast_dev_12152011, mgrast_dev_06072011, rast_rel_2008_09_30, rast_rel_2009_0925, 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, mgrast_rel_2008_0924, mgrast_rel_2008_1110_v2, rast_rel_2009_02_05, rast_rel_2011_0119, mgrast_rel_2008_0625, 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, rast_rel_2008_10_09, mgrast_dev_04012011, rast_release_2008_09_29, mgrast_rel_2008_0806, mgrast_rel_2008_0923, mgrast_rel_2008_0919, rast_rel_2009_07_09, rast_rel_2010_0827, mgrast_rel_2008_1110, myrast_33, rast_rel_2011_0928, rast_rel_2008_09_29, mgrast_rel_2008_0917, rast_rel_2008_10_29, mgrast_dev_04052011, mgrast_dev_02222011, rast_rel_2009_03_26, mgrast_dev_10262011, rast_rel_2008_11_24, rast_rel_2008_08_07, HEAD
update changes

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

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

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

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

#foreach my $fam_id ($figfams->all_families){
#    my $figfam = new FigFam($fig,$fam_id,$famsD);
#
#    $distribution->{scalar @{$figfam->pegs_of()}}++;
#}
#

my $association;
open (FH, "$famsD/families.2c");
while (my $line = <FH>){
    chomp ($line);
    my ($fam_id, $peg) = split(/\t/, $line);
    $association->{$fam_id}++;
}
close FH;

foreach my $key (sort {$association->{$b} <=> $association->{$a}} keys %$association){
    $distribution->{$association->{$key}}++;
}

foreach my $key (sort {$a <=> $b} keys %$distribution){
    print $key . "\t" . $distribution->{$key} . "\n";
}

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3