1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
|
#!/usr/bin/perl
use strict;
#
# Front-end for making intersection of more than two list files
# Copyright (C) Maarja Lepamets and University of Tartu
# You can modify and distribute this script and derivateive works freely
#
# Usage:
# MakeIntersection.pl LIST1 LIST2 LIST3...
#
# NB! Please set the path to executables
#
use File::Basename;
my $nlists = scalar @ARGV;
my $list1;
my $list2;
my $out = "intrs";
my $xpath = ".";
if (!(-x "$xpath/glistcompare")) {
printf stderr "No glistcompare binary at %s\n", $xpath;
printf stderr "Please set the xpath variable in script code\n";
exit(1);
}
my $k = 1;
while ($nlists != 1) {
if ($nlists == scalar @ARGV) {
print STDERR "mkdir -p ${out}_${k}\n";
system "mkdir -p ${out}_${k}";
my $i = 0;
while ($i < $nlists) {
$list1 = $ARGV[$i];
if ($i == $nlists - 1) {
print STDERR "cp $list1 ${out}_${k}/copy_".basename($list1) . "\n";
system "cp $list1 ${out}_${k}/copy_".basename($list1);
last;
}
$list2 = $ARGV[$i + 1];
print STDERR "$xpath/glistcompare $list1 $list2 -o ${out}_${k}/$i"."_".($i + 1)." -i\n";
system "$xpath/glistcompare $list1 $list2 -o ${out}_${k}/$i"."_".($i + 1)." -i";
$i += 2;
}
} else {
#system "rm -r -f ${out}_".($k - 2);
my $i = 0;
my $dir = "${out}_".($k - 1);
my @files = <$dir/*>;
$nlists = scalar @files;
if ($nlists == 2) {
$list1 = $files[$i];
$list2 = $files[$i + 1];
print STDERR "$xpath/glistcompare $list1 $list2 -o intrs -i\n";
system "$xpath/glistcompare $list1 $list2 -o intrs -i";
#system "rm -r -f ${out}_".($k - 1);
last;
}
print STDERR "mkdir -p ${out}_${k}\n";
system "mkdir -p ${out}_${k}";
while ($i < $nlists) {
$list1 = $files[$i];
if ($i == $nlists - 1) {
print STDERR "cp $list1 ${out}_${k}/copy_".basename($list1) . "\n";
system "cp $list1 ${out}_${k}/copy_".basename($list1);
last;
}
$list2 = $files[$i + 1];
print STDERR "$xpath/glistcompare $list1 $list2 -o ${out}_${k}/$i"."_".($i + 1)." -i\n";
system "$xpath/glistcompare $list1 $list2 -o ${out}_${k}/$i"."_".($i + 1)." -i";
$i += 2;
}
}
$nlists = int($nlists / 2 + 0.5);
$k += 1;
}
|
