def enrichment_no_t_task( exp, block, T, gs, patterns, base_filename, ): if settings.CELERY_DEBUG: import sys sys.path.append( '/Migration/skola/phd/projects/miXGENE/mixgene_project/wrappers/pycharm-debug.egg' ) import pydevd pydevd.settrace('localhost', port=6901, stdoutToServer=True, stderrToServer=True) gene_set = gs.get_gs() patterns = patterns.get_gs() e = EnrichmentInGeneSets(patterns.genes) enrich = e.getModuleEnrichmentInGeneSets(patterns.genes, gene_set.genes, pval_threshold=T) enrich = dict( (mod, (genes, map(lambda x: (gene_set.description[x[0]], x[0], x[1]), terms))) for (mod, (genes, terms)) in enrich.items()) ds = DictionarySet(exp.get_data_folder(), base_filename) ds.store_dict(enrich) return [ds], {}
def enrichment_no_t_task(exp, block, T, gs, cs, base_filename, ): if settings.CELERY_DEBUG: import sys sys.path.append('/Migration/skola/phd/projects/miXGENE/mixgene_project/wrappers/pycharm-debug.egg') import pydevd pydevd.settrace('localhost', port=6901, stdoutToServer=True, stderrToServer=True) gene_set = gs.get_gs() cs = cs.load_set() e = EnrichmentInGeneSets(cs) enrich = e.getModuleEnrichmentInGeneSets(cs, gene_set.genes, pval_threshold=T) enrich = dict((mod, (genes, map(lambda x: (gene_set.description[x[0]], x[0], x[1]), terms))) for (mod, (genes, terms)) in enrich.items()) ds = DictionarySet(exp.get_data_folder(), base_filename) ds.store_dict(enrich) return [ds], {}
def enrichment_task(exp, block, gs, H2, T, base_filename, ): if settings.CELERY_DEBUG: import sys sys.path.append('/Migration/skola/phd/projects/miXGENE/mixgene_project/wrappers/pycharm-debug.egg') import pydevd pydevd.settrace('localhost', port=6901, stdoutToServer=True, stderrToServer=True) gene_set = gs.get_gs() h2 = H2.get_assay_data_frame() e = EnrichmentInGeneSets(h2) ## compute enrichment in GO terms () enrich_bpGO = e.getEnrichmentInGeneSetsWithH(gene_set.genes, h2, T) # sort resultst accodring p-values # sorted_enrich_bpGO = sorted(enrich_bpGO.iteritems(), key=operator.itemgetter(1)) er_ratio = e.getEnrichmentRatioInGeneSetsWithH(gene_set.genes, h2, T, enrichment_threshold=0.05, N=10) enrich_bpGO['er_ratio'] = er_ratio ds = DictionarySet(exp.get_data_folder(), base_filename) ds.store_dict(enrich_bpGO) return [ds], {}