#! /usr/bin/env python
# -*- coding: utf-8 -*-

######################################################
# map the gene (ENSEMBL_ID) to its respective color  #
# and count the coverage of each GO_id in each color #
# output (tab-delimited file):                       #
# go_id  Red   Pink  Yellow   Black                  #
######################################################

import sys

# parsed dictionary from human_go_GF20130505.json
# d_go = d_go = {"GO:0022607": ["ENSG00000102054", "ENSG00000172426",.],.} 
from d_go import d_go

# dictionary containing the mapping
# d_uniqueExpressedColorMajority ={'ENSG00000227999': 'Black', ...}
from d_uniqueExpressedColorMajority import d_uniqueExpressedColorMajority

for go_id in d_go.iterkeys():
   # map the gene_id to its respective color
   color_list = [d_uniqueExpressedColorMajority[i] if i in d_uniqueExpressedColorMajority else "NA" for i in d_go[go_id]]
   # count the coverage of each go_id in each color
   color_count = dict(zip(color_list, map(color_list.count, color_list)))
   sys.stdout.write("%s\t%d\t%d\t%d\t%d\n"%(go_id, color_count['Red'],
                                             color_count['Pink'], 
                                             color_count['Yellow'],
                                             color_count['Black']))
예제 #2
0
#! /usr/bin/env python
# -*- coding: utf-8 -*-

######################################################
# map the gene (ENSEMBL_ID) to its respective color  #
# and count the coverage of each GO_id in each color #
# output (tab-delimited file):                       #
# go_id  Red   Pink  Yellow   Black                  #
######################################################

import sys

# parsed dictionary from human_go_GF20130505.json
# d_go = d_go = {"GO:0022607": ["ENSG00000102054", "ENSG00000172426",.],.} 
from d_go import d_go

# dictionary containing the mapping
# d_color ={'ENSG00000227999': 'Black', ...}
from d_color import d_color

for go_id, genes_id in d_go.iterkeys():
   # map the gene_id to its respective color
   color_list = [d_color[i] if i in d_color else "NA" for i in genes_id]
   # count the coverage of each go_id in each color
   color_count = dict(zip(color_list, map(color_list.count, color_list)))
   sys.stdout.write("%s\t%d\t%d\t%d\t%d\n"%(go_id, color_count['Red'],
                                             color_count['Pink'], 
                                             color_count['Yellow'],
                                             color_count['Black']))