#! /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']))
#! /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']))
