Dear Cristoforo, When I ran the script I realized that I am too stringent about my selection criteria.

In the second file if $current_line3 >= 3 then the I need to fake around plus/minus 8 of the corresponding position to this line which is ($current_line1). So the condition would be find from file 2 those that $current_line3 >= 3 and it is OK to that ($current_line1)plus minus 8 matches with file 1.

I thought if I put($current_line1)+/- 8 in an array. then say if from file1 $position1 exists in this array it is OK.

If I am not clear please let me know to explain it more. 


while(<INPUT2>){
         chomp;
          my @current_line  = split /\t/;
          # eliminate unqualified lines early
          next unless $current_line[2] == 1 && $current_line[3] >= 3;
          #$from = $current_line[1]-8;
          #$to = $current_line[1]+8;
              # for ($from .. $to){
              # push (@range, $_);
              #      }
              #}
        if ($file1{ $current_line[0]}) {
                ($from, $to) = @{ $file1{ $current_line[0] } };
                if ($from <= $current_line[1] && $current_line[1] <= $
+to) {

                print join("\t", $_, "***",$current_line[1]), "\n";


           $true_positives++;
            }
      }
}

[download]

In reply to Re^2: Reading two files, cmp certain cols by sesemin
in thread Reading two files, cmp certain cols by sesemin

Title:
Use:  <p> text here (a paragraph) </p>
and:  <code> code here </code>
to format your post, it's "PerlMonks-approved HTML":


  • 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:
    & &amp;
    < &lt;
    > &gt;
    [ &#91;
    ] &#93;
  • Link using PerlMonks shortcuts! What shortcuts can I use for linking?
  • See Writeup Formatting Tips and other pages linked from there for more info.