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

View of /FigKernelScripts/FFB2_xy.pl

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.1 - (download) (as text) (annotate)
Fri Oct 29 17:10:34 2010 UTC (9 years 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, 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
Update to FF processing

#
# run usable_motifs on a set of per-base kmer possibilites
#

use strict;
use Proc::ParallelLoop;

@ARGV == 3 or die "Usage: $0 input-dir output-dir range\n";

my $input_dir = shift;
my $output_dir = shift;
my $range = shift;

my @input = <$input_dir/all.*>;
my %input;

for my $f (@input)
{
    if ($f =~ /all.([A-Z])$/)
    {
	$input{$1} = $f;
    }
}

my @c = sort { $a cmp $b } keys %input;

pareach \@c, sub {
    my($c) = shift;
    my $file = $input{$c};
    print "$$ $c $file\n";
    for my $col (2, 3)
    {
	my $cmd = "FFB2_usable_motifs2 $output_dir/kmers.$col.$c $range $col $file";
	my $rc = system($cmd);
	if ($rc != 0)
	{
	    print "error $rc: $cmd\n";
	}
    }
}, { Max_Workers => 8 }

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3