The answers to
dave_the_m's questions will determine what sort of solution to look for. E.g. if the $data file is relatively small, you can load it into memory, then read each *.score.txt file exactly once to load up relevant info for each item in $data, then do stats on the info - something like this:
use strict;
use warnings;
my $data = "some/file.name";
open( F, $data ) or die "$data:$!\n";
my %targets;
while(<F>) {
next if ( /Header/ );
chomp;
my ( $chr, $start, $end ) = ( split( /\t/ ))[0,2,3];
push @{$targets{$chr}}, { table => $_, start => $start, end => $en
+d };
}
for my $chr ( keys %targets ) {
open( R, "$chr.score.txt" ) or die "$chr.score.txt: $!\n";
while (<R>) {
chomp;
my ( $pos, $score ) = ( split( /\t/ ))[1,2];
for my $range ( @{$targets{$chr}} ) {
if ( $pos >= $$range{start} and $pos <= $$range{end} ) {
push @{$$range{scores}}, $score;
}
}
}
}
# do statistics on contents of %targets…
If $data contains too much stuff, and/or requires too much stuff from the *.score.txt files to be held in memory for each $chr, then maybe you have to create output files for each $chr (or each start-end range in $data), so that it'll be quick/easy to do stats on those files.
Posts are HTML formatted. Put <p> </p> tags around your paragraphs. Put <code> </code> tags around your code and data!
Titles consisting of a single word are discouraged, and in most cases are disallowed outright.
Read Where should I post X? if you're not absolutely sure you're posting in the right place.
Please read these before you post! —
Posts may use any of the Perl Monks Approved HTML tags:
- a, abbr, b, big, blockquote, br, caption, center, col, colgroup, dd, del, details, div, dl, dt, em, font, h1, h2, h3, h4, h5, h6, hr, i, ins, li, ol, p, pre, readmore, small, span, spoiler, strike, strong, sub, summary, sup, table, tbody, td, tfoot, th, thead, tr, tt, u, ul, wbr
You may need to use entities for some characters, as follows. (Exception: Within code tags, you can put the characters literally.)
| |
For: |
|
Use: |
| & | | & |
| < | | < |
| > | | > |
| [ | | [ |
| ] | | ] |
Link using PerlMonks shortcuts! What shortcuts can I use for linking?
See Writeup Formatting Tips and other pages linked from there for more info.