This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
genetica:bonngwas_all [2013/01/31 13:42] tgif |
genetica:bonngwas_all [2020/08/04 10:58] (current) |
||
---|---|---|---|
Line 48: | Line 48: | ||
my $ifile = shift; | my $ifile = shift; | ||
- | |||
die "Must supply input filename!\n" | die "Must supply input filename!\n" | ||
- | |||
(my $bn) = $ifile =~ / | (my $bn) = $ifile =~ / | ||
my $pfile = $bn.' | my $pfile = $bn.' | ||
my $ofile = $bn.' | my $ofile = $bn.' | ||
- | |||
my $tl = int qx/wc -l $ifile | awk {' | my $tl = int qx/wc -l $ifile | awk {' | ||
- | |||
- | |||
open(IPF, "< | open(IPF, "< | ||
- | |||
- | |||
my @pairs; | my @pairs; | ||
- | |||
print " | print " | ||
- | |||
my $l = 0; | my $l = 0; | ||
- | |||
while(< | while(< | ||
(my $a, my $b) = / | (my $a, my $b) = / | ||
Line 76: | Line 66: | ||
my @aps = uniq get_proxies($a, | my @aps = uniq get_proxies($a, | ||
my @bps = uniq get_proxies($b, | my @bps = uniq get_proxies($b, | ||
- | |||
foreach my $ap (sort @aps){ | foreach my $ap (sort @aps){ | ||
foreach my $bp (sort @bps){ | foreach my $bp (sort @bps){ |