From 97fdefab064f63642fa3ece05b807d29b459df31 Mon Sep 17 00:00:00 2001 From: Hui Lan Date: Wed, 4 Dec 2019 19:03:19 +0800 Subject: brain: add python and R code to local repository. --- Code/brain_number_genes_edges.py | 48 ++++++++++++++++++++++++++++++++++++++++ 1 file changed, 48 insertions(+) create mode 100644 Code/brain_number_genes_edges.py (limited to 'Code/brain_number_genes_edges.py') diff --git a/Code/brain_number_genes_edges.py b/Code/brain_number_genes_edges.py new file mode 100644 index 0000000..d1ace68 --- /dev/null +++ b/Code/brain_number_genes_edges.py @@ -0,0 +1,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)) -- cgit v1.2.1