From df8386f75b0538075d72d52693836bb8878f505b Mon Sep 17 00:00:00 2001 From: KatolaZ Date: Mon, 19 Oct 2015 16:23:00 +0100 Subject: First commit of MAMMULT code --- structure/metrics/avg_edge_overlap.py | 47 +++++++++++++++++++++++++++++++++++ 1 file changed, 47 insertions(+) create mode 100644 structure/metrics/avg_edge_overlap.py (limited to 'structure/metrics/avg_edge_overlap.py') diff --git a/structure/metrics/avg_edge_overlap.py b/structure/metrics/avg_edge_overlap.py new file mode 100644 index 0000000..e68fd8b --- /dev/null +++ b/structure/metrics/avg_edge_overlap.py @@ -0,0 +1,47 @@ +#### +## +## Compute the average edge overlap of a multiplex, i.e. the average +## number of layers in which an edge is present +## +## + +import sys + + +if len(sys.argv) < 2: + print "Usage: %s [...]" % sys.argv[0] + sys.exit(1) + +max_N = -1 + +all_edges = {} + +layer_ID = -1 + +for layer in sys.argv[1:]: + layer_ID += 1 + with open(layer, "r") as lines: + for l in lines: + if l[0] == "#": + continue + s, d = [int(x) for x in l.strip(" \n").split(" ")[:2]] + if s > d: + tmp = s + s = d + d = tmp + if (s,d) in all_edges: + all_edges[(s,d)].append(layer_ID) + else: + all_edges[(s,d)] = [layer_ID] + +K = len(all_edges.keys()) +M = len(sys.argv) - 1 + +numer = 0 + +for k in all_edges.keys(): + numer += len(set(all_edges[(k)])) + + +print 1.0 * numer / K, 1.0 * numer / (K * M) + -- cgit v1.2.3