I want to first submit all of the data to the cluster and then do a 30+% and 30-% selection of the data on the bootstrapped set. How could I implement this?#!/usr/bin/perl -w ###################################################################### +########### # + # # Program to submit bootstrap jobs to cluster + # # + # # Usage: qsubbootstrap.pl <output dir> <lower#> <upper#> <ARACNE argum +ents> # # Example: qsubbootstrap.pl adjdir 1 100 -i data.exp -l tflist.txt -s +probe.txt # # -p 1e-5 -e 0 + # # + # # Description: This program submit the bootstrap jobs to cluster and + # # store all the output bootstrapping networks in the directory <outp +ut dir>. # # The directory will be created if it doesn't exist. + # # Users also need to specify the range of bootstrapping samples to s +ubmit. # # + # ###################################################################### +########### use Fcntl; if (scalar(@ARGV) < 4) { die "Incorrect number of arguments!\nUsage: qsubbootstrap.pl <output + dir> <lower#> <upper#> <ARAC arguments> \nExample: qsubbootstrap.pl +adjdir 1 100 -i data.exp -l tflist.txt -p 1e-5 -e 0\n"; } else { ($dir,$r1,$r2,@otherarg)=@ARGV; $otherarg="@otherarg"; } $otherarg =~ m/-i\s+(\S+)\.\S+\s/; $inputfile = $1; $dir=$dir."\/" unless ($dir=~m/\/$/); @folderfound=<./*/>; $allfolder="@folderfound"; system "mkdir ".$dir unless ($allfolder=~m/$dir/); foreach $bs ($r1..$r2){ $bs3 = sprintf("%03i", $bs); $shfile = "aracne".$bs3."\.sh"; sysopen(OUTFILE, $shfile, O_WRONLY|O_TRUNC|O_CREAT, 0666) || die $!; select OUTFILE; print "\#\$ -S \/bin\/bash -cwd -j y -o soe \n\n \.\/arac $otherarg +-r $bs -o ".$dir.$inputfile."_r".$bs3.".adj"; close OUTFILE; system "qsub -p -100 $shfile"; system "rm ".$shfile; }
In reply to Bootstrapping Implementation #2 by hoffmann
| For: | Use: | ||
| & | & | ||
| < | < | ||
| > | > | ||
| [ | [ | ||
| ] | ] |