summaryrefslogtreecommitdiff
path: root/Code/brain_number_genes_edges.py
blob: d1ace68b2ff9e771da6a78a725db3afd84f96f52 (plain)
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
import os, operator, sys

def make_dict(fname):
    f = open(fname)
    gene_dict = {}
    target_dict = {}
    tf_dict = {}
    edge_dict = {}
    for line in f:
        line = line.strip()
        lst = line.split('\t')
        target = lst[0].split()[0].upper()
        tf     = lst[1].split()[0].upper()
        if not tf in tf_dict:
            tf_dict[tf] = 1
        else:
            tf_dict[tf] += 1
            
        if not target in target_dict:
            target_dict[target] = 1
        else:
            target_dict[target] += 1

        if not tf in gene_dict:
            gene_dict[tf] = 1
        else:
            gene_dict[tf] += 1
        if not target in gene_dict:
            gene_dict[target] = 1
        else:
            gene_dict[target] += 1
        k = tf + '.' + target
        if not k in edge_dict:
            edge_dict[k] = 1
        else:
            edge_dict[k] += 1
            
    f.close()
    return tf_dict, target_dict, gene_dict, edge_dict


# main
edge_fname = '/home/hui/network/v03/Data/history/edges/edges.txt'
tf_dict, target_dict, gene_dict, edge_dict = make_dict(edge_fname)
print('Number of TFs: %d' % len(tf_dict))
print('Number of targets: %d' % len(target_dict))
print('Number of genes: %d' % len(gene_dict))
print('Number of edges: %d' % len(edge_dict))