Hi,
This is my first day working with perl, so I really dont know anything. It took me whole day to come up with some code, but still it doesn't work :(
I have two text files which are tab delimited like this:
file 1 (test)
chr1 30
chr2 20
chr2 50
chr3 80
file 2 (reference)
chr1 40
chr1 50
chr2 60
chr2 80
chr3 100
I want to compare file1 which is the testfile agains file2 which is the reference file. When column 1 of both files is the same, I want to get the difference between the two files for column 2, and it must return the smallest value for that particular chr.
So the output should look like this:
output
chr1 10
chr2 40
chr2 10
chr3 20
My code looks like this now:
#!/usr/bin/env perl
use strict;
use warnings;
my ( @cols, $p1, $p2, $p3, $p4, @sec, @cols2 );
@ARGV or die "No input file specified";
open my $first , '<',$ARGV[0] or die "Unable to open input file: $!";
open my $second,'<', $ARGV[1] or die "Unable to open input file: $!";
print scalar <$first>;
<$second>; #...throw away first line...
while (<$first>) {
@cols = split /\s+/;
$p1 = $cols[0];
$p2 = $cols[1];
#print $p1;
while (<$second>) {
@cols2 = split /\s+/;
$p3 = $cols2[0];
$p4 = $cols2[1];
if ($p3 eq $p1) {
print "yes";
}
}
}
But this doesn't work.. Could somebody please help me?
Thanks!
Posts are HTML formatted. Put <p> </p> tags around your paragraphs. Put <code> </code> tags around your code and data!
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, 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, sup, table, tbody, td, tfoot, th, thead, tr, tt, u, ul, wbr
Outside of code tags, you may need to use entities for some characters:
| |
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.
|
|